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;
}
double f(double x, double n)
{
double y;
y
= 2.0 * pow(x
, n
) + 3.0 * x
+ 5.0; return y;
}
dm9pZCBiaWJ1bihkb3VibGUgbik7CmRvdWJsZSBmKGRvdWJsZSB4LCBkb3VibGUgbik7CgppbnQgbWFpbigpCnsKCWJpYnVuKDIuMCk7CgliaWJ1big2LjApOwoJcmV0dXJuIDA7Cn0KCnZvaWQgYmlidW4oZG91YmxlIG4pCnsKCWRvdWJsZQl4ID0gMS4wOwoJZG91YmxlCWR4ID0gMC4wMTsKCWRvdWJsZQlmZCwgY2QsIGJkOwoKCWZkID0gKGYoeCArIGR4LCBuKSAtIGYoeCwgbikpIC8gZHg7CgljZCA9IChmKHggKyBkeCwgbikgLSBmKHggLSBkeCwgbikpIC8gKDIuMCAqIGR4KTsKCWJkID0gKGYoeCwgbikgLSBmKHggLSBkeCwgbikpIC8gZHg7CglwcmludGYoIm49JWZcbiIsIG4pOwoJcHJpbnRmKCLliY3pgLI6JWZcbiIsIGZkKTsKCXByaW50Zigi5Lit5aSuOiVmXG4iLCBjZCk7CglwcmludGYoIuW+jOmAgDolZlxuIiwgYmQpOwp9Cgpkb3VibGUgZihkb3VibGUgeCwgZG91YmxlIG4pCnsKCWRvdWJsZQl5OwoKCXkgPSAyLjAgKiBwb3coeCwgbikgKyAzLjAgKiB4ICsgNS4wOwoJcmV0dXJuIHk7Cn0K