Attempted to add matrix expression for rate but not succesful. ODE solver succesfully used