y1(t):=dsolve({de2,y(0)=yo},y(t),numeric);