%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% forceydot.m
% Forced ydot.m
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function ydot=forceydot(t,y,d, force)
g=9.8;length=1;
ydot(1) = y(2);
ydot(2) = -((g/length) + (force * cos(8*pi*t)) ) * sin( y(1) ) - (d * y(2));