Выведение массива. Построение графика. Паскаль. Делфи.

[Rabbit777]

Не могу справиться с ошибками при выводе массива. Необходимо также построить график по двум конечным величинам(широта и долгота). Код ниже.

unit uRabbit777;

interface

uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls;

type
  TForm1 = class(TForm)
    Memo1: TMemo;
    Memo2: TMemo;
    Button1: TButton;
    Label1: TLabel;
    GroupBox1: TGroupBox;
    Edit1: TEdit;
    Label2: TLabel;
    Label3: TLabel;
    Edit2: TEdit;
    Label4: TLabel;
    GroupBox2: TGroupBox;
    Edit3: TEdit;
    Label5: TLabel;
    Edit4: TEdit;
    Label6: TLabel;
    Label7: TLabel;
    Label8: TLabel;
    Edit5: TEdit;
    Edit6: TEdit;
    Edit7: TEdit;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  Form1: TForm1;
const
  K_GAUSS=0.01720209895; // см - https://ru.wikipedia.org/wiki/Постоянная_Гаусса
  EPS=1e-6;              // погрешность, радианы
implementation
uses Math;

{$R *.dfm}
function MeanMotion(a:Double):Double;  // среднее движение , радиан/сутки
begin
 Result:=K_GAUSS/(a*Sqrt(a));
end;
function KeplerEqu(M,e:Double):Double; // уравнение Кеплера M-->E
var E1,E2:Double;
begin
 E2:=M;
 repeat
  E1:=E2;
  E2:=M+e*Sin(E1);
 until Abs(E1-E2)<EPS;
 Result:=E2;
end;
function TrueAnomV(Ex,e:Double):Double;  // истиная аномалия
begin
 Result:=2*ArcTan(Tan(Ex/2)*Sqrt((1+e)/(1-e)));
end;
function RadVect(a,e,Ex:Double):Double;  // радиус-вектор
begin
 Result:=a*(1-e*Cos(Ex));
end;
//
procedure TForm1.Button1Click(Sender: TObject);
var a,e,t0,n,i,u,w,Node,l_lon,b_lat,Period,t,M,ExAn,v,r,Lon:Double;
    k,z:Integer;
begin
 Label7.Caption:='';
 Memo1.Lines.Clear;
 Memo2.Lines.Clear;
 a:=StrToFloatDef(Edit1.Text,0);
 e:=StrToFloatDef(Edit2.Text,0);
 i:=StrToFloatDef(Edit5.Text,0);
 w:=StrToFloatDef(Edit6.Text,0);
 v:=StrToFloatDef(Edit7.Text,0);
 if (a<0)or(e<0)or(e>=1) then Exit;
 n:=MeanMotion(a);
 Period:=2*PI/n;
 Label7.Caption:='период = '+FloatToStr(RoundTo(Period,-2))+' секунд';
 t0:=StrToFloatDef(Edit3.Text,0);
 k:=StrToIntDef(Edit4.Text,0);
 for z:=0 to k do
  begin
   t:=t0+z*Period/k;
   M:=t*n;
   if M<0 then M:=M+2*PI;
   if M>2*PI then M:=M-2*PI;
   ExAn:=KeplerEqu(M,e);
   u:=w+v;
 l_lon:=Node+ArcTan2(Cos(i)*Sin(u),Cos(u));
 b_lat:=ArcSin(Sin(u)*Sin(i));
   Memo1.Lines.Add(FloatToStr(l_lon));
   Memo2.Lines.Add(FloatToStr(b_lat));
  end;
end;

end.

[type_Oleg]

Это кросспост. Уж третий раз одно и то же.