matlab如何实现二分法方程求根

% -------------- inputs -------------------f = @(x) 3*x^2-x-3;a = 0;b = 2;% tolerance / max iterTOL = 1e-4; NI = 50;% -------------------------------------------------------% STEP 1: initializationi = 1;fa = f(a);converge = false; % convergence flag% STEP 2: iterationwhile i<=NI% STEP 3: compute p at the i's stepp = a+(b-a)/2;fp = f(p);% STEP 4: check if meets the stopping criteriaif (abs(fp)<eps || (b-a)/2 < TOL) % eps is Matlab-machine zeroconverge = true; % bisection method converged!break; % exit out of while loopelse% STEP 5i = i+1;% STEP 6if fa*fp > 0a = p; fa = fp;elseb = p;endendend