epsi=1E-6;
[x,y]=BLsolve1(epsi,100);
hold off ;
plot(x,y,'.-b','LineWidth',1.5);
hold on ;
xe = [xi:h:xf] ;
plot(xe,yesatta(xe,epsi),'-r','LineWidth',0.5);
legend('Numerica');
axis([0,1,-2,1]);