
        writeln(fTABL_Sol, '');                                             // ------------------------- fTABL_Sol
        writeln(fTABL_Sol, '  ISTEP = ' + IntToStr(ISTEP));
        writeln(fTABL_Sol, '  IIND  = ' + IntToStr(IIND));
        writeln(fTABL_Sol, '  JIND  = ' + IntToStr(ISTEP));
        writeln(fTABL_Sol, '');
        writeln(fTABL_Sol, '  QMOD  = ' + FloatToStrF(QMOD, ffexponent, 6, 2));
        writeln(fTABL_Sol, '');
        //writeln(fTABL_Sol, '    D(I)       V(I,' + IntToStr(NIC[1]) + ')       V(I,' + IntToStr(NIC[2]) + ')       V(I,' + IntToStr(NIC[3]) + ')');
        Str1 := '    D(I)       V(I,' + IntToStr(NIC[1]) + ')';
        if NIC[2] <> 0 then Str1 := Str1 + '       V(I,' + IntToStr(NIC[2]) + ')';
        if NIC[3] <> 0 then Str1 := Str1 + '       V(I,' + IntToStr(NIC[3]) + ')';
        writeln(fTABL_Sol, Str1);
        Str2 := ' +---------+------------+';
        if NIC[2] <> 0 then Str2 := Str2 + '------------+';
        if NIC[3] <> 0 then Str2 := Str2 + '------------+';
        writeln(fTABL_Sol, Str2);
        IH := 1;
        for I := 1 to N5 do              //STEP ITABL
        begin
          if I < IH then continue;
          {writeln(fTABL_Sol, '  ' + FloatToStrF(D[IH],         ffexponent, 3, 2) +
                             '  ' + FloatToStrF(V[IH, NIC[1]], ffexponent, 6, 2) +
                             '  ' + FloatToStrF(V[IH, NIC[2]], ffexponent, 6, 2) +
                             '  ' + FloatToStrF(V[IH, NIC[3]], ffexponent, 6, 2));}
          Str1 := '  ' + FloatToStrF(D[IH], ffexponent, 3, 2);
          if V[IH, NIC[1]] >= 0 then Str1 := Str1 + '  ' + FloatToStrF(V[IH, NIC[1]], ffexponent, 6, 2)
          else                       Str1 := Str1 + '  ' + FloatToStrF(V[IH, NIC[1]], ffexponent, 5, 2);
          if NIC[2] <> 0 then
            if V[IH, NIC[2]] >= 0 then Str1 := Str1 + '  ' + FloatToStrF(V[IH, NIC[2]], ffexponent, 6, 2)
            else                       Str1 := Str1 + '  ' + FloatToStrF(V[IH, NIC[2]], ffexponent, 5, 2);
          if NIC[3] <> 0 then
            if V[IH, NIC[3]] >= 0 then Str1 := Str1 + '  ' + FloatToStrF(V[IH, NIC[3]], ffexponent, 6, 2)
            else                       Str1 := Str1 + '  ' + FloatToStrF(V[IH, NIC[3]], ffexponent, 5, 2);
          writeln(fTABL_Sol, Str1);
          IH := IH + ITABL;
        end;
        if ITABL * INT(N5 / ITABL) = N5 then
        begin
          I := N5;
          {writeln(fTABL_Sol, '  ' + FloatToStrF(D[I],         ffexponent, 3, 2) +
                             '  ' + FloatToStrF(V[I, NIC[1]], ffexponent, 6, 2) +
                             '  ' + FloatToStrF(V[I, NIC[2]], ffexponent, 6, 2) +
                             '  ' + FloatToStrF(V[I, NIC[3]], ffexponent, 6, 2));}
          Str1 := '  ' + FloatToStrF(D[I], ffexponent, 3, 2);
          if V[I, NIC[1]] >= 0 then Str1 := Str1 + '  ' + FloatToStrF(V[I, NIC[1]], ffexponent, 6, 2)
          else                      Str1 := Str1 + '  ' + FloatToStrF(V[I, NIC[1]], ffexponent, 5, 2);
          if NIC[2] <> 0 then
            if V[I, NIC[2]] >= 0 then Str1 := Str1 + '  ' + FloatToStrF(V[I, NIC[2]], ffexponent, 6, 2)
            else                       Str1 := Str1 + '  ' + FloatToStrF(V[I, NIC[2]], ffexponent, 5, 2);
          if NIC[3] <> 0 then
            if V[IH, NIC[3]] >= 0 then Str1 := Str1 + '  ' + FloatToStrF(V[I, NIC[3]], ffexponent, 6, 2)
            else                       Str1 := Str1 + '  ' + FloatToStrF(V[I, NIC[3]], ffexponent, 5, 2);
          writeln(fTABL_Sol, Str1);
        end;
        writeln(fTABL_Sol, '');    writeln(fTABL_Sol, '  NDEL  = ' + IntToStr(N));
      end;

    3: // - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 3
      begin    //ShowMessage('MARK_ = 3');
        Form4.Left    := Form1.Left + 417;
        Form4.Top     := Form1.Top  + 577;
        Form4.Button1.Caption := ' Solution # ' + IntToStr(ISOLV);
        Form4.ShowModal;

        writeln(fTABL_Qsl, '');                                             // ------------------------- fTABL_Qsl
        writeln(fTABL_Qsl, '  ISOLV = ' + IntToStr(ISOLV));
        writeln(fTABL_Qsl, '');
        writeln(fTABL_Qsl, '  QMOD  = ' + FloatToStrF(QSOLV, ffexponent, 6, 2));
        writeln(fTABL_Qsl, '');
        //writeln(fTABL_Qsl, '    D(I)       Y(I,' + IntToStr(NIC[1]) + ')       Y(I,' + IntToStr(NIC[2]) + ')       Y(I,' + IntToStr(NIC[3]) + ')');
        {Str1 := C1;
        for I := 2 to NDIM do
          Str1 := Str1 + dC1 + IntToStr(I) + dE1;
        writeln(fTABL_Qsl, Str1);
        Str2 := D1;
        for I := 2 to NDIM do
          Str2 := Str2 + dD1;
        writeln(fTABL_Qsl, Str2);}
        Str1 := '    D(I)       Y(I,1)';
        for I := 2 to NDIM do
          Str1 := Str1 + '       ' + 'Y(I,' + IntToStr(I) + ')';
        writeln(fTABL_Qsl, Str1);
        Str2 := ' +---------+------------+';
        for I := 2 to NDIM do
          Str2 := Str2 + '------------+';
        writeln(fTABL_Qsl, Str2);
        for I := 1 to N5 do
        begin
          Str1 := '  ' + FloatToStrF(D[I],    ffexponent, 3, 2);
          for J := 1 to NDIM do
            if Y[I, J] >= 0 then Str1 := Str1 + '  ' + FloatToStrF(Y[I, J], ffexponent, 6, 2)
            else                 Str1 := Str1 + '  ' + FloatToStrF(Y[I, J], ffexponent, 5, 2);
          writeln(fTABL_Qsl, Str1);
        end;
        writeln(fTABL_Qsl, '');
        writeln(fTABL_Qsl, '  NDEL  = ' + IntToStr(N));
      end;
