euler预测校正法解微分方程.txt

数值分析程序

#include <iostream.h>
#include <math.h>

double F(double x,double y)
{
	return (y-(2*x/y));
}

void EulerF_V(double *x,double *y,int n,double h)
{
	double ty;
	int i;
	for(i=1;i<=n;i++)
	{
		ty=y[i-1]+h*F(x[i-1],y[i-1]);
		y[i]=y[i-1]+h*F(x[i],ty);
	}
	for(i=0;i<=n;i++)
		cout<<y[i]<<endl;
}

void main()
{
	double a,b,*x,*y,h;
	int n,i;

	cout<<"a,b:"<<endl;    
	cin>>a>>b;

	cout<<"n"<<endl;
	cin>>n;

	h=(b-a)/n;

	x=new double[n+1];
	y=new double[n+1];
	
	x[0]=a;
	cout<<"y0"<<endl;
	cin>>y[0];

	for(i=1;i<=n;i++)
		x[i]=x[i-1]+h;

	EulerF_V(x,y,n,h);

	delete[] x;
	delete[] y;
}