1 function u = solve_poisson(integerN, f)
3 ## Numerically solve the poisson equation,
7 ## in one dimension, subject to the boundary conditions,
12 ## over the interval [0,1]. It is assumed that the function `f` is at
13 ## least continuous on [0,1].
16 [xs,h] = partition(integerN, a, b);