PlotEssai04.m

function PlotEssai04(FileName)
T=load(FileName);
loglog(T(:,1),T(:,2),'rd--',T(:,1),(T(1,2)/T(1,1))*T(:,1),'k-.',T(:,1),T(1,2)/T(1,1)^2*T(:,1).^2,'k-')
xlabel('h')
legend('quadPM','O(h)','O(h^2)')
I=strfind(FileName,'.');
if isempty(I)
  SaveFig=[FileName];
else
  SaveFig=[FileName(1:I(end)-1)];
end
fprintf('Save figure in file %s.png ...\n',SaveFig);
print('-dpng',[FileName])