clc;
f=@(x)?sqrt(x^2+9)-tan(x);
a=0;
b=pi/2-eps;
tol=1e-4;
while?abs(a-b)>tol
if?feval(f,a)*feval(f,(a+b)/2)>0
a=(a+b)/2;
elseif?feval(f,a)*feval(f,(a+b)/2)<0
b=(a+b)/2;
else
break;
end
x=(a+b)/2;
fval=feval(f,x);
disp('近視根為')
disp(x)
return