def m_quadratic_sum(A, B, max_it=50):
 gamma1 = solve_discrete_lyapunov(A, B, max_it)
 return gamma1

def m_quadratic_sum ( A , B , max_it = 50 ) : gamma1 = solve_discrete_lyapunov ( A , B , max_it ) return gamma1
computes the quadratic sum .