dsolve({diff(x(t),t)+(a*x(t)+b*y(t))*f(t) = g(t), diff(y(t),t)+(c*x(t)+d*y(t))*f(t) = h(t)})