Jul-07-2019, 01:19 PM
F(z) := [z[0]*(a-b*z[1]), z[1]*(c*z[0]-d)];
methodeEuler(F, Z0, h, N) := {
local Z, n;
Z := [Z0];
pour n de 0 jusque N faire
Z[n+1] := Z[n] + h*F(Z[n]);
ffaire;
return Z;
}:;
methodeEuler(F, Z0, h, N) := {
local Z, n;
Z := [Z0];
pour n de 0 jusque N faire
Z[n+1] := Z[n] + h*F(Z[n]);
ffaire;
return Z;
}:;