In Matlab?

The following is equation of a motion of a simple pendulum d^2 theta/dt^2 + mg/L sin (theta) = 0 with m=1, g=9.8 and L=0.5, Linearize the model about the equilibrium point, theta = 0. Simulate both the nonlinear system and the linearized version for theta(0) = 0.1 and theta(0) = 0. Compare the two responses. Repeat for theta (0) = pi/4 and theta(0) = 0