#include<stdio.h>
#include<math.h>
#define EPS 1e-6

float f(float x)
{
   return x * (x * x  + 1) - 1;
}

int main()
{
   float x0, x1;
   int i, n=100;

   printf("Punkt startowy : ");
   scanf("%f", &x0);

   for (i = 0; i < n; i++)
   {
      x1 = x0 - f(x0) / (3*x0*x0 + 1);
      if (fabs(x1 - x0) < EPS) break;
      x0 = x1;
   }

   if (i == n) printf("Nie znaleziono rozwiazania\n");   
   else 
   {
      printf("x = %f\n", x1);
      printf("iterations %d\n", i+1);
   }
}