Need help with part c!

The American League of Major League Baseball has decided to replace the Designated Hitter (the player who bats for the pitcher) with a Robotic Designated Hitter (RDH) for the 2021 season. Careful measurements of the adopted model of robotic hitter reveal that the distribution of ball speed off the bat is a Gaussian of mean 100 miles per hour and standard deviation 15 mph, and the distribution of launch angle is a Gaussian of mean 45 degrees and standard deviation 10 degrees. The American League is concerned that the new robotic hitters will so outperform human players that they will face a backlash from fans. Your job is to quantify how good a home-run hitter this RDH actually is by computing the expected ratio of at-bats to home runs, the AB/HR-ratio. Top human hitters have AB/HR-ratios of 10 or so (Babe Ruths was 11.76; note that better hitters have smaller AB/HR ratios).

To compute the expected AB/HR-ratio for a RDH, model the ball as it leaves the bat as a projectile moving under the combined forces of gravity and air resistance. If the motion in confined to the x y plane, with y being the height above the ground, the motion of the ball is governed by the equations

dv/dt = -g - v (C_d A v)/2m and

dr/dt = v

where the state of the ball is given by r = (x,y) and v = (v_x,v_y). is the unit vector in the y direction.

[Note Feb 21, 2023: corrected a typo above where dr/dt = dv/dt instead of dr/dt = v was given above.]

Assume the baseball has a mass of 0.145 kg and a diameter of 7.4 cm. The additional parameters are g, the acceleration of gravity at the Earths surface (9.81 m s²), the density of air (1.2 kg m³), and C_d is the balls drag coefficient, which you can take to be equal to the constant 0.35. Finally, A is the cross-sectional area of the ball and v = | v |. Assume the ball is hit from a starting height of 1 m.

a) Solve the projectile problem above by writing a Python function that implements the Euler, Euler-Cromer and Midpoint methods. Required input should be the initial ball speed (in m/s), the launch angle (in degrees), the time-step used to numerically integrate the equations of motion (in seconds, called in your textbook), and which of the three methods (Euler, Euler-Cromer or Midpoint) should be used for the integration. It would be useful to be able to turn off air resistance; handle this with an argument to your function. Your routine should return the horizontal range of the ball as an array(in metres).

As a test, run your routine for the initial conditions of velocity = 15 m/s, theta=45 degrees, tau=0.1 using the Euler method. Keep in mind that the trajectory shown for Eulers method in this figure is computed without air resistance.

Plot a figure showing the outputs from the Euler, Euler-Cromer and Midpoint methods with a time-step of 0.1.

b) Determine the AB/HR-ratio for the proposed RDH. To do this, run your projectile program for a number of simulated at-bats with the starting v and drawn from the random distributions given above. Remember that a normally-distributed random variable of mean and standard deviation , N(, ²), can be obtained by scaling the output of the numpy.random.randn function, N(0, 1), via sigma * np.random.randn(...) + mu

You can assume any projectile with a range of 400 feet or greater is a home run as this is a typical distance between home plate and the centre field.

c) Real stadiums have fences and their height varies a lot: for example, the Green Monster left field fence in Bostons Fenway Park is 11.3 m high, while the outfield fence in London's Labatt Park is only about 2 m high. Show how the AB/HR-ratio for the RDH varies as fences from 0.5 to 15.0 m are placed at the 400 foot mark. If you want to keep the AB/HR-ratio greater than 10, how high do you need to place the fence, or is it always higher than 10 ? (Hint: start by modifying your routine of Part 1 to instead return the height of the ball at a distance of 400 feet note that this could be negative).