1. function [a,b]=findbracket(f,x0) dx=1e-3; a=x0; b=x0; while 1 fprintf('Bracket = [%12.8f , %12.8f ]\ ',a,b); a=a-dx; if f(a)*f(b)<0, break; end b=b+dx; if f(a)*f(b)<0, break; end dx=2*dx; end fprintf('Bracket = [%12.8f , %12.8f ]\ ',a,b); 2. function [r,h] =schroderbisection(a,b,f,fp,fpp,t) if f(a)*f(b)>0, error('Equal signs on brackets'); end p=a; fprintf(' n method a b p f(p) \ '); fprintf('------------------------------------------------------------------\ '); for n=1:100 h=h-(f*fp)/(fb*fp-f*fpp); if ~(a <= p & b >= p) p=(a+b)/2; method='Bisect'; else method=' schroder'; end if f(p)*f(b)<0 a else b=p; end fprintf('%2d %s %12.8f %12.8f\ ',n,method,a,b,p,f(p)); if abs(f(p))