二分法改进源代码
//二分法求非线性方程组近似解
#include<stdio.h>
#include<math.h>
double F(double x)
{
double y=x*x*x+2*x*x+10*x-20;//可自定义
return y;
}
int main()
{
double xleft=0;
double xright=0;
double xmid=0;
double yleft=F(xleft);
double yright=F(xright);
double ymid=F(xmid);
printf("输入左右区间\n");
scanf("%lf%lf",&xleft,&xright);
xmid=(xleft+xright)/2;
double wc=0;//误差限
printf("输入想要的精度:\n");
scanf("%lf",&wc);
int maxdd=0;//迭代最大次数
int count=0; //迭代次数
printf("输入最大迭代次数:\n");
scanf("%d",&maxdd);
while(1)
{
count++;
if(F(xmid)*F(xleft)<0)//迭代过程
{
xright=xmid;
xmid=(xleft+xright)/2;
}
else if(F(xmid)*F(xright)<0)
{
xleft=xmid;
xmid=(xleft+xright)/2;
}
else if(F(xmid)==0)
{
printf("精确解为:%lf\n",xmid);
break;
}//迭代过程
if(F(xmid)<=wc&&F(xmid)>=(-wc))
{
printf("迭代次数为:%d\n",count);
printf("二分法迭代结果为:%lf\n",xmid);
break;
}
else if(count>maxdd)
{
printf("超出最大迭代次数\n");
break;
}
}
return 0;
}
结果展示
牛顿法改进源代码
//牛顿法求非线性方程组近似解
#include<stdio.h>
#include<math.h>
double F(double x)//原函数
{
double y=x*x*x+2*x*x+10*x-20;//自定义
return y;
}
double F_(double x)//原函数的导函数
{
double y_=3*x*x+4*x+10;//根据所要迭代的函数自行求解
return y_;
}
int main()
{
double x;//x初值
printf("输入x初值:\n");
scanf("%lf",&x);
double wc;//误差限
printf("输入精度:\n");
scanf("%lf",&wc);
int maxdd=0;//最大迭代次数
int count=0;
printf("输入最大迭代次数\n");
scanf("%d",&maxdd);
while(1)
{
count++;
x=x-F(x)/(F_(x)); //牛顿法迭代公式
if(F(x)<=wc&&F(x)>=(-wc))
{
printf("迭代次数为:%d\n",count);
printf("迭代结果为:%lf\n",x);
break;
}
else if(count>maxdd)
{
printf("超出最大迭代次数\n");
break;
}
}
return 0;
}
结果展示
割线法改进源代码
//割线法求非线性方程组近似解
#include<stdio.h>
#include<math.h>
double F(double x)
{
double y=y=x*x*x+2*x*x+10*x-20;//自定义
return y;
}
int main()
{
double x1;//左区间初始值
double x2;//右区间初始值
double x3; //迭代后的值
double s;
int i;
int n=0;
printf("输入左右两个初始值x:\n");
scanf("%lf%lf",&x1,&x2);
double wc;//误差限
printf("输入精度:\n");
scanf("%lf",&wc);
int maxdd=0;//最大迭代次数
int count=0;
printf("输入最大迭代次数\n");
scanf("%d",&maxdd);
while(1)
{
count++;
if(F(x1)<F(x2))//迭代过程
{
s=(F(x2)-F(x1))/(x2-x1);
x3=x1-F(x1)/s;
x2=x1;
x1=x3;
}
else if(F(x1)>F(x2))
{
s=(F(x2)-F(x1))/(x2-x1);
x3=x1-F(x1)/s;
x1=x2;
x2=x3;
}
else if(F(x1)==F(x2))
{
printf("该函数在该区间上无法用割线法迭代\n");
break;
}//迭代过程
if(F(x3)<=wc&&F(x3)>=(-wc))
{
printf("迭代次数为:%d\n",count);
printf("迭代结果为:%lf\n",x3);
break;
}
else if(count>maxdd)
{
printf("超出最大迭代次数\n");
break;
}
}
return 0;
}
结果展示
牛顿法牛批!(又能水一篇文章,爽死)
————如有错误,欢迎指正(#^.^#)