Замена цикла repeat на цикл while

Program lab1;
Uses Crt;
Const
N=21; Time=600; Lx=0.3; Lyam=0.9; p=2000; Cp=920; T0=50; Te1=30; Te2=70; k1=1000;
k2=500; Ebs2=0.2; Ebs1=0.5; Sigma=5.669e-8; delta=0.00001;
Type vec = array [1..N] of real;
Var
alfa, beta, x,T, Tn1:vec;
A,B,C,F,tau, h,ta, TN, d : real;
i:integer;
Begin
ClrScr;
h:=Lx/(N-1);
x[1]:=0;
For i:=2 to N do
x[i]:=x[i-1]+h;
tau:=Time/1000;
For i:=1 to N do
T[i]:=T0;
While (ta<time) do begin
ta:=ta+tau;
alfa[1]:=2*tau*lyam/(2*tau*(lyam+h*k1)+p*Cp*h*h);
for i:=1 to N do
Tn1[i]:=T[i];
repeat (*)
d:=T[1];
beta[1]:=(p*Cp*h*h*Tn1[1]+2*tau*k1*h*Te1+2*tau*Ebs1*Sigma*h*(sqr(sqr(Te1))-sqr(sqr(d))))/(2*tau*(lyam+h*k1)+p*Cp*h*h);
For i:=2 to N-1 do begin
A:=lyam/sqr(h);
B:=(2*lyam/sqr(h))+p*Cp/tau;
C:=lyam/sqr(h);
F:=-p*Cp*Tn1[i]/tau;
alfa[i]:=A/(B-C*alfa[i-1]);
beta[i]:=(c*beta[i-1]-F)/(B-C*alfa[i-1]);
End;
Repeat
TN:=T[N];
//T[N]:=(p*Cp*sqr(h)*Tn1[N]+2*tau*(lyam*beta[N-1]+k2*h*Te2+Ebs2*Sigma*h*(sqr(sqr(Te2))-sqr(sqr(TN)))))/(p*Cp*sqr(h)+2*tau*(lyam*(1-alfa[N-1])+k2*h));
T[N]:=(Lyam*beta[N-1]+h*k2*Te2+Ebs2*Sigma*h*(Sqr(Te2)*Sqr(Te2)-Sqr(TN)*Sqr(TN)))/(h*k2+Lyam*(1-alfa[N-1]));
Until abs(T[N]-TN)<delta;
For i:= N-1 downto 1 do
T[i]:=alfa[i]*T[i+1]+beta[i];
until abs(d-T[1])<=delta; (*)
End;
For i:=1 to N do
Writeln (x[i]:3:3,' ', T[i]:8:6);
End.