1.cpp

来自「用C++写改进的欧拉法 计算方法上机实验」· C++ 代码 · 共 19 行

#include<iostream.h>
#include<math.h>
#define f(x,y) (pow(x,3)+pow(y,3)+1)
void main()
{   
  double x0,y0,h,n,j,k,t;
   cout<<"输入x0,y(x0):";
   cin>>x0>>y0;
   cout<<"输入步长h,计算次数n:";
   cin>>h>>n;
   j=x0,k=y0;
   for(int i=0;i<n;i++)
   {
	   t=k+h*f(j,k);
	   k+=(f(j,k)+f(j+h,t))*h/2;
	   j=j+h;
	   cout<<"y"<<i+1<<"="<<k<<endl;
   }
}