Write Matlab separate function files that will numerically solve Ordinary differential equations using the Midpoint and Heun's method. The function's should be in the form:

function [x,y] = midpt(dydx,tspan,y0,h)

function [x,y] = heun(dydx,tspan,y0,h)