x(u,v) = sin (3 * v) * cos (4 * u) * (v * (sin(2 * u))^2 + sin (v)),
y(u,v) = sin (3 * v) * cos (4 * u) * (v * (sin(2 * u))^2 + sin (v)),
z(u,v) = cos (3 * v) * (v * (sin(2 * u))^2 + sin (v)).
with (plots):
setoptions3d(style = patchnogrid,
scaling = constrained,
projection = .7,
ambientlight = [ .3, .3, .3 ],
light = [65, 50, .9, .9, .9]);
f := (x, y) -> [ sin(3 * y) * cos(4 * x) * (sin(2 * x)^2 * y + sin(y)),
sin(3 * y) * sin(4 * x) * (sin(2 * x)^2 * y + sin(y)),
cos(3 * y) * (sin(2 * x) * (sin(2 * x)^2 * y + sin(y)))];
d := [Pi/4, Pi/2, 0, Pi * .45, 40, 40, 40]:
plot3d(f(x,y), x = d[1] .. d[2], y = d[3] .. d[4], grid = [d[5], d[6]],
color = [1, .5, .4], orientation = [-20, 50]);