【C++较复杂版解法】

参考代码:

#include <iostream>
#include <math.h>
#include <iomanip>

using namespace std;
int main()
{
 double a,b,x1,x2;
 cin>>a;
 x1=a/2.0;
 x2=(x1+a/x1)/2.0;
 while(x1-x2>0.00001)
 {
    x1=(x2+a/x2)/2;
    x2=(x1+a/x1)/2;
   if(x1-x2<=0.00001)
   {
   break;
   }
 }
 if(x1-x2<=0.00001)
 {
  cout<<fixed<<setprecision(3)<<x1<<endl;
 }
 return 0;
}