TypeError: unsupported operand type(s) for

import numpy as np
def _print_mat(x):
    for row in x:
        print(row)
n = int(4); 
kp = p = 1
ak = np.array ([[4, 3, 7, 2,],[1,3,5,6,],[5,7,1,8,],[6,7,8,1,]])
 
a = [[2,4,1,1, 1,0,0,0,],
     [0,2,1,0, 0,1,0,0,],
     [2,1,1,3, 0,0,1,0,],
     [4,0,2,3, 0,0,0,1,]]
 
e = xo = np.array([[0,0,0,0,],[0,0,0,0,],[0,0,0,0,],[0,0,0,0,]])
print('Вычисление обратной матрицы - XO')
print('Классическим методом Гаусса - Жордана')
print(a)
 
for i in range(0, n):
    aii = a[i],[i]
    for j in range(i+1 , n+i+1):
        a[i][j]=a[i][j]/aii
    a[i],[i]=1
    for k in range(0 , n):
        if k !=i:
            kp = a[k],[i]
            for j in range(i+1 , n+i+1):
                a[k , j]=a[k , j]-a[i , j]*kp
            a[k , i]=0
print('Преобразованная матрица A DLY XO')
print(a)
print('Обратная матрица XO')
for i in range(0,n):
    for j in range(0,(n)):
        xo[i , j]=a[i , n+j]
print(xo)
e = np.dot(xo,ak)
print('Единичная матрица E = XO * AK')
print(e)
print('конец программы')

import numpy as np
def _print_mat(x):
    for row in x:
        print(row)
n = int(4); 
kp = p = 1
ak = np.array ([[4, 3, 7, 2,],[1,3,5,6,],[5,7,1,8,],[6,7,8,1,]])
 
a = [[2,4,1,1, 1,0,0,0,],
     [0,2,1,0, 0,1,0,0,],
     [2,1,1,3, 0,0,1,0,],
     [4,0,2,3, 0,0,0,1,]]
 
e = xo = np.array([[0,0,0,0,],[0,0,0,0,],[0,0,0,0,],[0,0,0,0,]])
print('Вычисление обратной матрицы - XO')
print('Классическим методом Гаусса - Жордана')
print(a)
 
for i in range(0, n):
    aii = a[i][i]
    print(aii)
    for j in range(i+1 , n+i+1):
        a[i][j]=a[i][j]/aii
    a[i][i]=1
    for k in range(0 , n):
        if k !=i:
            kp = a[k][i]
            for j in range(i+1 , n+i+1):
                a[k][j]=a[k][j]-a[i ][j]*kp
            a[k][i]=0
print('Преобразованная матрица A DLY XO')
print(a)
print('Обратная матрица XO')
for i in range(0,n):
    for j in range(0,(n)):
        xo[i][j]=a[i][n+j]
print(xo)
e = np.dot(xo,ak)
print('Единичная матрица E = XO * AK')
print(e)
print('конец программы')

import numpy as np
def _print_mat(x, form):
    for row in x :
        print(*map(form.format, row))
n = int(4); 
kp = p = 1
ak = np.array ([[4, 3, 7, 2,],[1,3,5,6,],[5,7,1,8,],[6,7,8,1,]])
 
a = [[2,4,1,1, 1,0,0,0,],
     [0,2,1,0, 0,1,0,0,],
     [2,1,1,3, 0,0,1,0,],
     [4,0,2,3, 0,0,0,1,]]
 
e = xo = np.array([[0,0,0,0,],[0,0,0,0,],[0,0,0,0,],[0,0,0,0,]])
print('Вычисление обратной матрицы - XO')
print('Классическим методом Гаусса - Жордана')
form = '{:2d}'
_print_mat(a, form)
 
for i in range(0, n):
    aii = a[i][i]
    print(aii)
    for j in range(i+1 , n+i+1):
        a[i][j]=a[i][j]/aii
    a[i][i]=1
    for k in range(0 , n):
        if k !=i:
            kp = a[k][i]
            for j in range(i+1 , n+i+1):
                a[k][j]=a[k][j]-a[i ][j]*kp
            a[k][i]=0
print('Преобразованная матрица A DLY XO')
form = '{:.4f}'
_print_mat(a, form)
print('Обратная матрица XO')
for i in range(0,n):
    for j in range(0,(n)):
        xo[i][j]=a[i][n+j]
form = '{:2d}'
_print_mat(xo, form)
e = np.dot(xo,ak)
print('Единичная матрица E = XO * AK')
_print_mat(e, form)
print('конец программы')