Rotación de funciones en 2D respecto al eje Z (script en Matlab)

Código de Matlab:

% rotate 2D in z-axis

x = -5:0.1:5;

y = sin(3*x);

v = [x(:) y(:)];

for theta = 0:0.01:pi
 
   R = [cos(theta) sin(theta); -sin(theta) cos(theta)];
 
   vr = v*R;
 
   x = vr(:,1);
   y = vr(:,2);
 
   plot(x,y);
 
   axis([-8 8 -8 8]);
 
   pause(0.001);
 
end 

Video: https://youtu.be/RMj5En2riuk