c语言求微分方程库函数 c语言如何求微分

C语言编程，常微分方程初值问题的Taylor求解方法的研究与实现。拜托大

#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;i95;i++)printf("-");printf("\n");for(i=0;i

创新互联公司专注于企业全网整合营销推广、网站重做改版、官渡网站定制设计、自适应品牌网站建设、H5网站设计、成都做商城网站、集团公司官网建设、外贸网站制作、高端网站制作、响应式网页设计等建站业务，价格优惠性价比高，为官渡等各大城市提供网站开发制作服务。