Текст подпрограммы и версий
zp37r_p.zip  zp37e_p.zip  		Тексты тестовых примеров
tzp37r_p.zip  tzp37e_p.zip 

Подпрограмма:  ZP37R (модуль ZP37R_p)

Назначение

Возведение в степень полинома с вещественными коэффициентами.

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

ZP37R вычисляет М первых коэффициентов полинома 

       B(X)  =  XS*(B(1) + B(2)*X + ... + B(M)*X(M - 1) + ...)  ,

такого, что  В (Х) = (А (Х))R, где  А (Х) - заданный полином, степени NР 

       A(X)  =  A(1) + A(2)*X + A(3)*X2 + ... + A(NP + 1)*XNP  ,

у которого коэффициенты  А ( I ),  I = 1, 2, 3, ..., NР + 1 являются вещественными. Здесь S = К * R, где К - число первых коэффициентов полинома  А (Х), равных 0.

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

procedure ZP37R(NP :Integer; var A :Array of Real; M :Integer;
                var B :Array of Real; R :Real; var S :Real);


   Параметры
NP - заданная степень полинома  А (Х) (тип: целый);
A - вещественный вектор длины NР + 1, содержащий коэффициенты заданного полинома  А (Х) в порядке возрастания степеней;
M - заданное число искомых коэффициентов полинома  В (Х) (тип: целый);
B - вещественный вектор длины М, содержащий вычисленные коэффициенты полинома  В (Х) в порядке возрастания степеней;
R - заданный показатель степени (тип: вещественный);
S - вещественная переменная, значение которой в результате работы подпрограммы полагается равным К * R. 
   Версии
ZP37E - возведение в степень полинома с вещественными коэффициентами с расширенной (Extended) точностью. При этом векторы А и В должны иметь тип Extended. 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию:  нет

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

Unit tzp37r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, ZP37R_p;

function tzp37r: String;

implementation

function tzp37r: String;
var
NP,M,_i :Integer;
R,S :Real;
B :Array [0..3] of Real;
const
A :Array [0..4] of Real = ( 24.0,-50.0,35.0,-10.0,1.0 );
begin
Result := '';  { результат функции }
R := 2.0;
NP := 4;
M := 4;
ZP37R(NP,A,M,B,R,S);
Result := Result + Format('%10d %10d %20.16f %20.16f ',
 [NP,M,R,S]) + #$0D#$0A;
Result := Result + #$0D#$0A;
for _i:=0 to 4 do
 begin
  Result := Result + Format('%20.16f ',[A[_i]]);
  if ( ((_i+1) mod 4)=0 )
   then Result := Result + #$0D#$0A;
 end;
Result := Result + #$0D#$0A;
Result := Result + #$0D#$0A;
for _i:=0 to 3 do
 begin
  Result := Result + Format('%20.16f ',[B[_i]]);
  if ( ((_i+1) mod 4)=0 )
   then Result := Result + #$0D#$0A;
 end;
Result := Result + #$0D#$0A;
UtRes('tzp37r',Result);  { вывод результатов в файл tzp37r.res }
exit;
end;

end.


Результаты:

       В(1) = 576. ,    В(2) = - 2400. ,   В(3) = 4180. ,  B(4)  =  - 3980. ,
       S  =  0