procedure coeff;
var
  h1,h2,hn1,hn2,hii,hi,hi1:double;
  ii,i1,i2,n1,n2,n3:integer;
begin
  h1:=x[2]-x[1];
  h2:=x[3]-x[2];
  n1:=n-1;
  n2:=n-2;
  hn1:=x[n]-x[n1];
  hn2:=x[n1]-x[n2];
  c[1]:=h1/3+1/sqr(h1)*(ro[1]+ro[2]);
  d[1]:=h1/6-1/h1*(1/h1+1/h2)*ro[2]-1/sqr(h1)*ro[1];
  e[1]:=ro[2]/(h1*h2);
  b[2]:=d[1];
  a[3]:=e[1];
  a[n]:=ro[n1]/(hn1*hn2);
  b[n]:=hn1/6-1/hn1*(1/hn1+1/hn2)*ro[n1]-ro[n]/sqr(hn1);
  c[n]:=hn1/3+1/sqr(hn1)*(ro[n1]+ro[n]);
  d[n1]:=b[n];
  e[n2]:=a[n];
  f[1]:=(y[2]-y[1])/h1-y1;
  f[n]:=yn-(y[n]-y[n1])/hn1;
  for i:=2 to n1 do
  begin
    II:=I-1;
    i1:=i+1;
    hii:=x[i]-x[ii];
    hi:=x[i1]-x[i];
    f[i]:=((y[i1]-y[i])/hi-(y[i]-y[ii])/hii);
  end;
 
  for i:=2 to n2 do
  begin
    ii:=i-1;
    i1:=i+1;
    i2:=i+2;
    hii:=x[i]-x[ii];
    hi:=x[i1]-x[i];
    hi1:=x[i2]-x[i1];
    d[i]:=hi/6-1/hi*((1/hii+1/hi)*ro[i]+(1/hi+1/hi1)*ro[i1]);
    b[i1]:=d[i];
  end;
  n3:=n-3;
 
  for i:=2 to n3 do
  begin
    i1:=i+1;
    i2:=i+2;
    hi:=x[i1]-x[i];
    hi1:=x[i2]-x[i1];
    e[i]:=ro[i1]/(hi*hi1);
    a[i2]:=e[i];
  end;
...

 
{************************************************************}
{                                                            }
{                             KSpline                        }
{                Copyright (c) 2009  KSoftWare               }
{                                                            }
{                                                            }
{  Разработчик: Николай Товеровский                          }
{  Модифицирован: 31 марта 2009                              }
{  Описание: Этот модуль предоставляет удобный               }
{            и безопасный доступ к методам                   }
{            модуля FSplineUnit, который позволяет           }
{            аппроксимировать функции сплайном.              }
{************************************************************}
unit KSpline;
 
interface
uses
  SysUtils, KArrayTypes, KArrayRoutines;
 
type
 
 {TSplineController — основной класс, инкапсулирующий работу с
  модулем FSplineUnit.}
  TSplineController = class
  private
    function GetApproximatingArray: TDoublePointArray;
    function GetSecondDerivativeArray: TDoublePointArray;
    function GetFirstDerivativeArray: TDoublePointArray;
    function GetInputArray: TDoublePointArray;      
  protected
   {Этот метод позволяет извлекать данные из глобальных массивов FSplineUnit.}
    function GetDataFromGlobalArrays(const GlobalXArray: array of Double;
      const GlobalYArray: array of Double): TDoublePointArray;
  public
   {Применить сплайн к функции.
 
    На вход подается массив и жесткость сплайна.
    Чем выше жесткость, тем более гладкой получится кривая.
 
    В результате получится набор массивов, которые доступны через
    свойства такие как ApproximatingArray (см. ниже).
 
    Заметьте, что жесткость можно задавать для каждой точки отдельно,
    но эта возможность скрыта внутри.}
    procedure ApplySpline(const InputArray: TDoublePointArray;
      Rigidity: Double = 50000);
 
   {Через эти свойства можно получить доступ к результатам работы
    аппроксимации (см. ApplySpline).}
    property InputArray: TDoublePointArray read GetInputArray;
    property ApproximatingArray: TDoublePointArray read GetApproximatingArray;
    property FirstDerivativeArray: TDoublePointArray read GetFirstDerivativeArray;
    property SecondDerivativeArray: TDoublePointArray read GetSecondDerivativeArray;
  end;