%
% --------------- rose2 ------------------
%
function r = rose2(x,y);
rho = sqrt(x.^2 + y.^2);
theta = atan2(y,x);
r = rho - .5*abs(sin(5*theta));
r = (rho > 1e-8).*r;