График что та не то

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, TeEngine, Series, ExtCtrls, TeeProcs, Chart;
  const mf=500;
type
vector=array[1..mf] of real;
  TForm1 = class(TForm)
    Button1: TButton;
    Memo1: TMemo;
    Edit1: TEdit;
    Label1: TLabel;
    Edit2: TEdit;
    Label2: TLabel;
    Edit3: TEdit;
    Label3: TLabel;
    Label4: TLabel;
    Edit5: TEdit;
    Label5: TLabel;
    Edit6: TEdit;
    Label6: TLabel;
    Edit7: TEdit;
    Label7: TLabel;
    Edit8: TEdit;
    Label8: TLabel;
    Chart1: TChart;
    Series1: TLineSeries;
    Edit4: TEdit;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
 
implementation
 
{$R *.dfm}
function lamda(x:real):real;
{Функция вычисления коэффициента теплопроводности по формуле (31)}
begin
lamda:=5500/(560+x)+0.942*(1E-10)*x*sqr(x);
end;
procedure TForm1.Button1Click(Sender: TObject);
var
{Раздел описания переменных, которые мы будем использовать в программе}
  i,j,N:integer;
  T,alfa,beta:vector;
  ai,bi,ci,fi:real;
  ro,ch,tau,h,kap:real;
  Th,T0,Tc,L,t_end,time:real;
//f,g:text;
begin
        N:=StrToInt(Edit1.Text);
        t_end:=StrToFloat(Edit2.Text);
        L:=StrToFloat(Edit3.Text);
        ro:=StrToFloat(Edit4.Text);
        kap:=StrToFloat(Edit5.Text);
        T0:=StrToFloat(Edit6.Text);
        Th:=StrToFloat(Edit7.Text);
        Tc:=StrToFloat(Edit8.Text);
{определяем расчетный шаг сетки по пространственной координате}
h:=L/(N-1);
{определяем расчетный г=шаг сетки по времени}
tau:=t_end/100;
{определяем поле температуры в начальный момент времени}
for i:=1 to N do
T[i]:=T0;
{проводим интегрирование нестациораного уравнения теплопровдности}
for i:=1 to N do
T[i]:=T0;
time:=0;
Chart1.Series[0].Clear;
//Chart1.Series[1].Clear;
//Chart1.Series[2].Clear;
//Chart1.Series[3].Clear;
while time<t_end do {используем цикл с предусловием}
begin
{увеличиваем переменную времени шаг tau}
time:=time+tau;
{определяем начальные прогоночные коэффициенты на основе левого граничного условия}
alfa[1]:=0;
beta[1]:=Th;
{цикл с параметром для определения прогончных коэффициентов по формуле (8)}
for i:=2 to N-1 do
begin
{ai, bi, ci ,fi - коэффициенты канонического представления системы уравнений с трехдиагнальной матрицей}
ai:=0.5*(lamda(T[i])+lamda(T[i+1]))/sqr(h);
ci:=0.5*(lamda(T[i-1])+lamda(T[i]))/sqr(h);
bi:=ai+ci+ro*kap/tau;
fi:=-ro*kap*T[i]/tau;
{alfa[i], beta[i] - прогоночные коэффициенты}
alfa[i]:=ai/(bi-ci*alfa[i-1]);
beta[i]:=(ci*beta[i-1]-fi)/(bi-ci*alfa[i-1]);
end;
{определяем значение температуры на правой границе на основе правого граничного условия}
T[N]:=Tc;
{используя соотношения (7) определяем неизвестное поле температуры}
for i:=N-1 downto 1 do 
T[i]:=alfa[i]*T[i+1]+beta[i];
 Memo1.Lines.Clear;
  for i:=1 to n do
   begin
     Chart1.Series[0].Add(T[i],'',clred);
      Memo1.Lines.Add(FloatToStr(T[i]));
 end;
end; {цикл с предусловием окончен}
{выводим результат в файл}
end;
 
end.