1. Approximating ODEs [40 points] Imagine you dug a tunnel straight through the center of the earth to the other side, and jumped in. You will experience acceleration due to gravity (g), but this will decrease as you approach the center of the planet (ignore the molten core) and then start increasing again as you start to move to the surface on the other side. You will also experience a drag force proportional to your velocity in the direction opposite your motion. The governing equation of motion is expressed by the following ODE for velocity: dv m- dt msign(r)g(r) sign(v) cav, where vis velocity (m/s); t = time (s); g(r) is the changing gravitational constant (m/s) as a function of height above the Earth's center (r, in m); sign(x) is a function that returns -1, 0, and 1 for negative, zero, and positive values of x, respectively (you can use the Python function numpy. sign() for this); ca is the drag coefficient (kg/m); and m is mass (kg). Acceleration due to gravity is given by x 4 3 g(r) = Gpr, where G = 6.674 10-11 mkgs-2 is the gravitational constant and p = average density of the Earth. 5515 kg/m is the In both cases, calculate the position and velocity of the jumper for t = 0 to 3 105 s (a x bit over 80 hours) given the following parameters: m = 65.8 kg, and ca=0.25 kg/m. Assume that the initial height is the radius of the Earth, 6.371 x 106 m, and zero initial velocity. Hint: Start by writing the equation in terms of position. Could we maybe turn that into..a system? a) Try using the midpoint method (i.e., modified Euler method) with a step size of 0.5 s. (Perhaps just solve to an end time of 100s.) Can you obtain a reasonable solution? What does this tell you?