Экспоненциальная интегральная функция с одним аргументом
Вычислите экспоненциальные интегралы для числовых входных данных. Поскольку эти числа не являются символическими объектами, получаются результаты с плавающей запятой.
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 δ i при пересечении этого разреза. Вычислите экспоненциальные интегралы в -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] Гаутши, У. и У. Ф. Гейхилл «Экспоненциальный интеграл и связанные функции». Справочник по математическим функциям с формулами, графиками и математическими таблицами. (М. Абрамовиц и И. А. Стегун, ред.). Нью-Йорк: Дувр, 1972.