用二分法求方程X^5-3X+1=0在（0,1）上的近似解,精确到C=0.001,写出算法,画出流程图,并写出算法语句

取定范围,用整型,先乘1000就行了
var
x,rh,lf:longint;
function f(p:longint);
begin
f:=p^5-3p+1;
end;
begin
lf:=0;
rh:=1000;
while rh>lf+1 do
begin
x：=（lf+rh) div 2;
if f(x/1000)>0 then rh:=x
else lf:=x;
end;
write('x=',lf);
readln;
end.
fp的