I am new to matlab and am trying to solve a boundary value problem of the form Y"=1/L^2(Y-f(x)) where my f(x) comes from a griddedInterpolant using the xmesh for the pr