t = [0:sqrt(25*pi)/1000:sqrt(25*pi)];
x= cos(t.^2);
y=sin(t.^2);
figure;
comet(x,y);
pause

s = [0:25*pi/1000:25*pi];
x= cos(s);
y=sin(s);
figure;
comet(x,y);