function xNew = psi_RK2(t, xOld, xCur, tau)

xNew = xOld + tau*f_Extrap(t+tau/2,xCur);

end

Published with MATLAB® R2014a