I've run into an issue with a system of coupled 2nd order nonlinear ODEs. I've converted the system to a system of 1st order ODEs like normal, but its not in the format that ode45 requires. ode45 requires
x_dot = f(t,x), where x is the state vector
My form is
x_dot = f(t,x,x_dot)
where the function also depends on the derivative. I searched all around to see if anyone else has had this issue, and only found one post on a forum. It just said to proceed normally with ode45 and use the old values (previous step) of x_dot to calculate the rhs.
This makes sense, but ode45 doesn't work that way (not that I know of). I only have access to the old state values and not their derivatives.