迭代求方程的根c++

迭代求方程的根c++

题目

求x³-x-1=0的根

迭代求根方法

  1. 将方程变形,求出根的收敛表达示g(x),即x=∛(x+1)。而不能变成x=x³-1,因为这个x的表达式是发散的。
  2. 选一个方程的近似根2,赋给变量x0。
  3. 将x0的值保存于变量x1,然后计算g(x1),并将结果存于变量x0;
  4. 当x0与x1的差的绝对值还大于指定的精度要求时,重复步骤2、3的计算。

代码

cpp 复制代码
#include<stdio.h>
#include<math.h>
typedef double(*fpt)(double);//函数指针 
int rootX(fpt g,double initRoot,double precision,double *root){
	double x0=initRoot,x1=0;
	int n=0;//迭代次数
	do{
		printf("x%d=%.5f\n",n,x0);
		x1=x0;
		x0=g(x1);
		n++;
	} while(fabs(x0-x1)>= precision);
	*root=x0;
	printf("x%d=%.5f\n",n,x0);
	return n;
}
double f(double x){
	return pow(x+1,1.0/3.0);
}
int main(){
	double x;
	int m=rootX(f,2,0.000005,&x);
	printf("迭代%d次,根为%.5f\n",m,x);
}

说明

pow 函数

  • base 是底数,可以是任意浮点数。
  • exponent 是指数,也可以是任意浮点数。
  • 返回值是 base 的 exponent 次幂,类型为 double。
cpp 复制代码
double pow(double base, double exponent);
相关推荐
_不会dp不改名_15 小时前
leetcode_21 合并两个有序链表
算法·leetcode·链表
进击的_鹏15 小时前
【C++11】initializer_list列表初始化、右值引用和移动语义、可变参数模版等
开发语言·c++
mark-puls15 小时前
C语言打印爱心
c语言·开发语言·算法
Python技术极客15 小时前
将 Python 应用打包成 exe 软件,仅需一行代码搞定!
算法
tongsound15 小时前
igh ethercat 实时性测试
linux·c++
睡不醒的kun15 小时前
leetcode算法刷题的第三十四天
数据结构·c++·算法·leetcode·职场和发展·贪心算法·动态规划
吃着火锅x唱着歌15 小时前
LeetCode 978.最长湍流子数组
数据结构·算法·leetcode
晚云与城15 小时前
今日分享:C++ deque与priority_queue
开发语言·c++
我星期八休息16 小时前
深入理解跳表(Skip List):原理、实现与应用
开发语言·数据结构·人工智能·python·算法·list
lingran__16 小时前
速通ACM省铜第四天 赋源码(G-C-D, Unlucky!)
c++·算法