Update:
Mathcad was able to determine the symbolic eigenvector matrix after all. I must have made an error in defining the problem the first time. However, the result was too large to display, so it can't be used for your purpose of seeing an entire symbolic solution.
As an alternative to a symbolic solution, you could obtain a parametric solution to the ODE system which would allow you to show how the plots change with different values of the parameters.