A = bsr_matrix((T_Arm[:n].flatten(), (row, col)), shape=(5 * n, 2 * n))
A.toarray()
array([[0, 1, 0, 0, 0, 0], [2, 3, 0, 0, 0, 0], [4, 5, 0, 0, 0, 0], [6, 7, 0,
    0, 0, 0], [8, 9, 0, 0, 0, 0], [0, 0, 10, 11, 0, 0], [0, 0, 12, 13, 0, 0
    ], [0, 0, 14, 15, 0, 0], [0, 0, 16, 17, 0, 0], [0, 0, 18, 19, 0, 0], [0,
    0, 0, 0, 20, 21], [0, 0, 0, 0, 22, 23], [0, 0, 0, 0, 24, 25], [0, 0, 0,
    0, 26, 27], [0, 0, 0, 0, 28, 29]], dtype=int64)
b = Erg[:n].flatten()
Solve a multitude of linear least square system efficiently