Optimizing a Function of One Variable

#include <iostream>
#include <iomanip>
#include "Brent.h"

class foo
{
public:
    double operator()(double x) {return -x*exp(-x);}
};


int main()
{
	foo f;
	double minloc; // location of the minimum
	double minval; // value at the minimum
	std::cout << std::setprecision(7); // display seven decimal places in output

	// First minimize over [-1, 2]
	minval = Minimize<foo>(f, -1, 2, 1e-6, minloc);
	std::cout << "Computed minimum location: " << minloc << "\n"
		      << "Exact value:               " << 1.0    << "\n"
			  << "Computed minimum value:    " << minval << "\n"
			  << "Exact minimum value:       " << f(1.0) << "\n\n";

	// Now minimize over [2, 3]
	minval = Minimize<foo>(f, 2, 3, 1e-6, minloc);
	std::cout << "Computed minimum location: " << minloc << "\n"
		      << "Exact value:               " << 2.0    << "\n"
			  << "Computed minimum value:    " << minval << "\n"
			  << "Exact minimum value:       " << f(2.0) << "\n\n";

	return 0;
}