Fichier ode_05_rungekutta.jps — Modifié le 11 Mars 2006 à 10 h 53

ode_05_rungekutta.pdf
Source
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