Подскажите пожалуйста, как реализовать в python интеграл, представленный на рисунке (или 18-19 строки в коде)?
Пытался через sympy.integrate, но потерпел крах - ValueError: Invalid limits given: (20.665,)
import numpy as np import sympy from scipy import interpolate b_2 = 6.67 l_k = 24.0 U_1 = [0.0, 4.3, 8.1, 11.4, 14.1, 16.3, 18.0, 19.3, 20.2, 20.7] V_1 = [0.0, 0.6, 1.8, 3.5, 5.7, 8.3, 11.3, 14.8, 18.6, 22.8] "Первая производная уравнения центроиды" def Cn_1(x): x_points = U_1 y_points = V_1 tck = interpolate.splrep(x_points, y_points) return interpolate.splev(x, tck, der = 1) "Расчет длины центроиды в зависимости от координаты X" def l_cn(x): return sympy.integrate((1 + Cn_1(x)**2)**0.5, x) "Полная длина центроиды" L_max = l_cn(l_k-0.5*b_2) print("L_max = ", L_max)