%{
Zadanie 3 z listy 4. Animacja zmian sfery jednostkowej S_p, dla p zmieniaj¹cego
siê od 1 do nieskonczonosci. n=2.
%}

X=zeros(161,1);
X_p=X;
Y=zeros(161,1);
Y_p=Y;
for i=1:161                 %149 to ilosc punktow na okregu jednostkowym
  X(i)=cos(2*pi*(i-1)/160);
  Y(i)=sin(2*pi*(i-1)/160);
end;
for i=1:200                 %200 to ilosc klatek filmu
  p=1/(1-(i-1)/200);        %p przebiegnie zakres 1 do 200 w 200 krokach
  for j=1:161
    norma=(abs(X(j))^p+abs(Y(j))^p)^(1/p);
    X_p(j)=X(j)/norma;
    Y_p(j)=Y(j)/norma;
  end;
  plot(X_p,Y_p);
  axis([-2 2 -2 2]);
  axis square;
  M(i)=getframe;
end
for j=1:161                 %ostatnia klatka animacji: p=nieskoñczonoœæ
  if abs(X(j))>abs(Y(j))
    norma=abs(X(j));
  else
    norma=abs(Y(j));
  end;
  X_p(j)=X(j)/norma;
  Y_p(j)=Y(j)/norma;
end;
plot(X_p,Y_p);
axis([-2 2 -2 2]);
axis square;
M(201)=getframe;
movie(M,1);
