[수치계산] 오차&비선형 방정식의 해 - 할선법

할선법 (secant method)

● 가위치법의 개선 방법

● 근을 포함하는 구간을 구할 필요 없이 임의의 두 점에서 시작

#include <stdio.h>
#include <math.h>
#define THRESHOLD 0.000001

double f(double x) {
    double ret = 0;
    ret = x - exp(-x);
    return ret;
}

int main() {
    double x1 = 0;
    double x2 = 1;
    double x3 = 0;
    int i = 1;

    printf("i \t x1 \t\t x2 \t\t x3 \t\t f(x3) \n");
    while (1) {
        x3 = x2 - ((f(x2) * (x2 - x1)) / (f(x2) - f(x1)));
        printf("%d \t %f \t\t %f \t\t %f \t\t %f \n", i, x1, x2, x3, f(x3));
        if (THRESHOLD > fabs(f(x3)))    break;
        x1 = x2;
        x2 = x3;
        i++;
    }
}

결과