fortran習(xí)題參考答案

工程分析程序計(jì)上作業(yè)（五）模塊化編程上機(jī)目的：掌握內(nèi)部例程、主程序、外部例程、模塊等功能的使用方法。、有個(gè)六邊形，其面積。為求面積，作了3條助線。如圖所示示，角形面積

(s)(s)()

，其中s

a2

，abc為個(gè)長(zhǎng)用內(nèi)部函數(shù)來計(jì)算每個(gè)三角形的面積。l=10l=14l

ll

lllprogrammain

limplicit::A(9)=(/10,30,16,13,21,14,20,36,28/)iis',sfunctionimplicitp,areaarea=sqrt(p*(p-x)*(p-y)*(p-z))/

&

、編一個(gè)子例程程序，，N1N2)。把整型N1到行求和，并把求和的結(jié)果放置到S，整型數(shù)到N2進(jìn)行求積，并把求積的結(jié)果放置T并用這個(gè)子程序來計(jì)算：y4)(3*4*5*6)programmainimplicitb=0;d=0;g=1;n=1i=1,8A(i)=idocallcallcallcallprint*,'y=(1+2+3+4)+(3+4+5+6+7+8)+(3*4*5*6)-(1*2*3)=',ysum(s,t,N1,N2)s,t,N1,N2,it=t*A(i)do、編函數(shù)子程序求兩個(gè)數(shù)的最大公約數(shù)。求最大公約數(shù)的算法如下：把兩數(shù)中大的那個(gè)數(shù)作為被除數(shù)，兩數(shù)相除得到一個(gè)余數(shù)。把余數(shù)去除除數(shù)得到新一輪的余數(shù)。不斷重復(fù)這一過程直到余數(shù)為0這時(shí)的除數(shù)就是兩個(gè)數(shù)的最大公約數(shù)。/

調(diào)用此函數(shù)，求1260，，72個(gè)數(shù)的最大公約數(shù)。programmainimplicitprint*,'請(qǐng)輸入三個(gè)數(shù)：c'read*,a,b,cif(a<b)thentemp=aifcallGDC(a,b,m)if(c<m)c=mifcallGDC(m,c,n)print*,'最大公約數(shù)為GDC(x,y,z)x,y,z,kwhile(mod(x,y)/=0)x=yy=kdo、編一個(gè)模塊程，提供以下服務(wù)：定義出常定出程序，實(shí)現(xiàn)求和

n

2

、i求階乘n!。并在主程序中計(jì)算如下結(jié)果：從盤上輸入整數(shù)n、實(shí)型數(shù)、、R，/

n2iR0n2iR0n!n

（實(shí)型）和

An

2

n

（實(shí)型modulemytaskimplicit::::functioni,n,subsub=0i=1,ndosubfunctionfact(n)i,ni=1,nfact=fact*idomytaskprogrammainimplicitnA,R,R0x,yprint*,'請(qǐng)輸入依次n,A,R,R0:'x=fact(n)/sub(n)print*,'n!/n^2(i=1,n)=',x/

、編外部函數(shù)，算

sinxx

x35!

的值，直到最后一項(xiàng)的絕對(duì)值小于

-6止。再編寫主程序，從鍵盤讀入x，用該外部函數(shù)并輸出sinx的算果。注意不能用Fortran的準(zhǔn)函數(shù)SIN(X)programmainimplicitprint*,'請(qǐng)輸入角度xp=x/180*3.14callsub(p)implicitk,yii=1while(abs(k**i/(fact(i)))>1.0e-6)y=y+(-1)**((i-1)/2)*k**i/(fact(i))i=i+2doprint*,'sin(x)=functionfact(m)j=1,mfact=fact*jdo/

hhsubroutinesub、編程序，用歐法數(shù)值求解一階微分方程。數(shù)值解法的基本思想：用差分方程代替微分方程，然后在若干個(gè)離散點(diǎn)上逐點(diǎn)求解差分方程得各離散點(diǎn)xx,…等函數(shù)0

dydx

f()

的近似值y,yy,…其02各離散點(diǎn)x,xx…間的距離稱為步長(zhǎng)。0dy歐拉法的基本思想是，一階微分可用向前差分代替，即dx分方程，可得ynn(xy)nh

1y'()n

。帶入微因此，

y

y*fx,n

，其中

。當(dāng)給定x,y和長(zhǎng)后即可按下面步驟求得數(shù)值解：0xy

(y1

,(y)1

f

y

請(qǐng)按上述方法編寫程序用Euler法解微分方程

2

2

，當(dāng)x=0時(shí)y。取00=0.1試求出當(dāng)x=0.1,0.2,0.3,0.4iiprogrammainimpliciti,ji=1,10x(i)=x(i-1)+0.