Я пытаюсь использовать функцию интерполяции Python, чтобы получить значение y для данного x, но я получаю сообщение об ошибке «raise ValueError («массивы x и y должны быть равны по длине вдоль оси интерполяции», хотя мои массивы имеют равные размер и форма (согласно тому, что я получаю, когда использую .shape в своем коде). Я совсем новичок в программировании, поэтому я не знаю, как проверить, что еще может отличаться в моих массивах. Вот мой код:
s = []
def slowroll(y, t):
phi, dphi, a = y
h = np.sqrt(1/3. * (1/2. * dphi**2 + 1/2.*phi**2))
da = h*a
ddphi = -3.*h*dphi - phi
return [dphi,ddphi,da]
phi_ini = 18.
dphi_ini = -0.1
init_y = [phi_ini,dphi_ini,1.]
h_ini =np.sqrt(1/3. * (1/2. * dphi_ini**2. + 1/2.*phi_ini**2.))
t=np.linspace(0.,20.,100.)
from scipy.integrate import odeint
sol = odeint(slowroll, init_y, t)
phi = sol[:,0]
dphi = sol[:,1]
a=sol[:,2]
n=np.log(a)
h = np.sqrt(1/3. * (1/2. * dphi**2 + 1/2.*phi**2))
s.extend(a*h)
x = np.asarray(s)
y = np.asarray(t)
F = interp1d(y, x, kind='cubic')
print F(7.34858263)
