| Listing 2
A modular minimizer
|
// Global minimizer; Finds lowest
minimum in entire range//
double minimize(double (*f)(double),
double x0, double x1, long n){
static double xmin = x0;
static double ymin = f(x0);
double y, dx;
dx = (x1 - x0)/n;
for(long i=0; i<n; i++){
y = f(x0);
if(y
<
ymin){
ymin = y;
xmin = x0;
}
x0 += dx;
}
return xmin;
}
|
|
Back
|
|
|
|
|
|
