Consider the boundary value problem with Neumann conditions d²u +u = sin(x) dx² du = 1 (2TT) = 2 dx where we care about finding a solution over the interval [0, 27]. If we attempt to solve this problem numerically using 4 equal-width subintervals, second-order symmetric ap- proximations to the second derivative d²u u(x + h) — 2u(x) + u(x − h) h² dx² and first order approximations for the derivatives du -(0)~ u(h) — u(0) h du dx (2πT) ~ u(2π) — u(2π – h) h dx 2π where h = what is the linear system that results?