Zad2.
a=50;
[X,Y] = meshgrid(-a:.5:a);
R = sqrt(X.^2 + Y.^2) + eps;
figure('Resize','off');
W=[1:-0.01:0 0.01:0.01:1];
[m,n]=size(W);
for k = 1:n
 Z = sin(R*W(k))./(R*W(k));
 surf(X,Y,Z)
 axis([-a a -a a -0.2 1]);
 M(k) = getframe;
end
movie(M,10,30)

zad3.
a=20;
[X,Y] = meshgrid(-a:0.01:a);
R = figure('Resize','off','units','normalized','outerposition',[0 0 1 1])
W=[1:-0.01:0 0.01:0.01:1];
[m,n]=size(W);
for k = 1:n
 Z = cos(W(k)*x)+(W(k)*y);
 mesh(X,Y,Z)
 axis([-a a -a a -0.2 1]);
 M(k) = getframe;
end
movie(M,10,30)

Zadanie 5

W=[1:-0.01:0 0.01:0.01:1];
[m,n]=size(W);
for k = 1:n
t=0:0.001:2*pi;


x = cos(t)+cos(7*t)/2+sin(17*t)/3;
y = sin(t)+sin(7*t)/2+cos(17*t)/3;
comet(x,y);
 axis equal
 M(k) = getframe;
end
movie(M,10,30)

Zadanie 6

W=[1:-0.01:0 0.01:0.01:1];
[m,n]=size(W);
for k = 1:n
t=0:0.001:20*pi;
a=4;
b=9
x = a*t-b*sin(t);
y = a-b*cos(t);

comet(x,y);
 axis equal
 M(k) = getframe;
end
movie(M,10,30)


Zadanie 7

W=[1:-0.01:0 0.01:0.01:1];
[m,n]=size(W);
for k = 1:n
t=0:0.001:10*pi;
a=0.1;
x = a*t*cos(t);
y = a*t*sin(t);

comet(x,y);
 axis equal
 M(k) = getframe;
end
movie(M,10,30)


Zadanie 8

W=[1:-0.01:0 0.01:0.01:1];
[m,n]=size(W);
for k = 1:n
t=0:0.001:10*pi;
a=10;
b=3;
c=6;
d=4;
A=20;
B=10;
x = A*cos(a*t+b);
y = B*cos(a*t+d);
z=b*t;
comet3(x,y,z);
 axis equal
 M(k) = getframe;
end
movie(M,10,30)