dsolve({diff(x(t),t) = a*x(t)+g*y(t)+beta*z(t), diff(y(t),t) = g*x(t)+b*y(t)+alpha*z(t), diff(z(t),t) = beta*x(t)+alpha*y(t)+c*z(t)})