#include <bits/stdc++.h>
using namespace std;
double f(double x){return x*x*x - 4;}
int main(){
double x1=1,x2=3,x0,f1,f2,f0,E=1e-8;
f1=f(x1); f2=f(x2);
cout<<"Iter\t x1\t\t x2\t\t x0\t\t f1\t\t f2\t\t f0\n";
for(int i=1;i<=3;i++){
x0=x1-(f1*(x2-x1))/(f2-f1); f0=f(x0);
cout<<i<<"\t"<<x1<<"\t"<<x2<<"\t"<<x0<<"\t"<<f1<<"\t"<<f2<<"\t"<<f0<<"\n";
if(f1*f0<0) x2=x0; else x1=x0;
f1=f(x1); f2=f(x2);
}
}
I2luY2x1ZGUgPGJpdHMvc3RkYysrLmg+CnVzaW5nIG5hbWVzcGFjZSBzdGQ7CmRvdWJsZSBmKGRvdWJsZSB4KXtyZXR1cm4geCp4KnggLSA0O30KaW50IG1haW4oKXsKICAgIGRvdWJsZSB4MT0xLHgyPTMseDAsZjEsZjIsZjAsRT0xZS04OwogICAgZjE9Zih4MSk7IGYyPWYoeDIpOwogICAgY291dDw8Ikl0ZXJcdCB4MVx0XHQgeDJcdFx0IHgwXHRcdCBmMVx0XHQgZjJcdFx0IGYwXG4iOwogICAgZm9yKGludCBpPTE7aTw9MztpKyspewogICAgICAgIHgwPXgxLShmMSooeDIteDEpKS8oZjItZjEpOyBmMD1mKHgwKTsKICAgICAgICBjb3V0PDxpPDwiXHQiPDx4MTw8Ilx0Ijw8eDI8PCJcdCI8PHgwPDwiXHQiPDxmMTw8Ilx0Ijw8ZjI8PCJcdCI8PGYwPDwiXG4iOwogICAgICAgIGlmKGYxKmYwPDApIHgyPXgwOyBlbHNlIHgxPXgwOwogICAgICAgIGYxPWYoeDEpOyBmMj1mKHgyKTsKICAgIH0KfQ==