
clc;   %zadanie od fizyków
clear;
ymin=-10;
ymax=10;
krok=0.001;
ilosc=(ymax-ymin)/krok+1;
Y=[ymin:krok:ymax];
s=0.0448;
r=2;
G=r*(0.491-0.1018*r-0.003678*r*r);
n=size(Y,2);
EX=[1:1:ilosc];
x=8.08;
for i=1:ilosc
    y=Y(i);
    EX(i)=1/(exp(y/s));
end;

xmin=6;
xmax=80;
xkrok=5;
xilosc=(xmax-xmin)/xkrok+1;
X=[xmin:xkrok:xmax];
m=size(X,2);

for j=1:xilosc
    W=[1:1:ilosc];
    x=X(j);
    for i=1:ilosc
        y=Y(i);
        W(i)=(((y*y*(2-y)*(1+(G/(s*s*x*x)))/2)-(1-y*y*y))/s)-1;
    end;
    plot(Y,W);
    axis([-5 10 -40 300]);
    hold on;
    M(j)=getframe;
end;
hold on;
movie(M);
hold on;
plot(Y,EX);
