#include #include double f(x,q) double x,q; { return(x-(1/(x*x))+q*(x-1-(1/((x-1)*(x-1))))); } double g(x,q) double x,q; { return(1+2/(x*x*x)+2*q*(1/((x-1)*(x-1)*(x-1)))); } void main() { int k; double x,q; printf("x=\n"); scanf("%lf",&x); printf("q=\n"); scanf("%lf",&q); for (k=0 ;k<10; ++k) { x=x-f(x,q)/g(x,q); printf("%17.12lf %15.8e\n",x,f(x)); } }