%
% chiama il programma p1p1 con il dato irregolare e disegna il grafico
% della soluzione
%
a=-1; b=1;
n=input('Numero di suddivisioni n? ');
h=(b-a)/n;
fig=input('Numero figura? ');
figure(fig)
set(fig,'position',[407 12 612 684])
clf
%
[x,s,u]=p1p12(a,b,n);
%
set(fig,'defaultlinelinewidth',2)
subplot(2,1,1)
fplot('esatta2(1,x)',[-1 1],'r')
hold on;
plot(x,s(1:n+1),'*--')
title('La derivata della soluzione','fontsize',15,'color','r')
%
subplot(2,1,2)
fplot('esatta2(0,x)',[-1 1],'r')
hold on;
plot(x,u,'*--')
title('La soluzione','fontsize',15,'color','r')
