Solution of ordinary differential equations MATLAB code