Помогите разобраться.

PROGRAM zam_capl;
uses crt;
label w1,w2,w3,w4;
const
n=100;            {количество точек}
z1=1.0;z2=0.1; {коэффициент теплопроводности/ плотность/ уд. Теплоёмкость}
s=0.001;s2=0.002; {коэффициент теплопроводности/плотность/уд. Теплота плавления}
t0=0;                   {температура замерзания воды  }
T1=50;lm=0;       {температура и радиус внутренней области капли}
l=100;T2=-50;     {температура и радиус внешней области капли}
VAR
Rl,dR,r,dt,dtm:real;
k,i,nt,nm,j:integer;
T,tt,td,TT1:array [0..n+1] of real ;
RR:array [0..200] of real ;
capl:text;
BEGIN
 assign(capl, 'c:\kapl50.dat');
 rewrite(capl);
 clrscr;
 dtm:=0.001;  dr:=l/n; nm:=trunc(n/l*lm);
FOR k:=1 TO n do begin
   T[k]:=T1; end;
   T[n]:=t0;
   T[n+1]:=T2;
   nt:=n; Rl:=dR*Nt; i:=0;
while i<100 do begin
    dt:=32;
    for k:=nm+1 to nt-1 do begin
       r:=l/n*k;
       td[k]:=z1/(dR*dR)*(T[k-1]*(1-dR/r)+T[k+1]*(1+dr/r)-2*T[k]);
       for j:=1 to 200 do begin
if T[j]+dt*td[j]>T[j-1]+dt*td[j-1] then dt:=dt/2 else goto w1;
if dt=0 then dt:=dtm ;
     w1:end;
  end;
for k:=1 to 200 do begin
  if T[nt-1]+dt*td[nt-1]<T[nt] then dt:=dt/2 else goto w2;
  if dt=0 then dt:=dtm ;
   w2:end;
for k:=nt+1 to n do begin
    r:=l/n*k;
    td[k]:=z2/(dR*dR)*(T[k-1]*(1-dR/r)+T[k+1]*(1+dr/r)-2*T[k]);
    for j:=1 to 500 do begin
    if T[j]+dt*td[j]>T[j-1]+dt*td[j-1] then dt:=dt/2 else goto w3;
         if dt=0 then dt:=dtm;
 w3:end; end;
 for j:=1 to 200 do begin
   if T[nt+1]+dt*td[nt+1]>T[nt] then dt:=dt/2 else goto w4;
   if dt=0 then dt:=dtm ;
w4:end;
   Rl:=Rl-(s*(T[nt]-T[nt-1])-s2*(T[nt+1]-T[nt]))*dt;
   nt:=trunc(Rl/dR);
   RR[J]:=RL;
   if T[nt]<>t0 then T[nt]:=t0;
   td[nt]:=0;
   T[nt+1]:=(T[nt+2]+t0)/2;
   T[nt-1]:=(T[nt-2]+t0)/2;
   for k:=0 to n do
begin
     T[k]:=T[k]+td[k]*dt;
end;
     i:=i+1;
     T[0]:=T[1]
  end;
   for i:=0 to n do begin
     writeln (capl,RR[n],' ',T[i]);
   end; close(capl);
end.