voila une exemple d'ulitisation de la differentiation automatique
pour la resolution d'un Laplacien par difference fini
Si vous utilisez la classe R2, ou bien vous enlevez les operator<< de la classe R2 (car il y a conflit avec les operator<< de ddouble.h)
ou bien vous enlevez les operator<< de la class ddouble.

$$- \mu \Delta u =1 dans \Omega $$
$$  u = 0 sur \Gamma $$

But: Le calcul de  du/d\mu   pour \mu=1

compilation:

make

execution: 
./GradConjugueDA


Visualisation:

gnuplot
gnuplot> splot "pp.gp" u 1:2:3 w l t "u"        
gnuplot> splot "pp.gp" u 1:2:4 w l t "du/dmu"  
quit


