#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) 回帖