埃特金加速迭代求解方程
来源:互联网 发布:电路模拟软件下载 编辑:程序博客网 时间:2024/05/21 12:50
#include <iostream>
using namespace std;
void main()
{
float x0,result,e;
float u(float x);
float aitken(float (*f)(float),float x0,float e);
cout<<"输入初始值:";
cin>>x0;
cout<<"输入精度:";
cin>>e;
result=aitken(u,x0,e);
cout<<"运算结果是:"<<result<<endl;
}
float aitken(float (*u)(float),float x0,float e)
{
float x1,x2;
x1=u(x0),x2=u(x1);
x2=x2-(x2-x1)*(x2-x1)/(x2-2*x1+x0);
while(x2-x0>=e || x0-x2>=e)
{
x0=x2;
x1=u(x0);x2=u(x1);
x2=x2-(x2-x1)*(x2-x1)/(x2-2*x1+x0);
}
return x2;
}
float u(float x) [url=file://要/]file://要[/url]解的方程为x=u(x)
{
return x*x-12;
}
using namespace std;
void main()
{
float x0,result,e;
float u(float x);
float aitken(float (*f)(float),float x0,float e);
cout<<"输入初始值:";
cin>>x0;
cout<<"输入精度:";
cin>>e;
result=aitken(u,x0,e);
cout<<"运算结果是:"<<result<<endl;
}
float aitken(float (*u)(float),float x0,float e)
{
float x1,x2;
x1=u(x0),x2=u(x1);
x2=x2-(x2-x1)*(x2-x1)/(x2-2*x1+x0);
while(x2-x0>=e || x0-x2>=e)
{
x0=x2;
x1=u(x0);x2=u(x1);
x2=x2-(x2-x1)*(x2-x1)/(x2-2*x1+x0);
}
return x2;
}
float u(float x) [url=file://要/]file://要[/url]解的方程为x=u(x)
{
return x*x-12;
}
- 埃特金加速迭代求解方程
- 求解方程根的近似解:埃特金加速迭代法
- ccsu 1027 一元三次方程求解 (牛顿迭代)
- 二分逼近/牛顿迭代——一元高次非线性方程求解
- Codevs 1038 一元三次方程求解 NOIP 2001(导数 牛顿迭代)
- 迭代求解
- 特判方程迭代
- 求解方程
- 方程求解
- 方程求解
- 求解方程
- 求解方程
- 方程求解
- 方程求解
- 迭代求解fibnacci数列!
- G_S迭代求解线性方程
- 迭代求解矩阵测试
- PageRank迭代求解方法
- 最小二乘法
- html表单
- 约当消去法求解方程组
- 行列式求值
- 埃特金插值
- 埃特金加速迭代求解方程
- 经典格式式求解一阶微分的初值问题
- 改进的欧拉格式求解一阶微分的初值问题
- 变步长复化梯形法求积
- 亚当姆斯预报_校正系统解一阶微分的初值问题
- 优化浏览器渲染
- 快速弦截法求解方程
- 复化辛甫生求积
- 变步长复化一点高斯求积