function yprime = mcmc_dampharm(t,y,b,c);
yprime = zeros(2,1);
yprime(1) = y(2);
yprime(2) = -b*y(2)-c*y(1);