Всем добрый день!
Я хочу написать в питоне следующий интеграл:
Причем, r0 мне нужен в виде изменяемой переменной (он может быть не только 10). Гамма же равно 3.5. И по r границы зависят от r0. Если r0 = 1, то границы интеграла по r 1 - 1000, если r0 = 10, то 10 - 1000, если r0 = 100, то 100 - 1000.
Правильно ли я это делаю? Для меня тут проблемное место, чтоб интегрирование велось по r/r0 именно.
Я хочу написать в питоне следующий интеграл:
Причем, r0 мне нужен в виде изменяемой переменной (он может быть не только 10). Гамма же равно 3.5. И по r границы зависят от r0. Если r0 = 1, то границы интеграла по r 1 - 1000, если r0 = 10, то 10 - 1000, если r0 = 100, то 100 - 1000.
Правильно ли я это делаю? Для меня тут проблемное место, чтоб интегрирование велось по r/r0 именно.
Python:
import sympy as sp
import math
r, r0 = sp.symbols('r r0')
r0_value = 10
f_r = r ** 2 * r ** (-3.5)
x = sp.symbols('x')
r_expr = x * r0
f_x = f_r.subs(r, r_expr)
f_x_simplified = sp.simplify(f_x)
r_lower = 10
r_upper = 1000
x_lower = r_lower / r0_value
x_upper = r_upper / r0_value
integral_result = sp.integrate(f_x_simplified, (x, x_lower, x_upper))
numeric_result = integral_result.subs(r0, 10).evalf()
print("Численный результат:", numeric_result)
Последнее редактирование: