function dydt=dampmass1(y)
% dydt=-0.348 - 3.011 *y -33.531 * y.^(3/2) - 3.858 * y.^2 +81.547 * y .^(5/2);
dydt=-0.348 - 3.011 *y -3.531 * y.^(3/2) - 3.858 * y.^2 +1.547 * y .^(5/2);