График в делфи (аппроксимация. метод наименьших квадратов)

unit Unit1;
interface
uses
  Windows, Messages, SysUtils, Classes, Graphics, Controls, Forms, Dialogs,
  StdCtrls, TeEngine, Series, ExtCtrls, TeeProcs, Chart;
type
  TForm1 = class(TForm)
    GroupBox1: TGroupBox;
    Label1: TLabel;
    Label2: TLabel;
    Memo_ish: TMemo;
    B_open: TButton;
    GroupBox2: TGroupBox;
    Memo_rez: TMemo;
    B_View: TButton;
    B_Save: TButton;
    OpenDialog1: TOpenDialog;
    SaveDialog1: TSaveDialog;
    OpenDialog2: TOpenDialog;
    B_close: TButton;
    B_clear: TButton;
    B_calc: TButton;
    Chart1: TChart;
    Series1: TFastLineSeries;
    Series2: TFastLineSeries;
 
 
 
    procedure B_openClick(Sender: TObject);
    procedure B_SaveClick(Sender: TObject);
    procedure B_viewClick(Sender: TObject);
    procedure B_closeClick(Sender: TObject);
    procedure B_clearClick(Sender: TObject);
    procedure B_calcClick(Sender: TObject);
  private
  public
  end;
 
var
  Form1: TForm1;
 
 
Implementation
 
{$R *.dfm}
 
 Const n=9;
 
    Var
 
       x,yp:array [1..10] of real;
        S_x3y,S_x3,S_x6,S_y,S_x,S_x2,S_x4,S_xy: real;
       d,d1,d2,d3:real;
        FName_ish:string;
        FName_rez:string;
        File_ish,File_rez:TextFile;
        A,C,B:real;
 
 
procedure TForm1.B_openClick(Sender: TObject);
 
    Var i:integer;
  begin
  if OpenDialog1.Execute then
     begin
     FName_ish:=OpenDialog1.FileName;
     Memo_ish.Lines.LoadFromFile(FName_ish);
     AssignFile(File_ish,FName_ish);
     Reset(File_ish);
       for i:=1 to n do
         read(File_ish,x[i],Yp[i]);
         readln(File_ish);
      end;
 
 
  end;
procedure TForm1.B_SaveClick(Sender: TObject);
   Var i:integer;
 begin
   if SaveDialog1.Execute then
      begin
       FName_rez:=SaveDialog1.FileName;
           AssignFile(File_rez,FName_rez);
           Rewrite(File_rez);
           writeln(File_rez,'Теоретические значения y:');
           for i:=1 to n do
           writeln (File_rez,Yp[i]:5:3);
           writeln(File_rez);
           writeln (File_rez,'========================');
 
 
           writeln(File_rez);
           writeln(File_rez,'========================');
           writeln(File_rez,'Искомое уравнение:');
           writeln(File_rez,'y='+FormatFloat('0.00', A)+' x*x*x +'
           +FormatFloat('0.00',C)+'x+'+FormatFloat('0.00',B));
           Memo_rez.Lines.Add('y='+FormatFloat('0.00',A)
          +' x*x*x+ '+FormatFloat('0.00',C)+'x+'+FormatFloat('0.00',B));
           CloseFile(File_rez);
        end;
  end;
 
 
procedure TForm1.B_closeClick(Sender: TObject);
 begin
  close;
 end;
 
procedure TForm1.B_clearClick(Sender: TObject);
  begin
  Memo_ish.clear;
  Memo_rez.clear;
  end;
procedure TForm1.B_ViewClick(Sender: TObject);
    begin
     if OpenDialog2.Execute then
        begin
         FName_ish:=OpenDialog2.FileName;
         Memo_rez.Lines.LoadFromFile(FName_ish);
        end;
    end;
 
procedure TForm1.B_calcClick(Sender: TObject);
var i: integer;
 begin
 
S_x:=0;
S_x3y:=0;
S_x3:=0;
S_x6:=0;
S_y:=0;
S_x2:=0;
S_x4:=0;
S_xy:=0;
for i:=1 to n do
 
  begin
S_x3y:=S_x3y+x[i]*x[i]*x[i]*yp[i];
S_x3:=S_x3+x[i]*x[i]*x[i];
S_x6:=S_x6+x[i]*x[i]*x[i]*x[i]*x[i]*x[i];
S_y:=S_y+yp[i];
S_x:=S_x+x[i];
S_xy:=S_xy+x[i]*yp[i];
S_x2:=S_x2+x[i]*x[i];
S_x4:=S_x4+x[i]*x[i]*x[i]*x[i];
  end;
 
d:=S_x6*(S_x2*n-S_x*S_x)-S_x4*(S_x4*n-S_x3*S_x)+S_x3*(S_x4*S_x-S_x3*S_x2);
d1:=S_x3y*(S_x2*n-S_x*S_x)-S_x4*(S_xy*n-S_y*S_x)+S_x3*(S_xy*S_x-S_y*S_x2);
d2:=S_x6*(S_xy*n-S_y*S_x)-S_x3y*(S_x4*n-S_x3*S_x)+S_x3*(S_x4*S_y-S_x3*S_xy);
d3:=S_x6*(S_x2*S_y-S_x*S_xy)-S_x4*(S_x4*S_y-S_x3*S_xy)+S_x3y*(S_x4*S_x-S_x3*S_x2);
A:=d1/d;
C:=d2/d;
B:=d3/d;
 
 
   for i:=1 to n do
   begin
yp[i]:=A*x[i]*x[i]*x[i]+C*x[i]+B;
 
 
   end;
end;
 
end.