Eular方法解常微分方程
来源:互联网 发布:单机游戏网络对战平台 编辑:程序博客网 时间:2024/04/27 14:40
/**
***Eular方法***
Eular公式:Y(n+1)=Y(n) + h * f( x(n) , y(n) )
h=x(n+1)-x(n)
精度(误差):I= Y(x(n+1)) - Y(n+1) = h*h/2 Y''(ε) 局部截断误差
由泰勒公式可知Eular方法仅为一阶方法,是很粗糙的计算方法
属性:差分方法
《数值分析简明教程》-2 Editon -高等教育出版社- page 97 无算法流程图
代码维护:2005.6.14 DragonLord
**/
#include<iostream.h>
#include<stdio.h>
#include<math.h>
/*
举例方程:
y'= y - 2*x / y ( 0<x<1 )
y(0) = 1
*/
double f(double x,double y,double h)//方程形式
{
double re;
if(x==0)re=1;
else
re=y-2*x/y;
return re;
}
int main()
{
double x0,y0,x1,y1,h;//h为步长
double y;
int N;
while(cin>>x0>>y0>>h>>N)
{
int n=0;
for(;n<N;n++)
{
x1=x0+h;
y=sqrt(1+2*x1); //精确值
y1=y0+h*f(x0,y0,h);
printf("%.1f %.4f %.4f/n",x1,y1,y);
x0=x1;
y0=y1;
}
}
return 0;
}
- Eular方法解常微分方程
- 改进Eular方法解常微分方程
- Matlab解常微分方程
- 常微分方程数值解上机
- 常微分方程
- 常微分方程
- 用 GSL 解常微分方程初值问题
- 用 GSL 解常微分方程初值问题
- 常微分方程数值解:欧拉公式
- matlab求解常微分方程数值解
- 常微分方程实验(3.2):解的延拓
- Matlab基础学习-----------Matlab求解常微分方程边值问题的方法
- 常微分方程之差分法
- 线性常系数微分方程
- 常微分方程数值解法
- 【数值分析】常微分方程数值解:欧拉公式
- 【数值分析】常微分方程数值解:欧拉公式
- 【数学建模知识】常微分方程数值解:欧拉公式
- 语法分析表产生器的代码之三:自定义类的实现函数文件implement_syntax.cpp
- 测试SLR文法例子
- C++高手进阶
- 预计去南京
- 初学C#两星期,一片茫然
- Eular方法解常微分方程
- 富人和穷人的八大差异
- C++数据类型转换
- 改进Eular方法解常微分方程
- PING命令 前辈们说。首次中国红客攻击印尼的时候用的就是PING
- 人生最重要的十个健康伴侣
- 四阶Runge-Kutta法解常微分方程
- 床位的摆放有五大禁忌
- ASP.net中md5加密码的方法(转贴)