#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Tue Sep 24 16:26:42 2024

@author: sam
"""

import numpy as np
from scipy import linalg as la



def qreigs(A,tol):
    for i in range(100000):
        Q,R = la.qr(A)
        
        A = R@Q
    lam = 0
    for i in range(len(A)):
        lam += A[i][i]
    return lam



print('lam =',qreigs(np.random.rand(5,5),1e-5))