
import numpy as np

# Rozwiązuję poniższy układ równań
print ('1a + 1b = 35')
print ('2a + 4b = 94')

A = np.matrix([[1,1],[2,4]])
B = np.matrix([[35],[94]])

# Znajduję macierz odwrotną do A
A_inverse = np.linalg.inv(A)
print (A_inverse)

# Szukam X
X = A_inverse * B
print (X)

print('\n' * 3)
print('Kolejny przykład:')
print('\n' * 3)

# Rozwiązuję poniższy układ równań
print ('1a + 1b = 35')
print ('2a + 4b = 94')

# Tworzę równania
a = np.array([[1, 1],[2,4]])
b = np.array([35, 94])

# Wyświetlam rozwiązanie
print (np.linalg.solve(a, b))

