數(shù)值分析-常微分?jǐn)?shù)值解的求法c語言(已改無錯(cuò)字)
2023-07-19 04:39:23 本頁面
【正文】 , y0=1。 while(1) { printf(請(qǐng)輸入下限， 上限， 步長(zhǎng) ， 初值：(用空格隔開)\n)。 scanf(%lf %lf %lf %lf,amp。a,amp。b,amp。h,amp。y0)。 double c = (ba)/h。 double *yx,*y。 yx=Exact(a,b,h)。 while(1) { printf(歐拉 改進(jìn)歐拉 經(jīng)典KR4階Adams預(yù)測(cè)校正 修改數(shù)據(jù) 退出\n)。 scanf(%d,amp。selet)。 system(cls)。 if(selet==1) { y=Euler(a,b,h,y0)。 printf(歐拉方法：\n)。 } if(selet==2) { y=Euler_Pro(a,b,h,y0)。 printf(改進(jìn)歐拉方法：\n)。 } if(selet==3) { y=K_R(a,b,h,y0)。 printf(經(jīng)典KR方法(4階龍格庫(kù)塔方法)：\n)。 } if(selet==4) { double **r。 r=Adams(a,b,h,y0)。 printf(4階Adams預(yù)測(cè)校正方法:\n)。 printf(x值 預(yù)估值 校正值 誤差\n)。 printf(% %lf \n,0,r[0][0])。 printf(% %lf \n,r[1][0])。 printf(% %lf \n,r[2][0])。 for(i=3。i=10。i++) printf(% %f %f %\n,a+i*h,r[i][0],r[i][1],fabs(r[i][0]r[i][1]))。 } if(selet==1||selet==2||selet==3) { printf(x值 計(jì)算值 準(zhǔn)確值 誤差\n)。 for(i=0。i=10。i++) printf(% %lf %lf %\n,(a+i)*h,y[
范文總結(jié)相關(guān)推薦
鄂ICP備17016276號(hào)-1