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

#define EPSILON 1e-7   
#define MAX_ITER 100

// 【確定】あなたのプログラムが示してくれた挙動に100%合致する式
double f(double lambda) {
    return pow(lambda, 4) + pow(lambda, 3) - 88.0 * pow(lambda, 2) - 132.0 * lambda + 603.0;
}

// その導関数
double df(double lambda) {
    return 4.0 * pow(lambda, 3) + 3.0 * pow(lambda, 2) - 176.0 * lambda - 132.0;
}

void solve_newton(int attempt, double lambda0) {
    double lambda = lambda0;
    int iter = 0;

    while (iter <= MAX_ITER) {
        double f_val = f(lambda);
        double df_val = df(lambda);

        if (fabs(df_val) < 1e-12) {
            return; 
        }

        double delta = f_val / df_val;
        lambda -= delta;
        iter++;

        if (fabs(delta) < EPSILON) {
            printf("試行 #%d (初期値 %5.1f) -> 収束成功 (%2d回) | 固有値解 λ = %.10f\n", 
                   attempt, lambda0, iter, lambda);
            return;
        }
    }
}

int main() {
    // あなたが教えてくれた「8.6、1.台、-5.台が2つ」をピンポイントで仕留める初期値
    // 特に、2つの -5.something を確実に引き離すための -5.2 と -6.0 です
    double initial_lambda[] = {10.0, 1.5, -5.2, -6.0};

    printf("=== 【1つ目の行列】ニュートン・ラフソン法 探索結果 ===\n");
    for (int i = 0; i < 4; i++) {
        solve_newton(i + 1, initial_lambda[i]);
    }

    return 0;
}