数值计算方法作业:迭代法求方程近似解

编程实现,用迭代求 f ( x ) = x 3 + 4 x 2 − 10 f(x)=x^3+4x^2-10 f(x)=x3+4x210在区间[1,2]内的一个实根,要精确到小数点后第三位为精确值的近似
注意:1。迭代格式的构建,2.迭代法的停止条件

首先先找出迭代函数,将 f ( x ) = 0 f(x)=0 f(x)=0转化为 x = 10 / ( x + 4 ) x=\sqrt{10/(x+4)} x=10/(x+4) ,然后进行迭代求解(转化迭代函数时要注意转化得到的迭代函数的区间要包含 f ( x ) f(x) f(x)的区间,否则在运算时会出错),具体代码如下:

C++代码

#include <bits/stdc++.h>
/****************
author:WANG_zibi
time:2020/2/27
title:迭代法求方程近似值
course:数值计算方法
****************/
using namespace std;
const int N = 1e6+5;
typedef long long ll;
typedef unsigned long long ull;
const int INF= 0x3f3f3f3f;
#define f(i, a, b) for (int i = (a); i <= (b); ++i)
const double eps = 1e-4;

double func(double x)
{
    return x*x*x+4*x*x-10.0;
}
double g(double x)
{
return sqrt(10.0/(x+4.0));
}
int main()
{
    double x0,x1,left,right;
    left=1,right=2;
	x0=g(left);
	x1= g(x0); 
	
	while(fabs(x0-x1)>eps)
	{
		x0=x1;
		x1=g(x0);
	 } 
	cout<<"答案是:"<<x1<<"\n" ;
    return 0;
}

MATLAB代码:

eps=0.0001;

x0=l;
x0=(10/(x0+4))^(1/2);
x1=(10/(x0+4))^(1/2);
while abs(x0-x1)>eps
    x0=x1;
    x1=(10/(x1+4))^(1/2);
end
root=x1
  • 5
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值