I found this here. I had to do this in Euler, just to see. Here is the command

>plot2d("sqrt(cos(x))*cos(200*x)+sqrt(abs(x))-0.7*(4-x^2)^0.01", ...
> -2,2,n=10000,frame=0,grid=0,color=red);


and this is the result:

How is it working? It is made of three nice functions.


>plot2d("sqrt(cos(x))+sqrt(abs(x))-0.7*(4 -x^2)^0.01",r=1.8,color=red);
>plot2d("-sqrt(cos(x))+sqrt(abs(x))-0.7*(4 -x^2)^0.01",add=true,color=blue);
>plot2d("sqrt(abs(x))-0.7*(4 -x^2)^0.01",add=true,color=black);
>


Which can be seen in the following plot.

The cos(200x) term just oscillates between the red and blue bounds. In Euler, the functions becomes invalid for |x| too large and stops plotting.

Nice!