Научная литература
booksshare.net -> Добавить материал -> Математика -> Боглаев Ю.П. -> "Вычислительная математика и программирование " -> 141

Вычислительная математика и программирование - Боглаев Ю.П.

Боглаев Ю.П. Вычислительная математика и программирование — Высшая школа, 1990. — 546 c.
ISBN 5-06-00623-9
Скачать (прямая ссылка): vychmatiprog1990.djvu
Предыдущая << 1 .. 135 136 137 138 139 140 < 141 > 142 143 144 145 146 147 .. 168 >> Следующая

450 4,1 .
1) REAL X,Y 2) DOUBLE PRECISION X,Y 3) REAL X X=0.4 X= 1.0+D—03 INTEGER К
Y=EXP(X) Y=DEXP(X) X = 3. 4
К = INT(X)
12.3.2. Вычисление специальных функций. Определение специальных функций и формулы для их вычисления можно найти в [22,23].
А1А0: интегральная показательная функция Ei(x) Программа вычисляет
00
Ei (х) = J (е “*/s) ds, л: > О,
X
FUNCTION А1А0(ХД)
Параметры входные:
X—вещественный аргумент функции,
I —целое число—индекс ошибки, до обращения следует положить 1 = 0;
выходные:
А1А0—значение Ei(x);
О—нет ошибок,
* Л 1 —аргумент А"<0, для которого функция Е1(х) не определена.
А1А1: гамма-функция Г(л:)
Программа вычисляет
г(*)=
J (е s/sx l)ds, д;>0,
п\пх~1 г
lim —7—-Т—;-----------г, любые х.
п-*О0х(х+\)...(х+п-\)
FUNCTION А1А1(Х,1)
Параметры входные:
X—вещественный аргумент функции,
I —индекс ошибки, до обращения положить 1 = 0;
выходные:
А1А1:—значение Г(х);
15*
451
1=1
О — нет ошибок,
1—аргумент слишком велик,
2—аргумент слишком большое отрицательное число, 3 — аргумент близок к нулю,
4—аргумент — целое отрицательное число.
А1А2: функция ошибок erf(xr) Программа вычисляет
erf(*)=J7=
ехр( — s2)ds,
FUNCTION A1A2(X,I)
Параметры входные:
X—вещественный аргумент функции,
I —индекс ошибки, до обращения положить 1 = 0;
выходные:
А1А2—значение erf(xr);
1 = 0 —нет ошибок.
А1АЗ: функция Бесселя Т0(д:)
Программа вычисляет FqM Ддя *>0
FUNCTION А1АЗ(Х,1)
Параметры входные:
X—вещественный аргумент функции,
I —индекс ошибки, до обращения положить 1 = 0;
выходные:
А1 АЗ: — значение Y0 (х);
Г 0 — нет ошибок,
1=1 1 —аргумент слишком велик,
^2—аргумент отрицателен, ГоМ не определена.
А1А4: функция Бесселя ^(х)
Программа вычисляет Ti(x) для д;>0.
FUNCTION А1А4(Х,1)
Параметры входные:
X—вещественный аргумент функции,
I —индекс ошибки, до обращения положить 1 = 0;
452 ' *
l=<
выходные:
A1A4: — значение yi(x);
0 — нет ошибок,
1 —аргумент слишком велик,
2 — аргумент отрицателен, yt (л;) не определена,
3 — аргумент близок к нулю.
А1А5: функция Бесселя /0(х)
Программа вычисляет /0(^) Ддя х^О, в интервале jc<0 следует использовать тождество J0( — x) = J0(x).
FUNCTION A1A5(X,I)
Параметры входные:
X—вещественный аргумент функции,
I —индекс ошибки, до обращения положить 1 = 0;
выходные:
А1А5:—значение /о(*)»
J_ J 0 — нет ошибок,
~~ I 1 —аргумент слишком велик.
А1А6: функция Бесселя /х(х)
Программа вычисляет ^(х) для л:^0, в интервале д;<0 следует использовать тождество Jx ( — л;) = — ^(х).
FUNCTION А1А6(Х,1)
Параметры входные:
X—вещественный аргумент функции,
I —индекс ошибки, до обращения положить 1 = 0;
выходные:
1 =
А1А6:—значение
0—нет ошибок,
1 —аргумент слишком велик.
А1А7: функция Бесселя К0(х) Программа вычисляет ^(д;) для jc>0.
FUNCTION A1A7(X,I)
453
Параметры входные:
X—вещественный аргумент функции,
I —индекс ошибки, до обращения положить 1 = 0;
выходные:
А1А7: — значение /Co(x);
J I 0 — нет ошибок,
I 1—аргумент меньше нуля,/С0(х)—не определена.
? v) А1А8: функция Бесселя Кх (х)
Программа вычисляет Кх (х) для х>0.
‘ . FUNCTION А1А8(Х,1)
Параметры входные:
X—вещественный аргумент функции,
I —индекс ошибки, до обращения положить 1 = 0;
выходные:
А1А8:—значение i<fi(x);
{0 — нет ошибок,
1 —аргумент меньше нуля, Кх (х)—не определена,
2 — аргумент близок к нулю.
А1А9: функция Бесселя /0(х).
Программа вычисляет /0(х) для х^О, в интервале х<0 следует использовать тождество /0( —х) = /0(х).
FUNCTION А1А9(Х,1)
Параметры входные:
X—вещественный аргумент функции,
I —индекс ошибки, до обращения положить 1 = 0;
выходные:
А1А9: — значение /0(х)
J J 0 — нет ошибок,
~ I 1 — аргумент слишком велик.
454
А1В0: функция Бесселя 1Х (х).
\
'-V
Программа вычисляет /i(x) для дг^О, в интервале jc<0 следует использовать тождество /х (— л;) = — /х (л;).
FUNCTION А1В0(Х,1)
Параметры входные:
X—вещественный аргумент функции,
I —индекс ошибки, до обращения положить 1 = 0;
выходные:
1 =
А1В0: — значение /i(x)
0 — нет ошибок,
1 — аргумент слишком велик.
12.3.3. Вычисление элементарных функций комплексного переменного. Программы А2А0 — А2А8 являются библиотечными функциями фортрана, обращение к которым осуществляется с соответствующим описанием типа аргумента и функции, например:
COMPLEX Z,W REAL W
Z = (1.7, —3.2) COMPLEX Z
W = CCOS(Z) Z = (1.7, — 3.2)
W = AIMAG(Z)
12.3.4. Аппроксимация функций
АЗ АО: одномерная интерполяция по схеме Эйткена
Программа вычисляет по таблично заданной функции У1=/(х() в узлах xh\^i^N, значение функции /(z) в точке гфх{ по схеме Эйткена с заданной точностью в.
SUBROUTINE A3A0(N,X,Y,Z,E,F,P,I)
Параметры входные:
N — количество узлов,
X—одномерный массив хь 1 < / < N,
Y—одномерный массив yi9
Z—значение z, где функция вычисляется интерполированием,
Е—заданная точность в,
F — одномерный рабочий массив Fi9 К i ^ N;
выходные:
Р—значение /(z),
I —индекс ошибки,
Го — нет ошибок,
1=4 1 —количество узлов N недостаточно для достижения заданной точности.
455
АЗА1: двумерная интерполяция бикубическими сплайнами
Предыдущая << 1 .. 135 136 137 138 139 140 < 141 > 142 143 144 145 146 147 .. 168 >> Следующая

Реклама

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed

Есть, чем поделиться? Отправьте
материал
нам
Авторские права © 2009 BooksShare.
Все права защищены.
Rambler's Top100

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed