#include<iostream>
#include<cmath>
using namespace std;
#define JD 0.001
#define MAXCount 20000
double Next(double pre,double x){
double ans=0;
ans=pre-(pow(pre,7)+0.5*pre-x)/(7*pow(pre,6)+0.5);
return ans;
}
bool Jd(double pre,double x){
double value=(pow(pre,7)+0.5*pre-x);
if(abs(value)<JD){
return true;
}
return false;
}
bool Wc(double Tyn,double pre,double x){
double value0=(pow(Tyn,7)+0.5*Tyn-x);
double value=(pow(pre,7)+0.5*pre-x);
if(abs(value0)<abs(value)){
return true;
}
return false;
}
int main()
{
double x=0;
cin>>x;
double PYn=0.5*x;
double Tyn=PYn;
int count=0;
while (!Jd(Tyn,x)&&count<MAXCount) {
PYn=Tyn;
Tyn=Next(Tyn,x);
double mark=0.5;
if(!Wc(Tyn,PYn,x)){
Tyn=(Tyn-PYn)*mark+PYn;
mark/=2;
}
++count;
}
cout<<Tyn;
}
全部评论
(0) 回帖