disp(' '), disp(' question 2) b)')
S=solequadiff(1000);
u=S(:,2); v=S(:,3);
fprintf(' u(0,5) est environ égal à %f\n',u(500))
fprintf('avec une erreur absolue égale à %.10f\n',norm(v-u,'inf'))

disp(' '), disp(' '), disp(' question 2) c)')
S=solequadiff(10);
x=S(:,1);
x=[0;x;1]; % on ajoute les abscisses extrêmes
u=S(:,2); v=S(:,3);
u=[1;u;0]; v=[1;v;0]; % on ajoute les ordonnées extrêmes
disp(' On trace le graphe de la solution sur la figure 1.')
figure(1), clf
plot(x,v,'r',x,u,'k+')
title('solution numérique (par des croix) et solution exacte (en rouge), pour n=10')
xlabel('x')
ylabel('u(x)')

disp(' '), disp(' '), disp(' question 2) d)')
n=[5:5:30]';
erreur=zeros(size(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 étudie graphiquement l''erreur sur la figure 2.')
figure(2), clf
loglog(n,erreur,'kx');
title('Variation de l''erreur')
xlabel('n')
ylabel('erreur')