Please show code

The equation of motion for a simple pendulum is +wsin(O)=( (*). For small angles the pendulum performs harmonic motion of frequency w . For this exercise, assume that angles not small, and w =2.0 . (a) Write a MATLAB script, which solves (*) for arbitrary initial conditions 0(0)=0, and 0(0)=0. plots both the position and the speed, and saves the time series. (b) Using (a), experiment with initial conditions in the range OELI, 5] and zero initial speed and select a few to show differences. Feel free to experiment with the time range to capture the character of the solution. What happens near 0=r ? Repeat with zero initial position, but with 0. [3.0,4.0] . (c) For the selections made in (b), plot the phase planet orbits on the same plot. is (d) Repeat (a) - (c) introducing viscous damping so that the equation of motion +0.10+2 sin (O)=0 (d) Discuss the results. Optional: Plot the solution of the most interesting" case in the 3D using position, velocity and time as coordinates