Текст подпрограммы и версий ( Фортран ) zp32r.zip zp32d.zip |
Тексты тестовых примеров ( Фортран ) tzp32r.zip tzp32d.zip |
Текст подпрограммы и версий ( Си ) zp32r_c.zip zp32d_c.zip |
Тексты тестовых примеров ( Си ) tzp32r_c.zip tzp32d_c.zip |
Текст подпрограммы и версий ( Паскаль ) zp32r_p.zip zp32e_p.zip |
Тексты тестовых примеров ( Паскаль ) tzp32r_p.zip tzp32e_p.zip |
Вычисление коэффициентов полинома, проходящего через заданные точки плоскости.
ZР32R вычисляет коэффициенты полинома степени NР
A(1) + A(2)*X + A(3)*X2 + ... + A(NP + 1)*XNP ,
проходящего через заданные точки плоскости (Х ( I ), Y ( I )), I = 1, 2, 3, ..., NР + 1.
Другими словами, коэффициенты полинома вычисляются по заданным соотношениям:
NP+1 Y( I ) = ∑ A( J )*X( I )(J - 1) , I = 1, 2, 3, ..., NP+1 J=1
SUBROUTINE ZP32R (NP, X, Y, A, RAB1, RAB2)
Параметры
NP - | заданная степень полинома (тип: целый); |
X, Y - | вещественные векторы длины NР + 1, содержащие координаты заданных точек плоскости; |
A - | вещественный вектор длины NР + 1, содержащий вычисленные коэффициенты полинома в порядке возрастания степеней; |
RAB1 - RAB2 | вещественные векторы длины NР + 1, используемые в подпрограмме как рабочие. |
Версии
ZP32D - | вычисление коэффициентов полинома, проходящего через заданные точки плоскости, с повышенной точностью. Параметры Х, Y, А, RАВ1, RАВ2 имеют тип DОUВLЕ РRЕСISIОN. |
Вызываемые подпрограммы
ZP30R - | вычисление значения полинома с вещественными коэффициентами от заданного вещественного аргумента. Используется в подпрограмме ZР32R. |
ZP31R - | вычисление коэффициентов полинома по заданным его вещественным корням. Используется в подпрограмме ZР32R. |
ZP30D - | вычисление значения полинома с вещественными коэффициентами от заданного вещественного аргумента с повышенной точностью. Используется в подпрограмме ZР32D. |
ZP31D - | вычисление коэффициентов полинома по заданным его вещественным корням с повышенной точностью. Используется в подпрограмме ZР32D. |
Замечания по использованию: нет
DIMENSION X(5), Y(5), A(5), RAB1(5), RAB2(5) DATA X(1), X(2), X(3), X(4), X(5) /0., 1., 2., 3., 4./ DATA Y(1), Y(2), Y(3), Y(4), Y(5) /24., 0., 0., 0., 0./ NP = 4 CALL ZP32R (NP, X, Y, RAB1, RAB2) Результаты: А(1) = 24.0 , А(2) = - 50.0 , А(3) = 35.0 , A(4) = - 10.0 , A(5) = 1.0