Forcing an ODE solver to preserve the norm Posted by Krastanov, at scicomp.stackexchange.com, 15 Jul 2015 I have an ODE of the form $$ \frac{dy}{dt} = -i H y \enspace .$$ where $y$ is a complex vector and $H$ is a time dependent…