fortran习题参考答案

工程分析程序设计上机作业(五)

模块化编程

上机目的：掌握内部例程、主程序、外部例程、模块等功能的使用方法。

1、有一个六边形，求其面积。为求面积，作了3条辅助线。如图所示：(提示，三角形面

* a b c

积=s(s-a)(s-b)(s-c)，其中s ，a、b、c为三个边长)。要求用内部函数来计算每个三角形的面积。

program mai n

implicit none

real :: A(9)=(/10,30,16,13,21,14,20,36,28/)

in teger i

real s

s=area(A (3) ,A (4) ,A(7))+area(A(2),A(7),A(8))+ &

area(A(1),A(8),A(9))+area(A (5),A (6) ,A(9))

prin t*,'The area is',s

contains

function area(x,y,z)

implicit none

real x,y,z

real p,area

p=(x+y+z)/2

area=sqrt(p*(p-x)*(p-y)*(p-z))

end function

end

2、编写一个子例程子程序，SUM(S , T , N1 , N2)。把整型数N1到N2进行求和，并把求

和的结果放置到S,把整型数N1到N2进行求积，并把求积的结果放置到T。并用这个子程序来计算：

y =(1 2 3 4) (3 4 5 6 7 8) (3*4*5*6) -(1*2*3)

program mai n

implicit none

in teger A(8)

in teger i,b,c,d,e,f,g,m, n,y

b=O;d=O;g=1; n=1 do i=1,8

A(i)=i

end do

call sum(b,c,1,4) call sum(d,e,3,8) call sum(f,g,3,6) call sum (m,n ,1,3)

y=b+d+g-n

prin t*,'y=(1+2+3+4)+(3+4+5+6+7+8)+(3*4*5*6)-(1*2*3)=',y

contains

subroutine sum(s,t,N1,N2) in teger s,t,N1,N2,i do i=N1,N2 s=s+A(i) t=t*A(i)

end do

end subrout ine

end

3、编写函数子程序GDC求两个数的最大公约数。求最大公约数的算法如下：把两个数中大的那个数

作为被除数，两数相除得到一个余数。把余数去除除数得到新一轮的余数。

不断重复这一过程直到余数为0,这时的除数就是两个数的最大公约数。

调用此函数，求1260, 198, 72三个数的最大公约数。program mai n

implicit none

in teger a,b,c, m,n ,temp

prin t*,'请输入三个数：’

prin t*,'a b c'

read*,a,b,c

if (a

temp=a

a=b

b=temp

end if

call GDC(a,b,m)

if (c

temp=c

c=m

m=temp

end if

call GDC(m,c, n)

print*,'最大公约数为',n

contains

subrouti ne GDC(x,y,z) in teger x,y,z,k do while(mod(x,y)/=0) k=mod(x,y) x=y y=k end do z=y

end subrout ine

end

n 4、编写一个模块程序，提供以下服务：定义出常量二、e。定义出子程序，实现求和、n2、

i=d 求阶乘n!。并在主程序中计算如下结果：从键盘上输入整数n、实型数A、R、R0，求! A f R " - — I

—一（实型）和一A—— e交。丿（实型）。"R 2

module mytask

implicit none

real,parameter :: Pl=3.1415926

real,parameter :: e =2.7182818

contains

function sub( n)

in teger i,n, sub

sub=0

do i=1, n

sub=sub+i*i

end do

end function sub

fun ctio n fact( n)

in teger i,n

real fact

fact=1

do i=1, n

fact=fact*i

end do

end fun ctio n fact

end module mytask

program mai n

use mytask

implicit none

in teger n

real A,R,R0

real x,y

print*,'请输入依次n,A,R,R0 :'

read*, n,A,R,R0

x=fact( n) /sub( n)

y=(A* n/(2*PI*R*R))*(R/R0)** n*e**(-(R/R0)** n)

print*,'n!/ 刀n A2(i=1,n)=',x

prin t*,'(A* n/(2*PI*R*R))*(R/R0)** n*e**(-(R/R0)** n)',y end program

10-6为止。再编写主程序，从键盘读入 意不能用Fortran 的标准函数SIN(X)。

program mai n

implicit none

real x,p

exter nal sub

print*,'请输入角度x :' read*,x

p=x/180*3.14

call sub(p)

end program subrouti ne sub(k) implicit none real k,y in teger i y=0 i=1

do while (abs(k**i/(fact(i)))>1.0e-6) y=y+(-1)**((i-1)/2)*k**i/(fact(i)) i=i+2

end do

prin t*,'si n(x) = ',y

contains

fun ctio n fact(m) in teger m,j,fact

fact=1 do j=1,m

fact=fact*j end

do 5、编写外部函数，计算 3 x 5 x 7 x

sin x = x 「 + ...的3! 5! 7! 直到最后一项的绝对值小于

x ,调用该外部函数并输出

sinx 的计算结果。注