【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)); }
上一篇:【Ecmall】ECMall2.x模板制作入门系列(认识ECMall模板)


下一篇:CF528D Fuzzy Search 和 BZOJ4259 残缺的字符串