please use mathematica for code NOT MATLAB (3) (20 points) The (dimensionless) equations of motion for a frictionless double pendulum system as shown below (in the figure on the left) with mi m2 and L1 L are The solutions are graphed below (on the right) for the initial conditions θι (0) 2, θ1(0) 1.02(0) 0, and 02(0)-0 for oS t s 50. (a) Reformulate the IVP as a first order system.2 (b) Generate approximate solutions using any method (Euler, improved Euler, Runge-Kutta, NDSolvo, ode45) and graph θί and θ2 as functions of t. Your graph should look like the one below (right). 81 L. 4π 02 01 40 50 (3) (20 points) The (dimensionless) equations of motion for a frictionless double pendulum system as shown below (in the figure on the left) with mi m2 and L1 L are The solutions are graphed below (on the right) for the initial conditions θι (0) 2, θ1(0) 1.02(0) 0, and 02(0)-0 for oS t s 50. (a) Reformulate the IVP as a first order system.2 (b) Generate approximate solutions using any method (Euler, improved Euler, Runge-Kutta, NDSolvo, ode45) and graph θί and θ2 as functions of t. Your graph should look like the one below (right). 81 L. 4π 02 01 40 50