We introduce a numerical method to solve stochastic optimal control problems which are linear in the control. We facilitate the idea of solving two-point boundary value problems with spline functions in order to solve the resulting dynamic programming equation. We then show how to effectively reduce the dimension in the proposed algorithm, which improves computational time and memory constraints. An example, motivated as an invest problem with uncertain cost, is provided, and the effectiveness of our method demonstrated.