Текст подпрограммы и версий ( Фортран )
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 

Подпрограмма:  ZP32R

Назначение

Вычисление коэффициентов полинома, проходящего через заданные точки плоскости.

Математическое описание

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