C语言，数值分析，二分法，牛顿法，割线法解线性方程解法改进

二分法改进源代码

//二分法求非线性方程组近似解 
#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;
}

 结果展示

 

牛顿法牛批！（又能水一篇文章，爽死）

                                                                                ————如有错误，欢迎指正(#^.^#)