改进尤拉.cpp

来自「计算方法上机改进尤拉」· C++ 代码 · 共 49 行

#include "stdlib.h"
#include<stdio.h>
#include<iostream.h>
double f(double x0,double y0)
{
	double y;
	    y=y0-2*x0/y0;
	          return y;
}
void imp(double x0,double y0,double h,int N)
{
 int n,flag;
   double x1,y1,yp,yc;
     n=1,flag=0;
 cout<<"结果是:"<<endl;
   cout<<" xn       yn    "<<endl;
  do{
	  if(flag)
	  {
		n+=1;
		 x0=x1;
		  y0=y1;
	  }
	   x1=x0+h;
	     yp=y0+h*f(x0,y0);
	       yc=y0+h*f(x1,yp);
	         y1=(yp+yc)/2;
	   printf(" %.1f    %.4f\n",x1,y1);
	        flag=1;
  }while(n!=N);
}
int main(int argc, char* argv[])
{
	int N;
 double x0,y0,h;
  cout<<"请输入x0=";
    cin>>x0;
  cout<<"请输入y0=";
    cin>>y0;
  cout<<"请输入h=";
     cin>>h;
cout<<"请输入N=";
     cin>>N;		 
  imp(x0,y0,h,N);	
	return 0;
}