Расчет параметра спиральной коаксиальной линии — Pascal(Паскаль)

program prog1;
var N,D,E,M,A,X,H,F:real;
begin
10: writeln('РАСЧЕТ ПАРАМЕТРА X СПИРАЛЬНОЙ КОАКСИАЛЬНОЙ ЛИНИИ');
20: write('ВВЕДИТЕ N, D, E');
    readln(N, D, E);
    M := N*N;
30: write('BBEДИТЕ FW =');
    readln(A);
    X := 0;
    H := 0.1;
40: IF ABS(H) < E THEN goto 70;
50: X := X + H;
    F := A-SQRT(1+M*sqr(PI*X)*(1-sqr(X/D))/2/LOG(D/X));
    IF F/H >= 0 THEN goto 50;
60: H := -H /10;
    GOTO 40;
70: write('X=', X);
    {GOTO 30; а это зачем }
end.