function z = frhs(t, y)
z = t - y;