Одноаргументный экспоненциал интегральная функция
Вычислите экспоненциальные интегралы для числовых входов. Поскольку эти числа не являются символическими объектами, вы получаете результаты с плавающей точкой.
s = [ei(-2), ei(-1/2), ei(1), ei(sqrt(2))]
s = -0.0489 -0.5598 1.8951 3.0485
Вычислите экспоненциальные интегралы для тех же чисел, преобразованных в символические объекты. Для большинства символических (точных) чисел ei возвращает неразрешенные символические вызовы.
s = [ei(sym(-2)), ei(sym(-1/2)), ei(sym(1)), ei(sqrt(sym(2)))]
s = [ ei(-2), ei(-1/2), ei(1), ei(2^(1/2))]
Использовать vpa чтобы аппроксимировать этот результат с 10-значной точностью.
vpa(s, 10)
ans = [ -0.04890051071, -0.5597735948, 1.895117816, 3.048462479]
Отрицательная действительная ось - это разрез ветви. Экспоненциальный интеграл имеет переход высоты 2 Вычислите экспоненциальные интегралы в -1, выше -1, и ниже -1 чтобы продемонстрировать это.
[ei(-1), ei(-1 + 10^(-10)*i), ei(-1 - 10^(-10)*i)]
ans = -0.2194 + 0.0000i -0.2194 + 3.1416i -0.2194 - 3.1416i
Вычислите первую, вторую и третью производные экспоненциального интеграла с одним аргументом.
syms x diff(ei(x), x) diff(ei(x), x, 2) diff(ei(x), x, 3)
ans = exp(x)/x ans = exp(x)/x - exp(x)/x^2 ans = exp(x)/x - (2*exp(x))/x^2 + (2*exp(x))/x^3
Вычислите пределы экспоненциального интеграла с одним аргументом.
syms x limit(ei(2*x^2/(1+x)), x, -Inf) limit(ei(2*x^2/(1+x)), x, 0) limit(ei(2*x^2/(1+x)), x, Inf)
ans = 0 ans = -Inf ans = Inf
Интеграл экспоненциала одним аргументом сингулярен в x = 0. Тулбокс использует это специальное значение: ei(0) = -Inf.
Отношение между ei и expint является
ei(x) = -expint(1,-x) + (ln(x)-ln(1/x))/2 - ln(-x)
Обе функции ei(x) и expint(1,x) имеют логарифмическую особенность в источник и разрез ветви вдоль отрицательной действительной оси. ei функция не непрерывна при приближении сверху или ниже этого разреза ветви.
[1] Gautschi, W. and W. F. Gahill «Exponental Integral and Related Function». Руководство по математическим функциям с формулами, графиками и математическими таблицами. (М. Абрамовиц и И. А. Штегун, эд.). Нью-Йорк: Дувр, 1972.