Смещенная синусоидальная интегральная функция
ssinint( возвращает смещенную синусоидальную интегральную функцию X)ssinint(X) = sinint(X) — pi/2.
В зависимости от его аргументов, ssinint возвращает результаты с плавающей запятой или точные символьные результаты.
Вычислите смещенную синусоидальную интегральную функцию для этих чисел. Поскольку эти числа не являются символическими объектами, ssinint возвращает результаты с плавающей запятой.
A = ssinint([- pi, 0, pi/2, pi, 1])
A = -3.4227 -1.5708 -0.2000 0.2811 -0.6247
Вычислите смещенную синусоидальную интегральную функцию для чисел, преобразованных в символьные объекты. Для многих символических (точных) чисел, ssinint возвращает неразрешенные символьные вызовы.
symA = ssinint(sym([- pi, 0, pi/2, pi, 1]))
symA = [ - pi - ssinint(pi), -pi/2, ssinint(pi/2), ssinint(pi), ssinint(1)]
Использовать vpa для аппроксимации символьных результатов числами с плавающей запятой:
vpa(symA)
ans = [ -3.4227333787773627895923750617977,... -1.5707963267948966192313216916398,... -0.20003415864040813916164340325818,... 0.28114072518756955112973167851824,... -0.62471325642771360428996837781657]
Постройте график сдвинутой синусоидальной интегральной функции на интервале от -4*pi кому 4*pi.
syms x fplot(ssinint(x),[-4*pi 4*pi]) grid on

Многие функции, такие как diff, int, и taylor, может обрабатывать выражения, содержащие ssinint.
Найдите первую и вторую производные смещенной синусоидальной интегральной функции:
syms x diff(ssinint(x), x) diff(ssinint(x), x, x)
ans = sin(x)/x ans = cos(x)/x - sin(x)/x^2
Найдите неопределенный интеграл смещенной синусоидальной интегральной функции:
int(ssinint(x), x)
ans = cos(x) + x*ssinint(x)
Найти расширение серии Тейлор ssinint(x):
taylor(ssinint(x), x)
ans = x^5/600 - x^3/18 + x - pi/2
[1] Гаутши, У. и У. Ф. Кэхилл. «Экспоненциальный интеграл и связанные функции». Справочник по математическим функциям с формулами, графиками и математическими таблицами. (М. Абрамовиц и И. А. Стегун, ред.). Нью-Йорк: Дувр, 1972.