I don't understand this problem at all, so if you could help out I'd appreciate it. Please write in python and explain how you did it.

In class you learned that the ODE describing the velocity gained by a falling body is dtdv=cgm(1v) For a body of mass m, acceleration of gravity g, and drag coefficient c. The closed form solution of the equation is v(t)=cgm(1e(c/m))=vT(1e(c/m)) where vT is the terminal velocity. We also learned in class that the velocity of the ball can be iteratively calculated by numerically solving the ODE using the formula v(ti1)=v(ti)+(gmcv(ti))t=v(ti)+(gmcv(ti))(ti1ti) where v(ti1) is velocity of the body at time ti1, and v(ti) at time ti.t is size of the time-step, and the above equation works if the time-step (t) is small enough. Write a program to calculate the velocity of the body at different time, as it accelerates to it's terminal velocity, using equation (3). You can assume g=10,m=10 and c=50. (a) Plot the variation of velocity with time using the mathematically derived equation for v (equation 2). (b) In the above figure, plot the variation of velocity with time calculated using equation 3. (c) What is the maximum time-step size t that allows accurate calculation of the velocity ? (d) Will the maximum allowable time-step size t change if the mass of the body (m) increases to m=100 ? Describe what you observe. (e) Theoretically, the falling body will take time t= to reach terminal velocity (vT). Using the program, calculate the time required by the body to attain 99.99 percent of the terminal velocity (vT)