def _eval_log_log_integrate(a, k, x):
    return np.log(a) / (k + 1) * x ** (k + 1)

def log_log_integrate(a, k, x1, x2):
    return _eval_log_log_integrate(a, k, x2) - _eval_log_log_integrate(a, k, x1
        )
partial_integrals = []
for a, k, x_lower, x_upper in zip(intercepts, slopes, x[:-1], x[1:]):
    partial_integrals.append(log_log_integrate(a, k, x_lower, x_upper))
total_integral = sum(partial_integrals)
Создайте интервал массива, пропорциональный функции или другому массиву