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