%% auteur : JJ Bataille qui n'aurait pu rien faire sans JP Vignault et
%% tous les autres....

%4 setxyrapport
100 setxunit
-3.5 3.5 setxrange
-1.5 2 setyrange
1 1 settkstep
.5 .5 setsubtkstep
2 2 setmkstep
ticks
subticks
marks
traceaxes
unites

%%%% L'équation différentielle y'=y+x avec la condition initiale %%%%y(0)=0 sur lintervalle [-3;3]%%%%

%%% La foction solution%%%%%%%%%
/f {setxvar
#rpn#   Sin (x)
} def

%%%%%%%%% Paramètres pour l'approx par la méthode d'Euler%%%%%%%
/nb 4 def	     %%%nbre d'intervalles de subdivision
/pas {2 nb div} def  %%%calcul du pas de la subdivision
/x0 {0} def          %% coordonnees du point a l'origine
/y0 {0} def
/a  {-3} def          %% intervalle de calcul = [a;b]
/b  {3} def

%%%%% L'exécutable qui calcule (x_k+1, y_k+1) en fonction de (x_k,y_k)%%%%%%

/F {
2 dict begin
   /y exch def
   /x exch def
#rpn# Cos (x)
end
} def

%%%%%%%% première execution et tracé en rouge %%%%%%%%%%
rouge
[a b {F} x0 y0 pas Euler ] ligne  

%%%%%%%%% deuxième exécution et tracé en vert %%%%%
/nb 20 def
vert
[a b {F} x0 y0 pas Euler ] ligne  

%%%%%%%%%% Tracé de C_f en bleu %%%%%%%%%%  
bleu
a b {f} Courbe

%%%%%%%%% Labels %%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%% Label 1 %%%
<latex>
$f:x\mapsto \sin x$\\
Solution de $\left\{\begin{array}{l}
y'=\cos x\\
y(0)=0 
\end{array}\right.\mbox{ sur } [-3;3]. $
</latex>

%% placement du label 1
noir
/fillstyle {jaune fill} def
boxit
%% #rpn# (xmin+xmax)/2 %% pour mettre au milieu du dessin
#rpn# xmin/2
1 cctexlabel
%%%%%%%% Label 2 %%%%%%%
<latex>
\begin{minipage}{8 cm}{
\begin{itemize}
\item $\mathcal{C}_f$ est  tracée en bleu.
\item en rouge et en vert ses approximations
par la méthode d'Euler avec des pas respectifs de 0,5 et de 0,1
\end{itemize}}
\end{minipage}
</latex>

%% placement du label 2
/fillstyle {blanc fill} def
boxit
#rpn# (xmax)/2
-.5 cctexlabel