%% skrypt do calkowania lotki-voltery
%% wymaga pliku polowy3.m
global alfa
alfa=10;

%% delta =0.008

%% epsilon = 0.1
clf
hold on
for i=1:10
    x=rand();
    y=rand();
%% runge-kutty metoda calkowania funkcji polowy3 od t =0 do t=20 
options = odeset('RelTol',1e-5,'AbsTol',1e-8);%MaxStep',1e-3);
[t,g]=ode45(@polowy3,[0,200],[1+x 1+y],options);
ind=find(t>150);
plot(g(ind,1),g(ind,2),'color',[1/i 1/i 0]);
end
hold off;