a=0; b=1;
n=input('Quanti sottointervalli? ');
e=input('Valore di epsilon? ');
i=input('Schema sdandard (1) o mass-lumping (0)? ');
[x,sol]=diffreaz(i,e,a,b,n,'f');
figure(1)
clf;
set(1,'defaultlinelinewidth',3)
set(1,'position',[407 12 612 684]);
fplot(['esatta(' num2str(e) ',t)'],[a b],'r')
hold on;
plot(x,sol,'*--')
h_legend=legend('esatta','calcolata',4);
h_text=findobj(h_legend,'type','text');
set(h_text,'FontUnits','points','FontSize',20);
