#include <cmath>
#include <cstdlib>
#include <iostream>
using namespace std;

int main() {
    double phi = (sqrt(5.0) - 1) / 2;
    double phi_nMinus1 = 1, phi_n = phi;

    cout << "n\tRecursion\tMultiplication\tPercent Error" << endl;
    for (int n = 1; n < 45; n++) {
        double phi_mult = pow(phi, n);
        double percentError = abs((phi_n - phi_mult ) / phi_mult * 100);
        cout << n << '\t' << phi_n << '\t' << phi_mult << '\t'
             << percentError << '\n';
        double phi_nPlus1 = phi_nMinus1 - phi_n;
        phi_nMinus1 = phi_nMinus1 = phi_n;
        phi_n = phi_nPlus1;
    }
}
