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

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

float fp(float x)
{
   return 3.0 * x * x + 1.0;
}

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

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

   while(fabs(f(x0)) > EPS && i < n)	
   {
      x1 = x0 - f(x0) / fp(x0);
      printf("%d %f\n", i+1, x1);
      x0 = x1;
      i++;
   }

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