function dydt = diffegy(t,y)
dydt = -y + t + 1;
