function dy=f(t,y)

k=1;
dy = (20-y)*k;
