r = (b**2+z**2)**(1/2)
rho_0 = 0.17
r_0 = 1.12
d = 0.54
def WS(A,r):
rho = rho_0/(1+np.exp((r-R(A))/d))
return rho
def R(A):
R = r_0*A**(1/3)-0.86*A**(-1/3)
return R
def collision(A_particle,B_particle,d):
t = lambda r: WS(A_particle,r)*WS(B_particle,r)[::-1]
T = scipy.integrate.quad(t,0,d)[0]
---> 21 t = lambda r: WS(A_particle,r)*WS(B_particle,r)[::-1]
22 T = scipy.integrate.quad(t,0,d)[0]
23 print('az ütközések száma: '+
IndexError: invalid index to scalar variable.
Если я не использую перевернутый массив, он работает нормально.
Я проверил их тип и длину, они совпадают, и я не знаю, как это исправить. Также пытался написать новую функцию, которая возвращает перевернутый массив, но это тоже не сработало.
def WSrev(A,r):
rho = rho_0/(1+np.exp((r-R(A))/d))
return rho[::-1]