%% skrypt do calkowania polowu ryb -ludwiga
%% wymaga pliku polowy2.m
global r y
r=2;
y=r/4-0.006;
%% delta =0.008
dx=0.1;
%% epsilon = 0.1
n0=1/2+sqrt(1/4-y/r);
%% runge-kutty metoda calkowania funkcji polowy2 od t =0 do t=20 "po dodatniej stronie":
[tp,yp]=ode45(@polowy2,[0,20],n0+dx);
%% runge-kutty metoda calkowania funkcji polowy2 od t =0 do t=20 "po ujemnej stronie":
[tn,yn]=ode45(@polowy2,[0,20],n0-dx);
plot(tp,yp,tn,yn,[0,20],[n0, n0], [0,20],[1-n0,1-n0]);