Программа для решения физической задачи

begin
  e0:=8.85*exp(ln(10)*(-12));
  C:=e1*e0*R*(((pi*R)/L)+ln((16*pi*R)/L)-1);
  while L/R<1 do begin
    razn:=Abs(R-L);
    if razn<prev_razn then begin
      mem:=R;
      prev_razn:=razn;
    end;
    with Form1.Series1 do AddXY(R,L);
    R:=R+0.0001;
  end;
  edt1.Text:=FloatToStrF(mem,ffFixed,4,2);
end;

 with Form1.Series1 do AddXY(R,L);

begin
**Series1.Clear;
**Re:=0;
**Reh:=0.0012;
**e0:=8.85*exp(Ln(10)*(-12));
**for i:=1 to 1000 do begin
****Re:=Re+Reh;
****C:=e1*e0*Re*(((pi*Re)/L)+ln((16*pi*Re)/L)-1);
****if (Re<L) then begin
****edt1.Text:=FormatFloat('0.0000',Re);
****R:=Re;
****end;
**end;
 
**Re:=0;
**Reh:=0.005;
**for i:=1 to 50 do begin
****Re:=Re+Reh;
****C:=e1*e0*Re*(((pi*Re)/L)+ln((16*pi*Re)/L)-1);
****Series1.AddXY(Re,C);
****end;
end;

const e0=8.85e-3;{в пФ/мм}
      e1=1;
      l=1;
      c=100;
      t=0.001;
function F(r:extended):extended;
begin
F:=e1*e0*r*(pi*r/l+ln(16*pi*r/l)-1);
end;
var x,r:extended;
begin
r:=l+t;
while abs(F(r)-c)>t do r:=r+t;
write('R=',r:0:3,'  C=',F(r):0:8);//C для проверки, должно быть около 100
readln
end.