Ok I found this post:
[differential equations - Julia's DifferentialEquations step size control - Stack Overflow]
and that helped me to figure out where I was going wrong… the solution was to exponentiate the state variables z[:,2:4] outside f_h!
Ok I found this post:
[differential equations - Julia's DifferentialEquations step size control - Stack Overflow]
and that helped me to figure out where I was going wrong… the solution was to exponentiate the state variables z[:,2:4] outside f_h!