Using scipy solve_ivp for equation with inital conditions, updating variables problem

Question

The code I have below attempts to map out the trajectory of a spaceprobe around the moon. I however, do not know how to update the variables for position and velocity which