Week 3 - Solving second order ODEs

Aim

Write a program in Matlab that will simulate the pendulum motion.Simulate the motion between 0-20 sec, for angular displacement=0,angular velocity=3 rad/sec at time t=0.

Theory

In the above equation,

g = gravity in m/s2,

L = length of the pendulum in m,

m = mass of the ball in kg,

b=damping coefficient.

L=1 metre,

m=1 kg,

b=0.05

Results

Not sucessfull

Solution of ODE

theta=theta1

dtheta/dt=dtheta1/dt=theta2

d^2theta/dt^2=d^2theta1/dt^2=dtheta2/dt

by this,

dtheta1/dt=theta2

dtheta2/dt=-(b/m)theta2-(g/l)sin(theta)

solution for theta

d/dt[theta;theta2]=[theta2;-(b/m)theta2-(g/l)sin(theta)]

g=9.81 m/s2.