Вы не добавили два последних уравнения вашей задачи в матрицу A
:
import numpy as np
A = np.array([[1,3,2,2],[2,1,1,2],[3,1,2,1],[2,0,1,3]])
B = np.array([1,0,1,0])
X2 = np.linalg.solve(A,B)
Дает:
array([-0.27272727, -0.18181818, 1.09090909, -0.18181818])
Это должно работать.