C语言实践之函数

牛顿法求方程近似解

牛顿法是一种求解方程近似根的方法。小蒜的老师给他出了一道题目，让他用牛顿法求下面这个方程的近似解。

2x^3 - 4x^2 + 3x - 6 = 0

根据牛顿法的定义，小蒜对 f(x) = 2x^{3} - 4x^{2} + 3x - 6进行了求导，得到了f'(x) = 6x^{2} - 8x + 3

如给出的代码所示，小蒜已经将 f(x) 和f′(x) 写成了两个 C 语言的函数。并且完成一个newton函数，实现了牛顿法求解的目的，并在main函数中进行了调用。newton函数中，x 的初值被设为了 1.5，按照牛顿法进行近似解迭代会到 10^{-6} 的误差范围时停止。你可以运行一下程序，看看小蒜实现的程序的计算结果（带入方程看看他计算得是否正确呢）。

#include <stdio.h>
#include <math.h>
#define EPSILON 1e-6

double f(double x) {
    return 2 * pow(x, 3) - 4 * pow(x, 2) + 3 * x - 6;
}

double f_prime(double x) {
    
    return 6 * pow(x, 2) - 8 * x + 3;
}
double h(double x) {
    return pow(x, 3) - 4 * pow(x, 2) + 3 * x - 6;
}

double h_prime(double x) {
    return 3 * pow(x, 2) - 8 * x + 3;
}

double newton(double (*fp)(double), double (*fp_prime)(double)) {
    double x = 1.5;
   while (fabs(fp(x)) > EPSILON) {
    x = x - fp(x) / fp_prime(x);
}
    return x;
}

int main() {
    printf("%g\n", newton(f,f_prime));
    printf("%g\n", newton(h, h_prime));
    return 0;
}

二分法求方程近似解

二分法是一种求解方程近似根的方法。对于一个函数 f(x)，使用二分法求 f(x) 近似解的时候，我们先设定一个迭代区间（在这个题目上，我们之后给出了的两个初值决定的区间 [−20,20]），区间两端自变量 x 的值对应的 f(x) 值是异号的，之后我们会计算出两端 x 的中点位置 x′ 所对应的 f(x') ，然后更新我们的迭代区间，确保对应的迭代区间的两端 x 的值对应的f(x) 值还会是异号的。

重复这个过程直到我们某一次中点值 x' 对应的 f(x′)<ϵ （题目中可以直接用EPSILON）就可以将这个 x'作为近似解返回给 main 函数了。

例如：