getf("sys_diff.sci")
//
// Exemple 1 : exponentielle
//
xset("window",0);xbasc();
// definition de la fonction 
deff("yprim=f1(t,y)","yprim=y");
Y0 = [-2,-1,1,2]; tmax = 3; pas = 0.01;
champty(f1,Y0,tmax,pas);
//
// Exemple 2 : Verhulst
//
xset("window",1);xbasc();
// definition de la fonction
deff("yprim=f2(t,y)","yprim=2*y*(1-y/4)");
Y0 = [0.5,1,1.5,2]; tmax = 4; pas = 0.01;
champty(f2,Y0,tmax,pas);
//
// Exemple 3
//
xset("window",2);xbasc();
// definition de la fonction
deff("yprim=f3(t,y)","yprim=2*y-y^2");
Y0 = [-0.01,0.5.1.5,4]; tmax = 2; pas = 0.01;
champty(f3,Y0,tmax,pas);
//
// Exemple 4
//
xset("window",3);xbasc();
// definition de la fonction
deff("yprim=f4(t,y)","yprim=(2+cos(t))*y-(y^2)/2-1");
Y0 = [0.39,1,2,3]; tmax = 5; pas = 0.001;
champty(f4,Y0,tmax,pas);









