\documentclass[12pt]{article}
\usepackage{fancyvrb,pst-all}
\makeatletter
\newcommand{\executGiac}[1]{
\immediate\write18{giac <#1 } }
\makeatother
\begin{VerbatimOut}{Xcass.in}
maple_mode(0);
Sortie:=fopen("xcass.tex");
read("xcass.out");
Resultat:=cat(latex(ans()[0]));
fprint(Sortie,Unquoted,Resultat);
fclose(Sortie);
\end{VerbatimOut}
\begin{VerbatimOut}{xcasf.in}
Sortie:=fopen("xcasf.tex");
courbe:=read("xcasf.out"):;
Resultat:=cat(latex(courbe));
fprint(Sortie,Unquoted,Resultat);
fclose(Sortie);
\end{VerbatimOut}
\newenvironment{Xcass}
{\VerbatimEnvironment\begin{VerbatimOut}{xcass.out}}
{\end{VerbatimOut}
\executGiac{Xcass.in}
\[\input{xcass}\]}
\newenvironment{xcasf}
{\VerbatimEnvironment\begin{VerbatimOut}{xcasf.out}}
{\end{VerbatimOut}
\executGiac{xcasf.in}
\input{xcasf}}
\begin{document}
La solution de l'équation différentielle $y"+y=0$ avec les conditions particulières $y(0)=y'(0)=1$
est fournie par XCAS grâce à la commande \texttt{desolve}. Il s'agit de la fonction
définie par~:
\begin{Xcass}
desolve([y'+y=0,y(0)=2],y)[0];
\end{Xcass}
Nous obtenons également son graphe grâce à la commande \texttt{plot}~:
\begin{xcasf}
plot(desolve([y''+y=0,y(0)=1,y'(0)=1],y)[0],x=-Pi..Pi,color=blue);
\end{xcasf}
Merci à Yves Delhaye et Bernard Parisse.
\end{document}