1. Zad2.
2. a=50;
3. [X,Y] = meshgrid(-a:.5:a);
4. R = sqrt(X.^2 + Y.^2) + eps;
5. figure('Resize','off');
6. W=[1:-0.01:0 0.01:0.01:1];
7. [m,n]=size(W);
8. for k = 1:n
9.  Z = sin(R*W(k))./(R*W(k));
10.  surf(X,Y,Z)
11.  axis([-a a -a a -0.2 1]);
12.  M(k) = getframe;
13. end
14. movie(M,10,30)
15.  
16. zad3.
17. a=20;
18. [X,Y] = meshgrid(-a:0.01:a);
19. R = figure('Resize','off','units','normalized','outerposition',[0 0 1 1])
20. W=[1:-0.01:0 0.01:0.01:1];
21. [m,n]=size(W);
22. for k = 1:n
23.  Z = cos(W(k)*x)+(W(k)*y);
24.  mesh(X,Y,Z)
25.  axis([-a a -a a -0.2 1]);
26.  M(k) = getframe;
27. end
28. movie(M,10,30)
29.  
30. Zadanie 5
31.  
32. W=[1:-0.01:0 0.01:0.01:1];
33. [m,n]=size(W);
34. for k = 1:n
35. t=0:0.001:2*pi;
36.  
37.  
38. x = cos(t)+cos(7*t)/2+sin(17*t)/3;
39. y = sin(t)+sin(7*t)/2+cos(17*t)/3;
40. comet(x,y);
41.  axis equal
42.  M(k) = getframe;
43. end
44. movie(M,10,30)
45.  
46. Zadanie 6
47.  
48. W=[1:-0.01:0 0.01:0.01:1];
49. [m,n]=size(W);
50. for k = 1:n
51. t=0:0.001:20*pi;
52. a=4;
53. b=9
54. x = a*t-b*sin(t);
55. y = a-b*cos(t);
56.  
57. comet(x,y);
58.  axis equal
59.  M(k) = getframe;
60. end
61. movie(M,10,30)
62.  
63.  
64. Zadanie 7
65.  
66. W=[1:-0.01:0 0.01:0.01:1];
67. [m,n]=size(W);
68. for k = 1:n
69. t=0:0.001:10*pi;
70. a=0.1;
71. x = a*t*cos(t);
72. y = a*t*sin(t);
73.  
74. comet(x,y);
75.  axis equal
76.  M(k) = getframe;
77. end
78. movie(M,10,30)
79.  
80.  
81. Zadanie 8
82.  
83. W=[1:-0.01:0 0.01:0.01:1];
84. [m,n]=size(W);
85. for k = 1:n
86. t=0:0.001:10*pi;
87. a=10;
88. b=3;
89. c=6;
90. d=4;
91. A=20;
92. B=10;
93. x = A*cos(a*t+b);
94. y = B*cos(a*t+d);
95. z=b*t;
96. comet3(x,y,z);
97.  axis equal
98.  M(k) = getframe;
99. end
100. movie(M,10,30)