# -*- coding: utf-8 -*-
"""
Created on Sun Jun 15 14:31:27 2025

@author: franc
"""

import numpy as np
import numpy.linalg as alg

A=1/6*np.array([[13,7,4],[7,13,4],[4,4,16]])
print(alg.det(A))
B=np.array([1,2,3])
print(alg.solve(A,B))
M=1/6*np.array([[13,0,0],[0,13,0],[0,0,16]])
N=M-A
Minv=6*np.array([[1/13,0,0],[0,1/13,0],[0,0,1/16]])


def itere(X0,n):
    l=[X0]
    X=X0
    for k in range(n+1):
        X=np.dot(Minv,np.dot(N,X)+B)
        l.append(X)
    return l
    
print(alg.eigvals(np.dot(Minv,N)))

X0=np.array([1,1,1])