clear all
%
% usa findif per fare il calcolo dell'errore (nella norma del vettore)
%
caso=2;
n=[2 4 8 16 32 64 128 256 512];
for i=1:length(n)
  [x sol]=findif(caso,0,1,n(i));
  normer(i)=norm(esatta(caso,x)-sol)/norm(esatta(caso,x));
end
for i=1:length(n)-1
  ordine(i)=-(log(normer(i+1)/normer(i))/log(n(i+1)/n(i)));
end
ordine
figure(1)
clf;
set(1,'defaultlinelinewidth',3)
set(1,'position',[407 12 612 684]);
loglog(n,normer)
