二分法求函数值的Pascal实现

用二分法求在（a,b)上单调的函数近似值

 

第八行的表达式可更改，第三行的kexi决定的精度，小数值计算可将第五行的extended更为real或double

PROGRAM EQUANTION (input,output);

CONST

kexi=0.0000001;

VAR

a,b,c:extended;

FUNCTION fx(x:extended):extended;

BEGIN

fx:=ln(x)/ln(2);

END;

BEGIN

writeln('Please input a and b’，'    Produced BY Chaos')

read(a,b);

c:=(a+b)/2;

IF fx(c)=0

THEN writeln('x=',c)

ELSE Repeat

IF fx(a)*fx(c)<0

THEN b:=c

ELSE a:=c;

UNTIL abs(a-b)<kexi;

IF fx(c)<>0

THEN writeln('x=',a);

END.