This is a question in python

The equations of motion for the golf ball are dtdvx=mFdrag,xmS0vydtdvy=mFdrag,y+mS0vxg, where y is the height of the ball above the ground, x is the horizontal direction in the same plane as the ball's initial velocity, and Fdrag is the atmospheric drag force given by FdragCAv2, similar to that in the projectile motion. is the angular velocity of the ball with back spin and spin axis perpendicular to both x and y.S0 is the coefficient of the Magnus force due to spin. (The Magnus force gives a large upward (vertical force) when a golf ball is hit with a significant amount of backspin.) S0 can be measured experimentally and S0/m0.25s1 for a typical case. It is believed to be a good approximation to assume that does not change significantly over the course of the trajectory. Thus we assume here a constant in our simulation. Assume an initial velocity of 70m/s (a typical value for an average player) and an initial angle of 9, write a program [50 pts] to solve the equations of motion using the Euler method. The golf ball has a diameter of 42.7mm and a mass of 45.9g. Use your program to study the ball's trajectories under the following hypothetical cases: (1) [10pts] Typical case: the drag coefficient C=0.5 for speeds up to v=14m/s, and C=7.0/v at higher velocities S0/m0.25s1