function [ yp ] = DampedPendulumSystem( t,y,alpha )
% right hand side of an ode system y'=f(t,y)
%  In this case the system is rescaled model damped pendulum
% the entries of y are theta and d theta/dtau
% y is a column vector with two entires as is yp;
% alpha is the dimensionless rescaled friction parameter
yp=[y(2);-sin(y(1))-alpha*y(2)];
end

