[수치계산] 오차&비선형 방정식의 해 - 고정점 반복법

고정점 반복법 (fixed - point iteration)

● 수렴? 발산?

- 고정점 반복법은 초기값의 선택에 발산할 수도 있음을 주의!

#include <stdio.h>
#include <math.h>

#define THRESHOLD 0.000001

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

void main() {
    double x = 0;
    double x2 = 0;
    int i = 1;

    printf("i \t xi \t\t g(x3) \n");
    printf("%d \t %f \t %f \n", i, x, g(x));

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

    return 0;
}