solve_poisson.m

   1 function u = solve_poisson(integerN, f)
   2   ##
   3   ## Numerically solve the poisson equation,
   4   ##
   5   ##   -u_xx(x) = f(x)
   6   ##
   7   ## in one dimension, subject to the boundary conditions,
   8   ##
   9   ##   u(0) = 0
  10   ##   u(1) = 1
  11   ##
  12   ## over the interval [0,1]. It is assumed that the function `f` is at
  13   ## least continuous on [0,1].
  14   ##
  15
  16   [xs,h] = partition(integerN, a, b);
  17 end