import random
import numpy as np

def GramSchmidt(*a):
    k=len(a[0])
    N=len(a)
    b = [[0] * k for i in range(N)]
    b[0]=a[0]
    for i in range(1,N):
        sum=a[i]
        for j in range(0,i):
            scolar_ab=0
            scolar_bb=0
            proj=[i for i in range(k)]
            for n in range(k):
                scolar_ab+=b[j][n]*a[i][n]
                scolar_bb+=b[j][n]*b[j][n]
            for n in range(k):
                proj[n]=(scolar_ab/scolar_bb)*b[j][n]
            for n in range(k):
                sum[n]-=proj[n]
        b[i]=sum
    return b

l1=np.array([1,-1,1,1])
l2=np.array([3,-2,1,2])
l3=np.array([-1,-2,1,2])
l4=np.array([1,-2,-1,2])
print('gs=',GramSchmidt(l1,l2,l3,l4))