Zad 1. clear; set(gcf,'NumberTitle','Off','MenuBar','None','Name','x^k'); W=2:2:20; // b) zmieniamy na W=1:2:21 [m,n]=size(W); x=-5:0.01:5; for k = 1:n y=x.^(W(k)); plot(x,y); axis equal M(k)=getframe; end movie(M,10,10) c) clear; set(gcf,'NumberTitle','Off','MenuBar','None','Name','x^k'); W=0:0.5:1 [m,n]=size(W); x=-2*pi:0.01:2*pi; for k = 1:n y=sin(x+W(k)*2*pi); plot(x,y); axis equal M(k)=getframe; end movie(M,10,10) 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)