改进欧拉公式.cpp

来自「计算方法的4个小程序,牛顿插值,lagrange插值,牛顿迭代,改进欧拉公式」· C++ 代码 · 共 29 行

#include<iostream.h>
double f(double x,double y)
	{
		return(-0.9*y)/(1+2*x);
	};
void main(void)
{
	double a,b,Y,h,yp=0;
	int n;
	cout<<"输入n"<<endl;
	cin>>n;
	cout<<"输入a"<<endl;
	cin>>a;
	cout<<"输入b"<<endl;
	cin>>b;
	cout<<"输入初始值y"<<endl;
	cin>>Y;
	h=(b-a)/n;
	cout<<"x[k]"<<'\t'<<"y[k]"<<endl<<endl;
	for(int i=0;i<=n;i++)
	{
		yp=(Y+i*h)+h*f((a+i*h),(Y+i*h));
		Y=Y+(h/2)*(f((a+i*h),(Y+i*h))+f((a+(i+1)*h),yp));
		cout<<a+i*h<<'\t'<<Y<<endl;
	}
	
}