Текст подпрограммы и версий ( Фортран ) zf16r.zip zf16d.zip |
Тексты тестовых примеров ( Фортран ) tzf16r.zip tzf16d.zip |
Текст подпрограммы и версий ( Си ) zf16r_c.zip zf16d_c.zip |
Тексты тестовых примеров ( Си ) tzf16r_c.zip tzf16d_c.zip |
Текст подпрограммы и версий ( Паскаль ) zf16r_p.zip zf16e_p.zip |
Тексты тестовых примеров ( Паскаль ) tzf16r_p.zip tzf16e_p.zip |
Вычисление нуля вещественной функции (трансцендентного уравнения) на заданном интервале методом дихотомии.
ZF16R вычисляет нуль вещественной функции y = F (x) на интервале (A, B), на котором F (x) меняет знак, с заданной точностью EPS методом дихотомии.
И.С.Березин, Н.П.Жидков. Методы вычислений, т.1. Изд - во "Наука", 1966.
SUBROUTINE ZF16R (F, A, B, EPS, ROOT, ITMAX, IERR)
Параметры
F - | имя вещественной подпрограммы - функции вычисления F (x) в любой точке интервала (A, B); |
A, B - | заданные нижняя и верхняя границы интервала, на котором F (x) меняет знак (тип: вещественный); |
EPS - | заданная абсолютная погрешность вычисления нуля функции (тип: вещественный); |
ROOT - | вещественная переменная, значение которой в результате работы подпрограммы полагается равным вычисленному нулю функции F (x); |
ITMAX - | целая переменная, значение которой перед началом работы подпрограммы должно быть положено равным максимальному числу операций, ориентировочно требуемых для обеспечения сходимости к нулю функции; |
IERR - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы подпрограммы; при этом |
IERR=65 - | когда нуль функции не может быть посчитан в пределах заданного числа итераций; |
IERR=66 - | когда функция не меняет знака на заданном интервале |
Версии
ZF16D - | вычисление нуля вещественной функции на заданном интервале методом дихотомии в режиме вычислений с удвоенной точностью. При этом параметры A, B, EPS и ROOT должны иметь тип DOUBLE PRECISION, а подпрограмма - функция вычисления F (x) должна быть оформлена как DOUBLE PRECISION FUNCTION |
Вызываемые подпрограммы
UTZF10 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы ZF16R |
UTZF11 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы ZF16D |
Замечания по использованию
В подпрограммах ZF16R и ZF16D имеется общий блок COMMON /ZF16RR/ ITER. Переменная ITER полагается равной количеству итераций, выполненных при вычислении нуля функции F (x). |
FUNCTION F(X) F = SIN(X) RETURN END EXTERNAL F COMMON /ZF16RR/ ITER A = - 2.57 B = 1.57 EPS = 1.E - 5 ITMAX = 30 CALL ZF16R (F, A, B, EPS, ROOT, ITMAX, IERR) Результаты: ROOT = - 0.000006386669 ITER = 19 IERR = 0