#!/usr/local/bin/tcc -run #include #include #define eps 1.0e-15 double f(double x) { return x*x*x-2.0*x*x+2.0; } int main() { int n=0; double x0=-2.0,x1=-1.8,x2,f0,f1,diff; do{ f1=f(x1);f0=f(x0); x2=x1-f1*(x1-x0)/(f1-f0); diff=fabs(f(x2)); x0=x1; x1=x2; printf("%25.15e %10.3e \n",x1,diff); }while(diff>=eps); printf("%25.15e %10.3e \n",x1,diff); }