我试图通过改变很多变量来解决这个问题，但没有任何效果。下面是我的代码：

import numpy as np
import matplotlib.pyplot as plt
import math as mt
from numba import njit

N=1000
J=1
h=0.5

plt.rcParams['figure.dpi']=100

plt.xlabel('T', fontsize=14)

ns=100000000

s=np.empty([N,0])
for i in range(N):
    s[i]=np.random.choice([-1,1]) 
    

E=0
M=0
for i in range(-1,N-1):
    E=E-J*(s[i]*s[i+1]+s[i]*s[i-1]) 
    M=M+s[i]
    

Energy=np.empty([0,0])

C=np.empty([0,0])

Magne=np.empty([0,0])

chi=np.empty([0,0])
@njit
def average(k , E, M, ns, x, Energy, C, Magne, chi, s):
    Ev=0.
    Ev2=0.
    Mv=0.
    Mv2=0.
    
    for z in range(1,ns+1):
        
        i=np.random.randint(-1, high=N-1)
        
        dE=2*J*(s[i]*s[i+1]+s[i]*s[i-1])
        dM=2*h*s[i]
        
        pace=1./(1+mt.exp(k*(dE+dM)))
        #pace=min(1,mt.exp(-k*dE))
        
        
        if np.random.random()<pace:
            s[i]=-s[i]
            E=E+dE
            M=M-dM/h

        if z>x:
            Ev=Ev+(E-h*M)
            Ev2=Ev2+(E-h*M)**2
            Mv=Mv +M
            Mv2=Mv2 +M**2
        
        
    Ev=Ev/(ns-x)
    Ev2=Ev2/(ns-x)
    varE=(Ev2-Ev**2)
    Energy=np.append(Energy,Ev)
    
    C=np.append(C,varE*(k**2))
    
    Mv=Mv/(ns-x)
    Mv2=Mv2/(ns-x)
    varM=(Mv2-Mv**2)
    Magne=np.append(Magne,Mv)
    
    chi=np.append(chi,varM*k)
    
    return E, M, Energy, C, s, Magne, chi

b=np.arange(0.3,1.,0.01)

x=0.75*ns

for k in b:
    
    E, M, Energy, C, s, Magne, chi =average(k,E,M,ns, x,Energy,C, Magne, chi, s)

         
plt.plot(b,Energy, '.', color='r',linestyle='--')
plt.plot(b,C, '.', color='b',linestyle='--')
plt.title('Energy and heat capacity')
plt.legend(['E','C'])
plt.show()

plt.rcParams['figure.dpi']=100

plt.xlabel(r'$\beta$', fontsize=14)

plt.plot(b,Magne, '.',linestyle='--', color='r')
plt.plot(b,chi, '.', color='b',linestyle='--')
plt.title('Magnetization and susceptibility')
plt.legend(['M',r'$\chi$'])
plt.show()

字符串
我一直收到错误消息：

Cannot unify float64 and array(float64, 1d, C) for 'Mv2.3', defined at c:\users\usuario\documents\python scripts\ex13hw7.py (65)

File "ex13hw7.py", line 65: def average(k , E, M, ns, x, Energy, C, Magne, chi, s): <source elided> if z>x: Ev=Ev+(E-h*M) ^

型
在c：\users\usuario\documents\python scripts\ex13hw7.py中键入任务期间（65）
有人知道问题出在哪里吗？？
我尝试更改变量名，并在函数中包含/删除参数。