Write an m-file that implements the Modified Euler's method. It should work for a system with any number of states. It should accept a function handle for the "right hand side" dynamics, a time vector, and an initial state. It should return the numerically integrated state values at all points in the time vector.