function f=f(y);
    f=0.08*y
    //*(1-y/400000);
endfunction;
f(60000);
h=0.1;
function eul=Euler(y);
    eul=y+h*f(y);
endfunction;
y0=60000;
f(y0);
y1=Euler(y0)
f(y1);
y2=Euler(y1)

//introduction du système de 2 équation (t'=1 et y'=-4y)
clear;
function f=un(t,y);f=1;endfunction;
function g=g(t,y);g=-4*y;endfunction;
function w=w(t,v);w=[un(v(1),v(2));g(v(1),v(2))];endfunction;
//tracé du champ de vecteurs associé
ttChamp=0:0.25:4;yyChamp=-5:0.5:5;
xset("window",1);fchamp(w,0,ttChamp,yyChamp);
//tracé de la solution de condition initiale 0 (sol particulière) et de celle 
//de condition initiale 5
tt=0:0.1:4;
yy0=ode(0,0,tt,g);
yy5=ode(5,0,tt,g);
plot(tt,yy0);plot(tt,yy5);



clear;
//introduction du système de 2 équation (t'=1 et y'=-4y+3cos(t))
function f=un(t,y);f=1;endfunction;
function g=g(t,y);g=-4*y+3*cos(t);endfunction;
function w=w(t,v);w=[un(v(1),v(2));g(v(1),v(2))];endfunction;
//tracé du champ de vecteurs associé
ttChamp=0:0.5:15;yyChamp=-5:0.5:5;
xset("window",1);fchamp(w,0,ttChamp,yyChamp);
//tracé de la solution de condition initiale 12/17 (sol particulière) en vert
tt=ttChamp;
yy=ode(12/17,0,tt,g)
plot(tt,yy,color='g')

//tracé de la même solution par la méthode d'Euler
pas=0.5;
ttEuler(1)=0;yyEuler(1)=12/17;
for i=1:pas:15/pas
    ttEuler(i+1)=ttEuler(i)+pas;
    yyEuler(i+1)=yyEuler(i)+pas*g(ttEuler(i),yyEuler(i));
end;
plot(ttEuler,yyEuler,color='r')

