%
% twoPointBvpDdirichletFD6_convergencePlot
%  
  

  clear, home, close all, format compact                     


  Nvec = [10 20 30 40 50 60 70 80 90 100];

  for i = 1:length(Nvec)
  
        e(i) = twoPointBvpDdirichletFD6(Nvec(i));
        h(i) = 1/(Nvec(i)-1);

  end
  
  subplot(1,2,1)
  loglog(Nvec,e,'b.-'), xlabel 'N'
  subplot(1,2,2)
  loglog(h,e,'b.-'), xlabel 'h'


  % [x,y] = ginput(2)
  % rho = ( log(y(2)) - log(y(1)) )/( log(x(2)) - log(x(1)) )