midpoint method function for ODE

Just wondering how I could modify my Euler's function to get a Midpoint Method function I tried changing the equation inside the for loop to mirror the Midpoint algorithm y(i+1) = y(i) + dydt(t(i+1/2), y(i+1/2))*(t(i+1)-t(i)) but adding the halves result in a logical error