亚当姆斯预报校正系统2.cpp

数值分析实验内容一些计算源程序 数值分析实验内容一些计算源程序

#include <iostream>
using namespace std;
float f(float x,float y)  
{
 return y-2*x/y;
}

void main()
{
 float x0,y0,h,**result;
 int N;
 float f(float x,float y);
 void adams(float,float,int,float,float **result);
 cout<<"输入初始点：";
 cin>>x0>>y0;
 cout<<"输入解的个数和步长：";
 cin>>N>>h;
 result=new float *[N];
 for(int i=0;i<N;i++)
  result[i]=new float[2];
 adams(x0,y0,N,h,result);
 cout<<"运算结果是："<<endl; 
 for(i=0;i<N;i++)
  cout<<result[i][0]<<"     "<<result[i][1]<<endl;
}

void adams(float x0,float y0,int N,float h,float **result) 
{
 float k1,k2,k3,k4;
 float yp,yp1;
 result[0][0]=x0;result[0][1]=y0;
 for(int n=1;n<N && n<4;x0+=h,y0=result[n++][1])  //经典格式求前三点
 {
  result[n][0]=x0+h;
  k1=f(x0,y0);
  k2=f(x0+h/2,y0+h*k1/2);
  k3=f(x0+h/2,y0+h*k2/2);
  k4=f(x0+h,y0+h*k3);
  result[n][1]=y0+h*(k1+2*k2+2*k3+k4)/6;
 }
 for(;n<N;x0+=h,y0=result[n++][1])
 {
  result[n][0]=x0+h;
  yp=result[n-1][1]+h*(55*f(result[n-1][0],result[n-1][1])-59*f(result[n-2][0],
     result[n-2][1])+37*f(result[n-3][0],result[n-3][1])-9*f(result[n-4][0],
     result[n-4][1]))/24;
  yp1=f(result[n][0],yp);
  result[n][1]=result[n-1][1]+h*(9*yp1+19*f(result[n-1][0],result[n-1][1])-5*f(
      result[n-2][0],result[n-2][1])+f(result[n-3][0],result[n-3][1]))/24;
 }
}