改进的欧拉法.cpp

来自「数值分析程序。。。。。好多程序下载 喜欢的朋友可以下 随时欢迎」· C++ 代码 · 共 29 行

#include<iostream.h>
double f(double x,double y);
void main()
{
  double x0,y0,h,y1,x1,yp,yc;
  int n=1;
  x0=0;
  y0=1;
  h=0.1;
  for(n=1;n<=10;n++)
  {
    x1=x0+h;
	yp=y0+h*f(x0,y0);
	yc=y0+h*f(x1,yp);
	y1=(yp+yc)/2;
	y0=y1;
    x0=x1;
	cout<<x1<<"             " <<y1<<endl;
  }
}
double f(double x,double y)
{
  double  w;
  w=y-2*x/y;
  return w;
}