C语言实现：用弦截法求方程的根

最新推荐文章于 2024-03-21 13:25:57 发布

爱吃白菜的金小妞

最新推荐文章于 2024-03-21 13:25:57 发布

阅读量1.9k

点赞数 7

分类专栏： C语言编程文章标签： c语言开发语言

本文链接：https://blog.csdn.net/weixin_49191101/article/details/129750403

版权

C语言编程专栏收录该内容

7 篇文章 0 订阅

订阅专栏

文章介绍了如何运用弦截法求解立方方程x^3-5x^2+16x-80=0。首先，选取两个函数值异号的点x1和x2，然后通过连接两点作弦，找到弦与x轴的交点x，不断迭代更新x1和x2，直到函数值的绝对值小于给定的阈值0.0001，从而得到方程的近似根。代码示例展示了具体的实现过程。

摘要由CSDN通过智能技术生成

问题描述

用弦截法求方程 x 3 -5x 2 +16x-80=0 的根
在这里插入图片描述

思路分析

设f(x)=x3-5x2+16x-80，用弦截法求方程根的方法如下：
①取相差不大的两个不同点x1、x2，如果f(x1)与f(x2)符号相反，则(x1,x2)区间内必有一个根。否则改变x1、x2的值使得f(x1)与f(x2)异号为止。

②连接f(x1)与f(x2)两点，此线(即弦)交x轴于x点(计算见图5-5)。再求出f(x)。

③若f(x)与f(x1)同号，则将x作为新的x1。若f(x)与f(x2)同号，则将x作为新的x2。

④重复第②与第③步，直到|f(x)|<e，其中e是一个很小的正数，如0.0001，此时可以认为： f(x)≈0。

代码

#include<stdio.h>
#include<math.h>
float f(float x);
float xpoint(float x1,float x2);
float root(float x1,float x2);
void main(){
	float x1,x2,f1,f2,x;//输入x1和x2并保证其函数值反号
	do{
		printf("Input x1,x2:\n");
		scanf("%f,%f",&x1,&x2); 
		f1=f(x1);     //计算x1的函数值并赋给f1 
		f2=f(x2);
	}while(f1*f2>=0);
	x=root(x1,x2);//调用求根函数
	printf("A root of equation is %8.4f",x); 
}

//定义f函数，求f(x)的函数值 
float f(float x){
	float y;
	y=((x-5.0)*x+16.0)*x-80;
	return y;
}

// 定义xpoint函数，求弦与x轴的交点
float xpoint(float x1,float x2){
	float x;
	x=(x1*f(x2)-x2*f(x1))/(f(x2)-f(x1));
	return x;
} 

//定义root函数，求近似根
float root(float x1,float x2){
	float x,y,y1;
	y1=f(x1);
	do{
		x=xpoint(x1,x2);
		y=f(x);
		if(y*y1>0) //新的f(x)与f(x1)同号时，用x代替x1，否则替换x2 
		{
			y1=y;
			x1=x; 
		}
		else{
			x2=x;
		}
	}while(fabs(y)>=0.0001);  //当f(x)约等于0时退出循环 
	return x; 
}