What might be causing the error when solving differential equations in MATLAB