function dydt=dampmass(t,y)
% dydt=-0.348 - 3.011 *y -3.531 * y^(3/2) - 33.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);