Вывести на экран таблицу значений функции y(x) и ее разложения в ряд

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Classes, Graphics, Controls, Forms, Dialogs,
  StdCtrls, Buttons, ExtCtrls,MatFu;
 
 
type
  TForm1 = class(TForm)
    Memo1: TMemo;
    RadioGroup1: TRadioGroup;
    BitBtn1: TBitBtn;
    BitBtn2: TBitBtn;
    RadioButton4: TRadioButton;
    RadioButton5: TRadioButton;
    procedure FormCreate(Sender: TObject);
    procedure BitBtn1Click(Sender: TObject);
    procedure BitBtn2Click(Sender: TObject);
    private
    { Private declarations }
  public
    { Public declarations }
  end;
 Type
   fun = function(x:extended):extended; // Объявление типа функция
var
  Form1: TForm1;
  xn,xk,h,s:extended;
 
implementation
 
{$R *.DFM}
 
procedure TForm1.BitBtn2Click(Sender: TObject);
begin
Close;
end;
 
procedure TForm1.FormCreate(Sender: TObject);
begin
  xn:=-2;
  xk:=-0.1;
  h:=-0.19;
  Memo1.Clear;
  RadioGroup1.ItemIndex:=0;
end;
 
 
procedure Tabl( f:fun;xn,xk,h:extended);  // Расчет таблицы
      var x,y: extended;
  begin
   x:=xn;
   repeat
     S:=f(x);
     Y:=f(x);
    Form1.Memo1.Lines.Add('x='+FloatToStrf(x,fffixed,8,3)+
                         ' y='+FloatToStrf(y,fffixed,8,3));
        x:=x+h;
   until (x>xk);
 
end;
 
procedure TForm1.BitBtn1Click(Sender: TObject);
var xn,xk,h : extended;
begin
 xn:=-2;
  xk:=-0.1;
  h:=-0.19;
   case RadioGroup1.ItemIndex of   // Выбор функции
   0 : Tabl(Y,xn,xk,h);
    1 : Tabl(S,xn,xk,h);
    end;
end;
 
end.
____________________________________________
unit Matfu;
interface
 
 Function Y (x:extended) : extended;
 Function S(x:extended ) : extended;
 Function Power (x:extended;n:integer): extended;
implementation
 
  Function Y(x:extended )   : extended;
 begin
 Result:=ln(1/(2+2*x+sqr(x)));
 end;
 
Function S(x:extended) : extended;
 var
  S:extended;
  I:Integer;
begin
  S:=0;
  for I := 0 to 16 do
 
    if i mod 2=0 then  S:=S+Power(x,2*i)/i
    else  S:=S+(-1)*(Power(x,2*i)/i);
 
 
  Result:=S;
end;
 
 Function Power;
  begin
    Result:=exp(n*ln(x));
  end;
 
end.