【正文】 { printf(年份 :~%.0f,x+h)。 printf(計算人口 (千萬 ):%f\n,y1)。 printf( )。 break。 double h,x,y[20],yy,yp,m,y1,g[20]。 m=xm/y01。 31 k=0。 if(fabs(yp)) { g[i+1]=y[i]+h/2*(f(x+h,y[i],b2,l)+f(x,y[i+1],b2,l))。 printf( )。 if(fabs(yyy[i+1])esp) { g[i+1]=y[i]+h/2*(f(x+h,y[i],b2,l)+f(x,y[i+1],b2,l))。 printf( )。 } x=x+h。 32 x=a。kN。 k4=f(x+h,y+h*k3,b2,l)。 printf(年份 :~%.0f,x)。 printf(計算人口 (千萬 ):%f\n,y1)。 y[0]=y0。i10。 x=(i+1)*h。 if(fabs(y[i+1]p[i+1])=) c=1。 p[i+1]=y[i]+h/2*f(x,y[i],b2,l)。 printf(計算人口 (千萬 ):%f\n,y1)。 case 1: while(fabs(y[i+1]p[i+1])) { h=h*。 if(fabs(y[i+1]p[i+1])=) { printf(估計人口 (千萬 ):%f,y[i+1])。 } }。 cout請輸入 y0endl。 cout請輸入 k1endl。 cout請輸入 kendl。i=max1。 c[i]=1。 //cinb1[i]。 } for(i=2。 //coutendl。i++) { //cout請輸入 f[i]endl。 f[i]=r(x1)*h*h。 for(i=2。 for(i=2。 printf(利用追趕法求得結(jié)果如下 :\n)。i=1。 } } void main() { int a=0,b=100,N=10,max,d。 double c[100],a1[100],b1[100]。 printf( *1:改進歐拉法 *\n)。 printf( *5:經(jīng)典龍格 庫塔法 *\n)。 do { printf(請選擇進行人口預測所使用的算法（ 7）：\n)。 cinb2。 cinxm。 printf(請輸入該國的平均人口出生率： \n)。 printf(請輸入該國的最大人口容載量 (千萬 ):\n)。 } if(d==3) { p()。 ciny0。 ImplicitEulerMethord(N,a,b,y0,xm,b2,l)。 printf(請輸入該國十年前的人口數(shù)量 (千萬 )： \n)。 l=b2/xm。 cinb2。 cinxm。 printf(請輸入該國的平均人口出生率： \n)。 printf(請輸入該國的最大人口容載量 (千萬 ):\n)。 } if(d==7) { cout請輸入結(jié)點個數(shù) max= 。 } printf(是否繼續(xù)選擇其他算法 (Y/N):\n)。y39。 } 38 數(shù)理系課程設計評分表 課程名稱： 數(shù)值分析 項 目 評 價 設計方案的合理性與創(chuàng)造性 設計與調(diào)試結(jié)果 設計說明書的質(zhì)量 答辯陳述與回答問題情況 課程設計周表現(xiàn)情況 綜合成績 教師簽名： 日 期： 。Y39。 printf(\n)。 coutendl。 l=b2/xm。 printf(請輸入該國十年前的人口數(shù)量 (千萬 )： \n)。 R_KMethord(N,a,b,y0,xm,b2,l)。 ciny0。 } if(d==5) { p()。 printf(請輸入該國的最大人口 容載量 (千萬 ):\n)。 printf(請輸入該國的平均人口出生率： \n)。 cinxm。 cinb2。 l=b2/xm。 printf(請輸入該國十年前的人口數(shù)量 (千萬 )： \n)。 GJEulerMethord(N,a,b,y0,xm,b2,l)。 ciny0。 if(d==1) { p()。 printf( *7:線性有限差分法 *\n)。 printf( *3:隱式歐拉法 *\n)。 printf( 對 logistic 人口預測模型運用數(shù)值方法進行人口預測 :\n\n\n)。 double xm。 printf( )。 printf(x[%d]=%f\n,max1,x[max])。i++) { y[i]=(f[i1]a1[i]*y[i1])/(b1[i]a1[i]*n[i1])。i++) { n[i]=c[i]/(b1[i]a1[i]*n[i1])。 f[max1]=f[max1]y1。 //coutendl。 } for(i=0。i++) { //cout請輸入 a1[i]endl。 x1=h*i。i=max。 //cinc[i]。 h=(k1k2)/k。 cout請輸入 k2endl。 cout請輸入 y1endl。 } } } void FiniteDifferenceMethod(double f[],double a1[],double b1[],double 34 c[],int max)//有限差分法 { int i,k。 printf(計算人口 (千萬 ):%f\n,y1)。 p[i+1]=y[i]+h/2*f(x,y[i],b2,l)。 } }。 if(fabs(y[i+1]p[i+1])) { printf(估計人口 (千萬 ):%f,y[i+1])。 switch(c) { case 1: while(fabs(y[i+1]p[i+1])=) { h=h/。 printf(年份 :~%.0f,x)。 y[i+1]=y[i]+h*f(x,y[i],b2,l)。 m=xm/y01。 double h,x,y[20],p[20],y1,m。 printf(估計人口 (千萬 ):%f,y)。 x=x+h。 k2=f(x+h/2,y+(h*k1)/2,b2,l)。 m=xm/y1。 double h,x,k1,k2,k3,k4,y,y1,m。 break。 printf( )。 break。 printf( )。 yp=h*(b2*l*y[i+1])1。i10。 y[0]=y0。 } x=x+h。 printf( )。 } yy=y[i+1]yy/yp。 printf(估計人口 (千萬 ):%f,y[i+1])。 yp=h*(b2*l*y[i+1])1。i10。 y[0]=y0。 printf(計算人口 (千萬 ):%f\n,y1)。 printf(年份 :~%.0f,x)。i10。 x=a。 29 printf(計算人口 (千萬 ):%f\n,y1)。 printf(年份 :~%.0f,x)。 x=x+h。 m=xm/y1。 double h,x,t1,t2,y,y1,m。 } double r(double x) { return x。最終，問題解決了 通過這次課程設計使我懂得了理論與實際相結(jié)合是很重要的，使 自己的實際動手能力和獨立思考的能力能得到很大的提高。 在實驗過程中，我都是 嚴格按著老師的要求，先理清問題的解題思路，然后一步一步的翻譯成計算機語言程序，上機操作，不斷調(diào)試。 i++=max x[max]=y[max]。 fabs(y[i+1]p[i+1]) 是 否 fabs(y[i+1]p[i+1])= 輸 出 ：y[i+11],y1 是 否 i++=max 是 否 結(jié)束 23