syms y(t) y(t) = dsolve(diff(y,t) == t - 2*t*y)