Соотношение Ньютона - C++

#include <stdio.h>
#include <conio.h>
#include <math.h>
 
int main()
{
        float xn, xl, x, eps, a, p;
        int count=0;
        printf("vvedit znachennia a\n");
        scanf("%f",&a);
        printf("vvedit znachennya p\n");
        scanf("%f",&p);
        printf("vvedit pogrewnost eps\n");
        scanf("%f",&eps);
        xn=a;
        x=pow(a,1./p);
 
        do {
                xl=xn;
                xn=(1./p)*((p-1)*xl+(a/pow(xl,p-1)));                
                count++;
        }while (fabs(xn-xl)>eps);
 
        printf("x=%f\n",x);
        printf("xn=%f\n",xn);
        printf("%d\n",count);
        getch();
        return 0;
}