clear all;
rad = pi/180;
a = 1*rad;
b = 1*rad;
omega = 10;

t = 0:pi/50:30*pi;
%x =  a*omega*cos(omega*rad*t);
%y =  b*omega*sin(omega*rad*t);
%z = -a*b*omega*(1-cos(2*omega*rad*t));
x = a*sin(omega*rad*t) ;
y = -b*cos(omega*rad*t);
z = a*b*omega*sin(2*omega*rad*t)/4.0 - a*b*omega*t/2.0;

%plot3(x,y,z)
plot(t,z);