clear allsyms x y z[y,z]=dsolve('Dy=3y-2z','Dz=2*y-z','x')
y =
C1exp(x) + (C2exp(x))/2 + C2xexp(x)
z =
C1exp(x) + C2x*exp(x)