// calculate square
double f(double x) {
	return x*x;
}