disp('question 2) b)')
getf('solequadiff.sci')
S=solequadiff(1000);
u=S(:,2); v=S(:,3);
printf(' \n u(0,5) est environ egal a %f\n...
avec une erreur absolue egale a %.10f\n',u(500),norm(v-u,'inf'))

disp('question 2) c)')
S=solequadiff(10);
x=S(:,1);
x=[0;x;1]; // on ajoute les abscisses extremes
u=S(:,2); v=S(:,3);
u=[1;u;0]; v=[1;v;0]; // on ajoute les ordonnees extremes
disp(' On trace le graphe de la solution sur la figure 1.')
xset("window",1), xbasc()
plot2d([x,x],[v,u],[5,-1],...
leg="solution exacte@solution discrete",rect=[0,0,1,1]);
xtitle('solution numerique et solution exacte pour n=10','x','u(x)')

disp('question 2) d)')
n=[5:5:30]';
erreur=zeros(n);
for l=1:length(n)
S=solequadiff(n(l));
u=S(:,2); v=S(:,3);
erreur(l)=norm(v-u,'inf');
end
disp('On etudie graphiquement l''erreur sur la figure 2.')
xset("window",2), xbasc()
plot2d("ll",n,erreur,-2); // courbe "log,log"
xtitle('Variation de l''erreur','n','erreur')