1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53
|
clear;
funcprot(0);
xdel(winsid()); //ferme les fenêtres graphiques
function [out]=iso_1(D,S_in,s)
for i=1:length(s)
out(i)=(D*(S_in-s(i)))/(s(i)/(0.2+s(i)))
end
endfunction
//tracer des isoclines
scf(3)
t2=linspace(0.2,2,100);
plot2d(t2,iso_1(0.8,2,t2),2,rect=[0,0,2.5,2])
plot2d(0.8*ones(1,100),linspace(0,2,100),3)
legends(['$\frac{ds}{dt}=0$';'$\frac{dx}{dt}=0$'],[2,3],opt="ur")
plot(0.8,1.2,'ro','markerfacecolor','r')
xstring(0.8,1.2,'$E_1$')
xstring(0.8,0,'$\lambda(D)$')
plot(2,0,'ro','markerfacecolor','r')
xstring(2,0,'$E_0$')
xtitle('Isoclines')
xlabel("s");
ylabel("x");
xrect(0,2,2.5,2) //pour avoir le cadre
x1=[1.4,0.8,0.8,0.5,0,0]
y1=[0.5485714,0.5,1.7,1.68,0.3,1.7];
x2=[1.4,0.9,0.7,0.5,0.1,0.1]
y2=[0.5485714+0.1,0.5,1.7,1.58,0.28,1.62];
X=[x1;x2];Y=[y1;y2]
xarrows(X,Y,3,5)
e=gce()
e.thickness=1
e.mark_size=0.1
e.arrow_size=1
e.mark_mode="on"
x11=[0,2.4]
y11=[0,0];
x22=[2,2]
y22=[0,0];
X=[x11;x22];Y=[y11;y22]
xarrows(X,Y,5)
e=gce()
e.thickness=2
e.mark_size=0.1
e.arrow_size=0.5
e.mark_mode="on"
xs2png(3, 'isoclines') |
Partager