Raíces de ecuaciones: Método de Steffensen, código C++ | Roots of equations: Steffensen's method, code C++

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
float fun (float x);
int main() {
float error=0.0005,x0,x1;
 printf ("Este codigo fue hecho por Jair Beltran USCO Garzon");
printf ("Digite el primer valor ");
scanf ("%f",&x0);
while (fabs (fun(x0))>error)
{x1=x0-((fun(x0)*fun(x0))/(fun(x0+fun(x0))-fun(x0)));
printf ("%f \n",x1);
x0=x1;
}
x0=x1;
printf ("la raiz es:%f \n",x1);
return 0;
}
float fun(float x)
{float f;
f=x*x-4;
return f;
}