% Uczymy siec zachowywac sie jak oscylator van der Pola
% w oscylatorze tym mamy dwie zmienne y(1) i y(2)
% jest on opisany rownaniem :
% dy(1)/dt = y(2)
% dy(2)/dt = (1-y(1)^2)*y(2)-y(1);
% czyli dla matlaba funkcja vdp 
% function dydt = vdp(t,y)
% dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];
%Robimy lekcje
%dla wielu roznych stanow poczatkowych obliczamy stany
%za 0.05 jedn. czasu
% najpierw wytworzymy troche syganlu i go obejrzymy:
 [t,y] = ode45(@vdp,[0 5],[1; 1]);
 plot(t,y)
 drawnow
 clear all
i=0;
for y1_0= -2:0.05:2
    for y2_0= -3:0.05:3
        i=i+1;
       [t,y] = ode45(@vdp,[0 0.05],[y1_0; y2_0]);
       %wejsciem do sieci jest stan poczatkowy ukladu
       X(i,:)=[y1_0 y2_0];
       % a pozadanym wyjsciem stan w chwili o 0.05 jedn. czasu pozniejszej
       Z(i,:)=[y(end,1)-y1_0 y(end,2)-y2_0];
   end
end
save  lekcje X Z




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%