Fichier ode_05_rungekutta.jps — Modifié le 11 Mars 2006 à 10 h 53
autocrop
-5 5 setxrange
-7 5 setyrange
-5 5 setyrange
30 setxunit
gsave
masque
quadrillage
/f {
2 dict begin
/y exch def
/x exch def
x 0 eq {
0
} {
y 3 exp y add x div
} ifelse
end
} def
.5 setlinewidth
/arrowscale {.5 dup} def
.4 setgray
%% syntaxe: f s1 s2 l champvecteur --> dessin du champ de vecteur
{f} .25 .25 .1 (-) champvecteur
/arrowscale {.75 dup} def
.7 setlinewidth
noir
tracerepere
.3 setlinewidth
%% on represente quelques solutions
bleu
-6 1 -1 {
/i exch def
[.1 xmax {f} i 1 .01 Rungekutta] ligne
[xmin -.01 {f} i 1 .01 Rungekutta] ligne
} for
rouge
1 1 6 {
/i exch def
[.1 xmax {f} i 1 .01 Rungekutta] ligne
[xmin -.01 {f} i 1 .01 Rungekutta] ligne
} for
grestore
#tex# $xy' - y = y^3$
0 -6 [2 dup] cctexlabel