public class RFuncLib{

public static double bracketRoot(double a, double b, double maxErr,Function f){

double fl = f.evaluate(a);

double fu = f.evaluate(b);

double upper = b;

double lower = a;

double x = 0.0;

double fx;

//System.out.println("delta = " + (upper - lower));

while((upper-lower) > maxErr){

x = (upper+lower)/2;

//System.out.println("x = " + (x));

fx = f.evaluate(x);

if(fx == 0) return(x);

else if((fx>0) && (fu>0)){ upper = x; fu = fx;}

else if((fx<0) && (fu<0)){ upper = x; fu = fx;}

else { lower = x; fl = fx;}

}

return(x);

}

public static double defIntegral(double a, double b, int nsteps,Function f){

int i;

double h = (b-a)/nsteps;

double x = a;

double sum = (f.evaluate(a)+f.evaluate(b))/2;

for(i=1;i<nsteps;i++){

x += h;

sum += f.evaluate(x);

}

return(sum * h);

}

}