fork(1) download
copy 
  1. #include <iostream>
  2. #include <math.h>
  3. using namespace std;
  4.  
  5. int main () 
  6. {
  7.     double a,xz,eps;
  8.     int p;
  9.     cin>>a>>p>>xz>>eps;
  10.     double x=pow(a,1.0/p);
  11.     double xn=xz,x_prev=-10000000;
  12.     int i=0;
  13.     while (fabs(xn-x_prev)>=eps)
  14.     {  
  15.         x_prev=xn;
  16.         xn=(x_prev/(p*p))*((p*p-1)+0.5*(p+1)*a/(pow(x_prev,p))-0.5*(p-1)*pow(x_prev,p)/(a));
  17.         i++;
  18.  
  19.     }
  20.     cout<<i<<' '<<xn<<' '<<x<<endl; 
  21.     return 0;
  22. }
Success #stdin #stdout 0s 3476KB
comments (?)
stdin
copy 
17 2 2 0.1
stdout
copy 
3 4.12311 4.12311
https://ideone.com/LJiEMx
language:
C++14 (gcc 8.3)
created:
9 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!