void bibun(double n);
double f(double x, double n);
 
int main()
{
	bibun(2.0);
	bibun(6.0);
	return 0;
}
 
void bibun(double n)
{
	double	x = 1.0;
	double	dx = 0.01;
	double	fd, cd, bd;
 
	fd = (f(x + dx, n) - f(x, n)) / dx;
	cd = (f(x + dx, n) - f(x - dx, n)) / (2.0 * dx);
	bd = (f(x, n) - f(x - dx, n)) / dx;
	printf("n=%f\n", n);
	printf("前進:%f\n", fd);
	printf("中央:%f\n", cd);
	printf("後退:%f\n", bd);
}
 
double f(double x, double n)
{
	double	y;
 
	y = 2.0 * pow(x, n) + 3.0 * x + 5.0;
	return y;
}
 

dm9pZCBiaWJ1bihkb3VibGUgbik7CmRvdWJsZSBmKGRvdWJsZSB4LCBkb3VibGUgbik7CgppbnQgbWFpbigpCnsKCWJpYnVuKDIuMCk7CgliaWJ1big2LjApOwoJcmV0dXJuIDA7Cn0KCnZvaWQgYmlidW4oZG91YmxlIG4pCnsKCWRvdWJsZQl4ID0gMS4wOwoJZG91YmxlCWR4ID0gMC4wMTsKCWRvdWJsZQlmZCwgY2QsIGJkOwoKCWZkID0gKGYoeCArIGR4LCBuKSAtIGYoeCwgbikpIC8gZHg7CgljZCA9IChmKHggKyBkeCwgbikgLSBmKHggLSBkeCwgbikpIC8gKDIuMCAqIGR4KTsKCWJkID0gKGYoeCwgbikgLSBmKHggLSBkeCwgbikpIC8gZHg7CglwcmludGYoIm49JWZcbiIsIG4pOwoJcHJpbnRmKCLliY3pgLI6JWZcbiIsIGZkKTsKCXByaW50Zigi5Lit5aSuOiVmXG4iLCBjZCk7CglwcmludGYoIuW+jOmAgDolZlxuIiwgYmQpOwp9Cgpkb3VibGUgZihkb3VibGUgeCwgZG91YmxlIG4pCnsKCWRvdWJsZQl5OwoKCXkgPSAyLjAgKiBwb3coeCwgbikgKyAzLjAgKiB4ICsgNS4wOwoJcmV0dXJuIHk7Cn0K