Jupyter python

Test

speed = 0.1 b = 0.15 xoft, yoft = trajectory(b, speed) assert (xoft[1]>0.145 and xoft[1]-1.97 and yoft[1] V(x, y) = xv e-(x+y?) (all terms dimensionless) with four maxima centred at x = +1 and y = +1, respectively. It is defined as a repulsive potential for any test particle (the pinball) moving through that potential. The equations of motion to solve hence are = -2.0 y x(1 - x) e-(x2+y2) dt2 and dx(t) dy(t) =-2.0 x y(1 12) e-(x2+x?). dt2 Task: Solve these two coupled ODE's in stages: (a) write a function trajectory(impactpar, speed) which takes the impact parameter, x(O), as input as well as the speed v, for initial conditions (impact dx(0) dy(0) parameter= x(0), Ux = y(0) and vy = in that order). Calculate the time array to solve for from np.linspace(0,maxtime, 300) where di maxtime works well as 10 / speed. The function shall return the trajectory, i.e. the arrays of solved x(t) and y(t) separately. Plot a trajectory of your choice, see below, with appropriate axes labels. di As initial values, fix vx = dx(0) = 0.0 and y0) = -2 (away from the potential in y). The range of sensible initial values for x(O), should be between dt -0.9