function error = twoPointBvpDdirichletFD6(N)

    x = linspace(0,1,N)';
   dm = dmSecond7Pt(x);

     alpha = 0.2;       %  u (0) = alpha
      beta = 0.9;       %  u (1) = beta
      
     f = exp(x);
     
     uExact = exp(x) + x*(beta - exp(1) - alpha + 1) + (alpha - 1);
     
      dm(1,:) = 0;     dm(1,1) = 1;     f(1) = alpha;      %  Dirichlet BC at x = 0
      dm(N,:) = 0;     dm(N,N) = 1;      f(N) = beta;      %  Dirichlet BC at x = 1

      uApprox = dm\f;
     
      error = norm(uExact-uApprox,inf);
%      plot(x,uApprox,'b--o',x,uExact,'r')