Текст подпрограммы и версий
zp42r_p.zip  zp42e_p.zip 
Тексты тестовых примеров
tzp42r_p.zip  tzp42e_p.zip 

Подпрограмма:  ZP42R (модуль ZP42R_p)

Назначение

Отыскание коэффициентов разложения заданного многочлена по многочленам Чебышева.

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

Пусть известны коэффициенты А = ( a1, a2, ..., aN ) разложения многочлена  Р (Х) по степеням  Х

                          N  
            P(X)  =  ∑    ak X k -1
                         k =1 

Подпрограмма вычисляет коэффициенты В = ( b1, b2, ..., bN ) разложения этого многочлена по многочленам Чебышева { Тk (Х) } :

                          N   
            P(X)  =  ∑    bk Tk -1(X)
                         k =1 

Использование

procedure ZP42R(var A :Array of Real; var B :Array of Real;
                var C :Array of Real; N :Integer);

Параметры

A - вектор заданных коэффициентов разложения многочлена по степеням аргумента (тип: вещественный);
B - вектор искомых коэффициентов разложения многочлена по многочленам Чебышева (тип: вещественный);
C - рабочий вектор размерности N на (N + 1)/2 (тип: вещественный);
N - размерность векторов  А и  В (тип: целый).

Версии

ZP42E - отыскание коэффициентов разложения заданного многочлена по многочленам Чебышева для данных расширенной (Extended) точности.

Вызываемые подпрограммы: нет

Замечания по использованию: нет

Пример использования

Unit TZP42R_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, ZP42R_p;

function TZP42R: String;

implementation

function TZP42R: String;
var
N,_i :Integer;
A :Array [0..2] of Real;
B :Array [0..2] of Real;
C :Array [0..5] of Real;
begin
Result := '';  { результат функции }
{ прототип оператора DАТА на FORTRANе  }
A[0] := 0.0;
A[1] := 1.0;
A[2] := 2.0;
N := 3;
ZP42R(A,B,C,N);
Result := Result + #$0D#$0A;
for _i:=0 to 2 do
 begin
  Result := Result + Format('%20.16f ',[B[_i]]);
  if ( ((_i+1) mod 3)=0 )
   then Result := Result + #$0D#$0A;
 end;
Result := Result + #$0D#$0A;
UtRes('TZP42R',Result);  { вывод результатов в файл TZP42R.res }
exit;
end;

end.

Результат:    В = ( 1., 1., 1. )