Linearly Independent Solutions
> lambda := 'lambda':
> charpoly( A, lambda );
> eigenvalues( A );
>
> eigenvectors( A );
>
Case 1: Real and Distinct Eigenvalues
Case 2: Complex Eigenvalues
Case 3: Repeated Eigenvalue
>
The following two-parameter family of solutions is the general solution to the system of ODEs
> SOLgen := c1 * Y1 + c2 * Y2;