u/Future_Fee_9765

On the CD booklet for Huge Chrome Cylinder Box Unfolding, there's some lines of code that seem to outline a sort of 3d plot algorithm, and I was wondering if there was a way to run the code?

Some of the syntax is a little different than what I'm used to seeing, if anyone has any ideas let me know. Transcribed to the best of my ability.

x(u.v)=sin(3v)cos(4u)(vsin^2(2u)+sinv). 
y(u.v)=sin(3v)cos(4u)(vsin^2(2u)+sinv).
z(u.v)=cos(3v)(vsin^2(2u)+sinv).

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]):