syms y(t) y(t) = dsolve(diff(y,t) == 3*t*(y^2) - t*y); y(t) = simplify(y)