Это часть гораздо большего кода, над которым я работаю. У меня есть набор из примерно 100 точек данных (p, l), из которых я получаю требуемое уравнение с помощью сплайновой интерполяции cubi c, где l - диапазон, а p - импульс. Мотивация для кода заключается в том, что для данного импульса p мне нужно уравнение для производных p ', p' ', p' '', взятых по l.
p=data[:,2]
l=data[:,10]
range_l = CubicSpline(p,l,bc_type='natural') #range as a function of p
fl = CubicSpline(l,p,bc_type='natural') #cubic spline to get f(l) from p and l
fl_1=CubicSpline.derivative(fl,nu=1) #using it to obtain the differential f'(l)
fl_2=CubicSpline.derivative(fl,nu=2) #using it to obtain the differential f''(l)
fl_3=CubicSpline.derivative(fl,nu=3) #using it to obtain the differential f'''(l)
print(fl(range_l(1.25)),fl_1(range_l(1.25)),fl_2(range_l(1.25)),fl_3(range_l(1.25)) )
Это код, который я написал, делает то же самое, но получаемые мной значения не имеют смысла. Я новичок в python, поэтому будьте добры.