function yp = hmotadif (t, y)
f = [0;25];
k = [0 -1;2.5 0.25];
yp = f-k*y;
