常微分方程的查分方法——改进的欧拉方法

最新推荐文章于 2023-06-27 19:35:35 发布

飞花舞者

最新推荐文章于 2023-06-27 19:35:35 发布

阅读量247

点赞数

分类专栏：数值分析

本文链接：https://blog.csdn.net/qq_44975450/article/details/117230471

版权

常微分方程数值解改进的欧拉方法迭代 C++

关键词由CSDN通过智能技术生成

数值分析专栏收录该内容

12 篇文章 4 订阅

订阅专栏

程序流程图如下：
在这里插入图片描述

//常微分方程的差分方法——改进的欧拉方法——求解初值问题
 /************Analysis**************
 *1.h:步长，N：步数，迭代次数，x0,y0为老值，x1,y1为新值
 *f(x)=y-2x/y
 *y(0)=1
 **************************/
#include<iostream>
using namespace std;
double f(double x,double y)
{
	return y-2*x/y;
}
int main()
{
	double x0,y0;
	double h;//step
	unsigned N;
	unsigned n=0;
	cout<<"-------------Inproved Eular Method----------"<<endl;
	cout<<"Please input initial value of x0 :"<<endl;
	cin>>x0;
	cout<<"Please input initial value of y0:"<<endl;
	cin>>y0;
	cout<<"Please input step value h:"<<endl;
	cin>>h;
	cout<<"Please input iterative times H:"<<endl;
	cin>>N;
	cout<<"------------------Output------------------"<<endl;
	
	double x1,y1;
	double yp;
	double yc;
	
	n=1;
	while(true)
	{
		x1=x0+h;
		yp=y0+h*f(x0,y0);
	    yc=y0+h*f(x1,yp);
		y1=(yp+yc)/2;	
		cout<<"x1:"<<x1<<"    "<<"y1:"<<y1<<endl;

		if(n!=N)
		{
			n=n+1;
			x0=x1;
			y0=y1;
		}
		else
			break;
	}


	cout<<"Hello,Boker..."<<endl;
	system("pause");
	return 0;
}