2022年fortran習(xí)題參考答案

1、fortran 習(xí)題參考答案 工程分析程序設(shè)計(jì) 上機(jī)作業(yè) 五 模塊化編程 上機(jī)目的 :把握內(nèi)部例程,主程序,外部例程,模塊等功能的使用方法； 1, 有一個(gè)六邊形 ,求其面積；為求面積 ,作了 3 條幫忙線；如以下圖 : 提示 , 三角形面積 = ss a s bs c ,其中 s ab c ,a,b,c 為三個(gè)邊長 ；要求 用內(nèi)部函數(shù)來計(jì) 2 算每個(gè)三角形的面積； l2=30 l 1=10 l 8=36 l 6=14 l 5=21 l 9=28 l 7=20 l 3=16 l 4=13 program main implicit none real : A9=/10,30,16,13,21,1

2、4,20,36,28/ integer i real s s=areaA3,A4,A7+areaA2,A7,A8+ & areaA1,A8,A9+areaA5,A6,A9 print*,The area is,s contains function areax,y,z implicit none real x,y,z real p,area p=x+y+z/2 area=sqrtp*p-x*p-y*p-z end function end 2, 編寫一個(gè)子例程子程序 ,SUMS,T,N1,N2 ；把整型數(shù) N1 到 N2 進(jìn)行求與 ,并把求與的結(jié)果 放置到 S,把整型數(shù) N1 到 N2 進(jìn)行求

3、積 ,并把求積的結(jié)果放置到 T；并用這個(gè)子程序來計(jì) 第 1 頁,共 6 頁fortran 習(xí)題參考答案 算: y 1 2 3 4 3 4567 8 3*4*5*6 1*2*3 program main implicit none integer A8 integer i,b,c,d,e,f,g,m,n,y b=0;d=0;g=1;n=1 do i=1,8 Ai=i end do call sumb,c,1,4 call sumd,e,3,8 call sumf,g,3,6 call summ,n,1,3 y=b+d+g-n print*,y=1+2+3+4+3+4+5+6+7+8+3*4*5*

4、6-1*2*3=,y contains subroutine sums,t,N1,N2 integer s,t,N1,N2,i do i=N1,N2 s=s+Ai t=t*Ai end do end subroutine end 3, 編寫函數(shù)子程序 GDC 求兩個(gè)數(shù)的最大公約數(shù)；求最大公約數(shù)的算法如下 :把兩個(gè)數(shù)中大 的那個(gè)數(shù)作為被除數(shù) ,兩數(shù)相除得到一個(gè)余數(shù)； 把余數(shù)去除除數(shù)得到新一輪的余數(shù)； 不斷 重復(fù)這一過程直到余數(shù)為 0,這時(shí)的除數(shù)就就是兩個(gè)數(shù)的最大公約數(shù)； 調(diào)用此函數(shù) ,求 1260,198,72 三個(gè)數(shù)的最大公約數(shù)； program main implicit none inte

5、ger a,b,c,m,n,temp print*, 請輸入三個(gè)數(shù) : print*,a bc read*,a,b,c if ab then temp=a a=b b=temp end if 第 2 頁,共 6 頁fortran 習(xí)題參考答案 call GDCa,b,m if c1 , 0e-6y=y+-1*i-1/2*k*i/facti i=i+2 end do print*,sinx = ,y contains function factm integer m,j,fact fact=1 do j=1,m fact=fact*j end do end function fact end s

6、ubroutine sub 6, 編寫程序 ,用歐拉法數(shù)值求解一階微分方程； 數(shù)值解法的基本思想 :用差分方程代替微分方程 ,然后在如干個(gè)離散點(diǎn)上逐點(diǎn)求解差分方 程,得到各離散點(diǎn) x0, x1, x2, 等處函數(shù) dy f x, y x 的近似值 y0, y1, y2, ；其中各 dx 離散點(diǎn) x0, x1, x2, 之間的距離稱為步長 h； 歐拉法的基本思想就是 ,一階微分 dy 可用向前差分代替 ,即 y 1 y n 1 y ；帶入微分 ndx h 方程 ,可得 yn 1 yn f x , y h因此 , yn 1 yn h * f xn , yn ,其中 h xn 1 xn ； 當(dāng)給定

7、x0, y0 與步長 h 后 ,即可按下面步驟求得數(shù)值解 : x0 y0 f x0 , y0 y1 x1 , f x1, y1 y2 x 2 f x2 , y y n請按上述方法編寫程序用 Euler 法求解微分方程 dy dx y 2 x 2 ,當(dāng) x0 =0 時(shí) ,y0 =1,0；取 h=0, 1,試求出當(dāng) xi=0, 1,0, 2,0, 3,0, 4 , ,1, 0 時(shí)的 yi 值； 第 5 頁,共 6 頁fortran 習(xí)題參考答案 program main implicit none real x0:10 real y0:10 real p0:10 integer2 i,j y0=1 , 0; x0=0 , 0 p0=fx0,y0 do i=1,10 xi=xi-1+0 , 1 end do do j=1,10 yj=yj-1+0 , 1*pj-1