Сеть Хопфилда

import numpy as np
n = 30
m = 3
           #|4|7|0|
E=np.array([[-1,-1,1],  
           [-1,1,1],
           [-1,1,1],
           [1,1,1],
           [1,1,1],
           [-1,-1,1],
           [-1,-1,-1],
           [1,-1,-1],
           [1,-1,-1],
           [1,1,1],
           [-1,-1,1],
           [1,-1,-1],
           [1,-1,-1],
           [-1,1,-1],
           [1,-1,1],
           [-1,-1,1],
           [1,-1,-1],
           [1,1,-1],
           [1,-1,-1],
           [1,-1,1],
           [-1,-1,1],
           [-1,-1,-1],
           [-1,1,-1],
           [-1,-1,-1,],
           [1,-1,1],
           [-1,-1,1],
           [-1,-1,1],
           [-1,1,1],
           [-1,-1,1],
           [1,-1,1]])
 
def h(x):
    if(x<0):
        res=-1
    else:
        res=1
    return res
 
def feedforward(x,w,n):
    y=np.zeros(n,float)
    for j in range(n):
        s=0
        for i in range(n):
            s+=w[i,j]*x[i]
        y[j]=h(s)
    return y
 
def ro(x,y):
    return np.sum(abs(x-y))/2
    
 
 
def teach(x,w,n):
    y1=x
    y2=feedforward(x,w,n)
    while (ro(y1,y2,n)>0):
        y1=y2
        y2=feedforward(y1,w,n)
    return y2
 
 
def ro(x,y,n):
    s=0
    for i in range(n):
        s+=np.abs(x[i]-y[i])
    return s
 
w=np.zeros((n,n), float)
for i in range(n):
    for j in range(n):
        w[i,j]=0.0
        if (i!=j):
            for k in range(m):
                w[i,j] += E[i,k]*E[j,k]
 
                
e1=E[0:n,0]
e2=E[0:n,1]
e3=E[0:n,2]
#x= [-1,-1,-1,1,1,-1,-1,1,1,1,-1,1,1,-1,1,-1,1,1,1,1,-1,-1,-1,-1,1,-1,-1,-1,-1,1]
#x= [-1,1,1,1,1,-1,-1,-1,-1,1,-1,-1,-1,1,-1,-1,,-1,1,-1,-1,-1,-1,1,-1,-1,-1,-1,1,-1,-1]    
#x= [1,1,1,1,1,1,-1,-1,-1,1,1,-1,-1,-1,1,1,-1,-1,-1,1,1,-1,-1,-1,1,1,1,1,1]    
y=teach(x, w, n)
r1=ro(e1,y,n)
r2=ro(e2,y,n)
r3=ro(e3,y,n)
print(r1)
print(r2)
print(r3)