Having issues towards the end where the functions to evaluate the symbolic matrices with given values at each time step of the ODE solver (using Euler method) being unable to properly substitute and evaluate the matrices with the given values.