syms y(t) y(t) = dsolve(diff(y,t) == -(y + 3*t^2)/(t + 2*y))