T=30;y0=0.3333333333;y1=y0;
N=evstr(x_dialog('valeur de N?','100'));
dt=T/N;
for i=1:N+1;
y0=y0+dt*(3*y0-1);
y1=y1+dt*(-3*y1+1);
end
disp([y0,y1])