dsolve({diff(x(t),t) = -3*x(t)+48*y(t)-28*z(t), diff(y(t),t) = -4*x(t)+40*y(t)-22*z(t), diff(z(t),t) = -6*x(t)+57*y(t)-31*z(t)})