Решение нелинейных уравнений.

uses
    CRT;
var
    E,x0,x1,y0,y1,x:real;
    imax,i,j,k,m:integer;
 
function F(x:real):real;
begin
    F := 1-exp(j*ln(abs(x)))-exp(k*ln(abs(pi*exp(m*ln(abs(x)))/4)));
end;
 
begin
    ClrScr;
        writeln('Enter j,k,m'); readln(j,k,m);
    E:=1.0E-10;
    x0:= 1; x1:= x0+(10-100)*E;
    i:=1; imax:=20;
    y0:=F(x0);
    repeat
        y1:=F(x1);
        x:=x1-(x1-x0)/(y1-y0)*y1;
                writeln('Xn+1 = ',x);
        if abs(x-x1)<E then break;
        x0:=x1; y0:=y1;
        x1:=x;  inc(i);
    until (i=imax);
 
    writeln('X = ',x+E:0:10,'; Iterates: ',i);
    readkey;
end.

program EvgeniiSp;
uses Crt;
const e = 0.0001;
 
//-------------------------
 function F(x: real): real;
  Begin
   F:=3*x-cos(x)-1;
  End;
//-------------------------
 function SGN(Y: real): integer;
  Begin
   if Y < 0
    then SGN:=-1
    else SGN:=1;
  End;
//-------------------------
 procedure Secant_method;
  var
   i: byte;
   x,y,a,b,ya: real;
  begin
   i:=0;
   b:=100;
   a:=-100;
   repeat
    Inc(i);                    //íà÷àëî èòåðàöèè
    x:=(A+B)/2;                // x - ñåðåäèíà îòðåçêà [A,B]
    Y:=F(x); YA:=F(A);         // çíà÷åíèÿ ôóíêöèè â ñåðåäèíå è íà êîíöå îòðåçêà
    if (SGN(Y)*SGN(YA)) > 0    // åñëè çíàêè ôóíêöèè â òî÷êàõ "A" è "x" ñîâïàäàþò
     then  A:=x  else  B:=x;   // òî,  ïåðåíîñ ãðàíèöû "A",  èíà÷å  - "Â"
   until (Abs(B-A) < e);
   writeln('X = ',x+E:0:10,'; Êîë. èòåðàöèé = ',i);
  end;
 
begin
 ClrScr;
 Secant_method;
end.