AB4-AM4解法

来源：互联网发布：qq飞车暗风暴数据编辑：程序博客网时间：2024/04/29 17:19
#include<stdio.h>#include<stdlib.h>#include<math.h>#define H 0.1int main(){float x,y,r=3;float k1,k2,k3,k4,t[100]={0},Y[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);       Y[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]=Y[i-1]+H*(55*fun1(x,Y[i-1])-59*fun1(x-H,Y[i-2])+37*fun1(x-2*H,Y[i-3])-9*fun1(x-3*H,Y[i-4]))/24;Y[i]=Y[i-1]+H*(9*fun1(x+H,t[i])+19*fun1(x,Y[i-1])-5*fun1(x-H,Y[i-2])+fun1(x-2*H,Y[i-3]))/24;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,Y[i],i,r,i,r-Y[i]);        printf("\n"); }system("pause");return 0;}float fun1(float x,float y){float f;f=-x*x*y*y;return f;}
0 0