clear all
close all
a=-pi/3;b=pi/4;nq=100;
u=@(x) cos(x);
d2u=@(x) -cos(x); % derivee seconde de u
nu=2;
f=@(x) -d2u(x)+nu*u(x);
ga=u(a);gb=u(b);

[x,U]=ep00.solveBVP01v2(a,b,f,nu,ga,gb,nq);

figure(1)
plot(x,U,'b',x,u(x),'r')
legend('sol. num.','sol. ex.')

figure(2) % représentation de l'erreur
plot(x,abs(U-u(x)))
title('error')
