# -*- coding: utf-8 -*-
"""
Created on Thu Dec 30 15:54:22 2021

@author: Sim
"""
import numpy as np
from numpy.linalg import solve, inv

A = np.array([[1,6,9], [8,2,4],[7,5,3]])
B = np.array([[1,2,3], [4,5,6], [7,8,9]])
b = np.array([40,24,26])

print(A@B, B@A)

x1 = inv(A)@ b
print(x1)

x2 = solve(A,b)
print(x2)

Ainv=inv(A); Binv = inv(B)
print(np.round(Ainv,3), np.round(Binv,3))