matlab用二分法求方程的近似根,保留小数点后四位有效数字.

function [x,fval]=myerfen()
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
end
x=(a+b)/2;
fval=feval(f,x);
disp('近视根为')
disp(x)
return