|
Текст подпрограммы и версий ( Фортран ) zp43r.zip zp43d.zip |
Тексты тестовых примеров ( Фортран ) tzp43r.zip tzp43d.zip |
|
Текст подпрограммы и версий ( Си ) zp43r_c.zip |
Тексты тестовых примеров ( Си ) tzp43r_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) zp43r_p.zip zp43e_p.zip |
Тексты тестовых примеров ( Паскаль ) tzp43r_p.zip tzp43e_p.zip |
Вычисление значения многочлена, заданного разложением по ортогональным многочленам.
Вычисляется величина
N
Y = ∑ Ck Qk -1(X) ,
k =1
где Х - значение аргумента, C = ( c1, c2, ..., cN ) заданные коэффициенты разложения многочлена по системе { Qk (x) } ортогональных многочленов, удовлетворяющих трехчленному рекуррентному отношению:
Q0(x) = 1 , Q1(x) = x - a1 ,
Qk +1(x) = ( x - ak +1 ) Qk(x) - bk Qk -1(x) , k = 1, 2, ..., N-2
с известными коэффициентами А = ( a1, a2, ..., aN - 1 ), B = ( b1, b2, ..., bN - 2 ).
SUBROUTINE ZP43R (A, B, C, X, Y, N)
Параметры
| A - | вещественный вектор длины N - 1, содержащий коэффициенты ak рекуррентного соотношения; |
| B - | вещественный вектор длины N - 2, содержащий коэффициенты bk рекуррентного соотношения; |
| C - | вещественный вектор длины N, содержащий коэффициенты разложения многочлена по системе ортогональных многочленов; |
| X - | значение аргумента (тип: вещественный); |
| Y - | вычисляемое значение многочлена (тип: вещественный); |
| N - | число коэффициентов разложения (тип: целый). |
Версии
| ZP43D - | вычисление значения многочлена, заданного разложением по ортогональным многочленам, для данных двойной точности. |
Вызываемые подпрограммы: нет
Замечания по использованию: нет
DIMENSION A(2), B(1), C(3)
DATA A(1), A(2), B(1) /0., 2*1./, C(1), C(2), C(3) /3*1./
X = 1.
N = 3
CALL ZP43R (A, B, C, X, Y, N)
Результат: Y = 1.