y1(t):=rhs(dsolve({de1,y(0)=yo},y(t)));