function dX = damped_pendulum_equat(t,X,b)

dX=zeros(2,1);

dX(1)=X(2);
dX(2)=-b*X(2)-sin(X(1));

