Ошибка при работе с матрицами

import numpy as np
 
C = np.matrix([[6.22, 1.42, -1.72, 1.91],[1.42, 5.33, 1.11, -1.82],
             [-1.72, 1.11, 5.24, 1.42],[1.91, -1.82, 1.42, 6.55]])
B = np.matrix([[7.53], [6.06], [8.05], [8.06]])
b = np.matrix([[0],[0],[0],[0]])
for i in range(0,4):
    b[i] = B[i]/C[i][i]
print(b)

b = np.matrix([[0.],[0],[0],[0]]) # матрицы видимо типизированные
for i in range(0,4):
    b[i] = B[i]/C[i,i]; # использовать запятую вместо двух [i], 
#  ибо похоже результатом [i] является матрица с одной строкой и C[0][0][0][0][0]... = С[0]

a = np.array([[6.22, 1.42, -1.72, 1.91], [1.42, 5.33, 1.11, -1.82],
               [-1.72, 1.11, 5.24, 1.42], [1.91, -1.82, 1.42, 6.55]]) # ваша матрица С
 
C = np.array([[a, 2*a],[3*a, 4*a]]) # Новая 4 мерная таблица
 
print('Table 4 dimensions: \n', C) # Показывает вам всю таблицу
 
print(C[1][0]) #  Показывает вам первую подтаблицу из 2ой dimension
print(C[1][1][3:]) # Показывает вам 3ию и ниже строки второй подтаблицы из 2ой dimension
print(C[1][0][2][3]) # Показывает конкретный элемент.

print(C[:,:,2,:])
print(C[1,0,:,:])
print(C[1,1,3:,:])
print(C[1,0,2,3])