M�thode d'EULER
1) On y va lentement et s�par�ment
d�finition des fonctions et r�solution exacte des �quations diff�rentielles
> |
restart:with(DEtools):with(plots):f:=(x,y)->(-(x^2)*y-y^2+2*x)/(1-x^3);eq:= D(y)(x)=f(x,y(x)); |
Warning, the name changecoords has been redefined
proc�dure de calcul des
par la m�thode d'Euler
> |
z:=proc(n,xo,yo,f,A) |
initialisations et calcul des
par appel de la proc�dure
d�termination des extr�mums des
> |
with(plots): |
trac�s simultan�s des 4 courbes
> |
with(plottools):seq_opt:=x=-0.05..xo,y=min(a1,b1,c1,a2,b2,c2)..max(a1,b1,c1,a2,b2,c2)+0.05,color=black: |
Warning, the names arrow and translate have been redefined
2) Le tout en un
> |
restart:with(plots):with(DEtools):with(plottools): |
Warning, the name changecoords has been redefined
Warning, the names arrow and translate have been redefined
> |
z:=proc(n1,n2,n3,xo,yo,f) |
On essaie
> |
> |
f:=(x,y)->(-(x^2)*y-y^2+2*x)/(1-x^3);eq:= D(y)(x)=f(x,y(x)); |