function yprime = f501(t, y)
yprime = -0.1*(y-10);
