【C/C++】实现牛顿迭代

 #include<bits/stdc++.h>
 using namespace std;

 /*(x+2)^2 +1
     -2(x-1)^2+7

 */
 double f(double x){
     return -x- pow(,x);
 }
 double df(double x){
     return - - pow(,x)*log();
 }

 //二分法
 void binary_search(double l,double r,double e){

     double mid = (l+r)/;
     int manN = ;
     int k = ;
     int flag  = ;
     while(){
         mid = (l+r)/;
         if(fabs(f(mid)) < e){
             flag = ;
             break;
         }
         if(k>=manN){
             flag = ;
             break;
         }
         if(f(l)*f(mid)<){
             r = mid;
         }else if(f(mid)*f(r)<){
             l = mid;
         }
         k++;
     }
     if(flag){
         printf("二分次数是：%d,解为：%.8lf\n",k,mid);
     }
 }

 void NewTon(double x0,double e,int N){
     //牛顿法求根，起始点x0开始
     //N是最大迭代次数
     int flag = ;  //1是奇异 2是找到了 3是失败
     double x1,result;
     int k = ;
     while(){
         if(df(x0) == ){
             flag = ;
             break;
         }
         x1 = x0 - f(x0)/df(x0);
         if(fabs(x1-x0) < e){
             result = x1;
             flag = ;
             break;
         }
         if(k == N){
             flag = ;
             break;
         }
         k++;
         //printf("%lf\n",x0);
         x0 = x1; 

     }
     if(flag == ){
         printf("找到了，并且迭代次数是：%d，根为：%.8lf\n",k,result);
     }else if(flag == ){
         printf("函数不可导");
     }else if(flag == ){
         printf("超出局部收敛范围，迭代失败！\n");
     }else{
         printf("return!");
     }
     return;
 }
 int main(){

 //    double x0;
 //    double e;
 //    int N = 100;
 //    //while(1){
 //        printf("输入起始点、精度、最大迭代次数：");
 //        cin>>x0>>e>>N;
 //        NewTon(x0,e,N);
 //
 //
 //    printf("-------------\n");
 //    printf("请输入二分法左右区间、精度：");
 //    double l,r,e1;
 //    cin>>l>>r>>e1;
 //    binary_search(l,r,e1);
     double x = 1.38616943;
     printf("%.10lf",f(x));

 }