Diagonalize the matrix (IN MATLAB)

A =[9 4 2 4;56 32 28 44;14 14 6 14;42 33 21 45]

To do this, compute the eigenvalues of A using the command eigs. Then, compute a bases for the eigenspaces corresponding to the eigenvalues using rref on the appropriate augmented systems.

Display the eigenvalues, the bases for their corresponding eigenvalues, the matrix P and the diagonal matrix D.