0001 function [P,T] = fin07bisect(a,b,tol,maxT) 0002 while (b-a)/2 > tol, 0003 c=(a+b)/2; 0004 [w,T] = fin07(0,4,0,4,0.1,25,25,c,3.85,0.005,0.75,1.0); 0005 if T > maxT, 0006 b=c; 0007 else 0008 a=c; 0009 end 0010 end 0011 P = (a+b)/2; 0012 [w,T] = fin07(0,4,0,4,0.1,25,25,P,3.85,0.005,0.75,1.0); 0013 end