PYTHON

(the part about using the function from Problem 1 is not necessary)

The procedure for finding the Drazin inverse using the Schur decomposition is given in Algo- rithm 14.1. Due to possible floating point arithmetic errors, consider all eigenvalues smaller than a certain tolerance to be 0. 2: 3: Algorithm 14.1 1: procedure DRAZIN(A, tol) (n,n) + shape(A) T1, Q1, ki + schur(A, |2| > tol) T2, Q2, k2 + schur(A,|2| tol) T2, Q2, k2 + schur(A,|2|