修正的AB4-AM4解法

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define H 0.1
int main()
{
	float x,y,r=3;
	float k1,k2,k3,k4,t[100]={0},Y[100]={0},R[100]={0};
	float fun1(float,float);
	int i,j,n;
	printf("请输入x,y的初值：\n");
	scanf("%f%f",&x,&y);
	printf("请输入运算次数：\n");
	scanf("%d",&n);
	printf("运算结果为：\n"); 
    for(i=4;i<n;i++)
	{
		if(i<5)
		{
            for(j=0;j<4;j++)
            {
      		    printf("x%-2d=%.2f,\ty%-2d=%f,\ty(%-2d)=%f",j,x,j,y,j,r);
		        printf("\n"); 
	            k1=fun1(x,y);
	            k2=fun1(x+H/2,y+H*k1/2);
	            k3=fun1(x+H/2,y+H*k2/2);
	            k4=fun1(x+H,y+H*k3);
       			R[j]=y;	
	            x=x+H;
	            y=y+H*(k1+2*k2+2*k3+k4)/6;
	            r=3/(1+x*x*x);
            }
            x=x-H;
		}
		t[i]=R[i-1]+H*(55*fun1(x,R[i-1])-59*fun1(x-H,R[i-2])+37*fun1(x-2*H,R[i-3])-9*fun1(x-3*H,R[i-4]))/24;
		Y[i]=R[i-1]+H*(9*fun1(x+H,t[i])+19*fun1(x,R[i-1])-5*fun1(x-H,R[i-2])+fun1(x-2*H,R[i-3]))/24;
		R[i]=251.0/270*Y[i]+19.0/270*t[i];
		x=x+H;
		r=3/(1+x*x*x);
        printf("x%-2d=%.2f,\ty%-2d=%f,\ty(%-2d)=%f,\tr(%-2d)=%f",i,x,i,R[i],i,r,i,r-R[i]);
        printf("\n"); 
	}
	system("pause");
	return 0;
}
	float fun1(float x,float y)
	{
		float f;
		f=-x*x*y*y;
		return f;
	}