Не программа хочит считать не пойму проблему

unit Unit1;

interface

uses
Winapi.Windows, Winapi.Messages, System.SysUtils, System.Variants, System.Classes, Vcl.Graphics,
Vcl.Controls, Vcl.Forms, Vcl.Dialogs, Vcl.StdCtrls, Vcl.Buttons;

type
TForm1 = class(TForm)
BitBtn1: TBitBtn;
Memo1: TMemo;
Label1: TLabel;
procedure BitBtn1Click(Sender: TObject);
private
{ Private declarations }
public
{ Public declarations }
end;
const
T0=293;
P0=100000;
N: integer=1000;
k: real=0.0012;
r: real=992;
Eps=0.1;
p01=9800 ; t01=293 ; w1=70 ; p02=5000 ; t02=353 ;
w2=100; d=0.3 ; l=0.2 ; D1=0.6 ; D2=0.9 ; t=0.5 ;
var
Form1: TForm1;
i, j: integer;
TAU,dt,P1, P2, T1, T2, temp, re, lyamda, ksi1, ksi2, ksi,
E, v1, v2, nyu, M1, M2, dmas1, dmas2, dm:real;
implementation

{$R *.dfm}

procedure TForm1.BitBtn1Click(Sender: TObject);
begin
i:=1;
j:=1;
dt:=t/n;
TAU:=0;
M1:=0;
M2:=0;
while i<=n do begin
TAU:=TAU+dt;
P1:=p01*cos(w1*TAU)+p0;
T1:=(T01)*cos(w1*TAU)+T0;
P2:=P02*cos(w2*TAU)+P0;
T2:=(T02)*cos(w2*TAU)+T0;
if (P1-P2>0) then begin
temp:=T1;
j:=1;
end
else begin
temp:=T2;
j:=2;
end;
nyu:=0.0178/(1+(0.0337*temp)+0.000221*temp*temp);
repeat
v1:=0;
v1:=v1+0.1;
re:=(v1*d)/nyu;
if (re<=2320) then lyamda:=64/re;
if (re>=3000) and (re<=20*d/k) then lyamda:=0.3164/(exp(0.25*ln(re)));
if (re>=20*d/k) and (re<=500*d/k) then lyamda:=0.11*exp(0.25*ln(k/d)+(68/re));
if (re>=500*d/k) then lyamda:=0.11*exp(0.25*ln(k/d));
If j=1 then begin
ksi2:=((d*d/(d2*d2))-1)*((d*d/(d2*d2))-1);
E:=0.57+0.043/(1.1-(d*d/(d1*d1)));
ksi1:=exp(2*ln((1-E)/E));
end;
If j=2 then begin
ksi2:=((d*d/(d1*d1))-1)*((d*d/(d1*d1))-1);
E:=0.57+0.043/(1.1-(d*d/(d2*d2)));
ksi1:=exp(2*ln((1-E)/E));
end;
ksi:=ksi1+ksi2;
v2:=sqrt(abs((2*(P1-P2))/(r*(((lyamda*l/d)+ksi)))));
until (abs((v1-v2)))M2 then label1.caption:='жтдкость течёт с право на лево'
else label1.caption:='жидкость течёт с лева на право' ;
dm:=abs(M1-M2);
memo1.Lines.Add('M1='+FloatToStr( M1) +#10+#13+ 'M2='+FloatToStr( M2) +#10+#13+'dm='+ FloatToStr(dm));
end;
end.