%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% dampydot.m
% Damped ydot.m
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function ydot=dampydot(t,y,d)
g=9.8;length=1;
ydot(1) = y(2);
ydot(2) = -(g/length)*sin(y(1)) - (d * y(2));