function dvt = pendu(t,y)
global c

dvt = [y(2,:); -c*y(2,:) - sin(y(1,:))];

