Linearly Independent Solutions
> lambda := 'lambda':
> charpoly( A, lambda );
> eigenvals( A );
>
> eigenvects( A );
Case 1: Real and Distinct Eigenvalues
Case 2: Complex Eigenvalues
Case 3: Repeated Eigenvalue