程序出来了 就是答案不对

def triangleFun(self):
        a=self.a
        b=self.b
        row=a.shape[0]
        col=a.shape[1]
        n=row-1
        for i in range(0,n):
            for j in range(i+1,n+1):
                if a[i,i]==0:
                    continue
                lam =float(a[j,i])/a[i,i]
                a[j,0:n+1]=a[j,0:n+1]-lam*a[i,0:n+1]
                b[j,0:n+1]=b[j,0:n+1]-lam*b[i,0:n+1]
        for p in range(n,-1,-1):
            for q in range(p-1,-1,-1):
                if a[p,p]==0:
                    continue
                lam=float(a[q,p])/a[p,p]
                a[q,0:n+1]=a[q,0:n+1]-lam*a[p,0:n+1]
                b[q,0:n+1]=b[q,0:n+1]-lam*b[p,0:n+1]
        for l in range(0,n+1):
            lam=a[l,l]
            if lam==0: