#include//y'=x+ydoublef1(doublex,doubley){returnx+y;}//y'=3y/(1+x)doublef2(doublex,doubley){return3*y/(1+x);}//y'=y*y;doubleftest(doublex,doubley){returny*y;}voidsolve(double(*func)(doublex,doubley),doubleminX,doublemaxX,doubley0,doubleh,doubleresult[][7],int*resultNum){doubleK1,K2,K3,K4;doubleXn_1,Yn_1;intn=0;result[n][0]=n;result[n][1]=minX;result[n][2]=0;result[n][3]=0;result[n][4]=0;result[n][5]=0;result[n][6]=y0;for(n=1;n*h<=maxX;n++){Xn_1=result[n-1][1];Yn_1=result[n-1][6];K1=(*func)(Xn_1,Yn_1);K2=(*func)(Xn_1+h/2,Yn_1+h/2*K1);K3=(*func)(Xn_1+h/2,Yn_1+h/2*K2);K4=(*func)(Xn_1+h,Yn_1+h*K3);result[n][0]=n;result[n][1]=minX+n*h;result[n][2]=K1;result[n][3]=K2;result[n][4]=K3;result[n][5]=K4;result[n][6]=Yn_1+h*(K1+2*K2+2*K3+K4)/6;}*resultNum=n;}voidprint(doubleresult[][7],intresultNum){inti;doublex;printf("%5s%15s%15s%15s%15s%15s%15s\n","n","Xn","K1","K2","K3","k4","Y");for(i=0;i<95;i++)printf("-");printf("\n");for(i=0;i
上一篇:嶽陽商道互聯信息技術有限公司怎麽樣？下一篇:襄陽天潤智慧城多少錢壹平方？