3637383940414243444546
Dbrent dbr = new Dbrent(); for (i=0;i<N;i++) { a=i; b=a+span; dbr.bracket(a,b,funcd); min=dbr.minimize(funcd); // Minimum of bessj0 localflag = localflag || (abs(funcd.df(min)) > sbeps); } globalflag = globalflag || localflag; if (localflag) {