fortran习题3答案

工程分析程序设计上机作业（四）

数组

上机目的：练习数组的声明、存储、操作，以及数组参数、动态数组的使用。

1、请声明一个大小为10的一维数组，它们的初值为A(1)=2, A(2)=4, A(3)=6, ……A(I)=2*I,

并计算数组中这10个数字的平均值。

program main

implicit none

integer A(10)

integer i,sum

sum=0

do i=1,10

A(i)=2*i

sum=sum+A(i)

end do

print*,'The average of the ten numbers is:',sum/10.

end

2、编写一个程序来计算费氏数列的前10项，并把它们按顺序保存在一个一维数组当中。

费氏数列(Fibonacci Sequence)的数列规则如下：

F(0)=0

F(1)=1

当n>1时

F(n)=f(n-1)+f(n-2)

program main

implicit none

integer i

integer(2) A(11)

A(1)=0; A(2)=1

do i=3,11

A(i)=A(i-1)+A(i-2)

end do

print*,A(2:11)

End

3、输入任意n个数存放在数组中（如5个数1、2、8、2、10），请在屏幕上打印如下方阵

1 2 8 2 10

10 1 2 8 2

2 10 1 2 8

8 2 10 1 2

2 8 2 10 1

program main

implicit none

integer ,allocatable ::A(:)

integer(2) i,j,n,temp

print*,'请输入数字的个数：'

read*,n

allocate(A(n))

print*,'请输入数字：'

read*,A(1:n)

print*,A

do i=2,n

temp=A(1)

do j=1,n-1

A(j)=A(j+1)

end do

A(n)=temp

print*,A

end do

deallocate(A)

end

4、打印杨辉三角形，打印的行数由键盘输入。1

1 1

1 2 1

1 3 3 1

1 4 6 4 1

program main

implicit none

integer(2),allocatable::A(:,:)

integer i,j,n

print*,'请输入要打印的行数：'

read*,n

allocate(A(n,n))

do i=1,n

A(i,1)=1;A(i,i)=1

end do

if(n>=2) then

do i=1,n-1

do j=1,i

A(i+1,j+1)=A(i,j)+A(i,j+1)

end do

end do

end if

do i=1,n

print*,A(i,1:i)

end do

deallocate(A)

end

5、用“冒泡算法”对一个数列A(n)进行排序：

若要排序的数有n个，则需要n-1轮排序。第j轮排序中，从第一个数开始，相邻两数比较，若不符合所要求的顺序，则交换两者的位置；直到第n-j个数为止，第一个数与第二个数比较，第二个数与第三个数比较，......，第n-j-1个与第n-j个比较，共比较n-1次。此时第n-j个位置上的数已经按要求排好，所以不参加以后的比较和交换操作。例如：第一轮排序：第一个数与第二个数进行比较，若不符合要求的顺序，则交换两者的位置，否则继续进行二个数与第三个数比较......。直到完成第n-1个数与第n个数的比较。此时第n个位置上的数已经按要求排好，它不参与以后的比较和交换操作；第二轮排序：第一个数与第二个数进行比较，......直到完成第n-2个数与第n-1个数的比较；......第n-1轮排序：第一个数与第二个数进行比较，若符合所要求的顺序，则结束冒泡法排序；若不符合要求的顺序，则交换两者的位置，然后结束冒泡法排序。

共n-1轮排序处理，第j轮进行n-j次比较。

算法描述：

如果共有n个数：

第1个数要进行n-1次两两比较

第2个数要进行n-2次两两比较

第j个数要进行n-j次两两比较，j=1, n-j

n个数总共要进行n-1次排序

DO I = 1, N-1

DO J = 1, N-J

如果A(J) > A(J+1) 交换A(J) 和A(J+1)

END DO

END DO

program main

implicit none

integer ,allocatable ::A(:)

integer(2) i,j,k,n,temp

print*,'请输入数字的个数：'

read*,n

allocate(A(n))

print*,'请输入数字：'

read*,A(1:n)

do i=1,n-1

do j=1,n-1

if (A(j)>A(j+1)) then

temp=A(j)

A(j)=A(j+1)

A(j+1)=temp

end if

end do

end do

print*,A

deallocate(A)

end

6、从A、B两个数列中，把同时出现在两个数列中的数据删去。例如：

A：2 5 5 8 9 12 18

B：5 8 12 12 14

操作完成后：

A：2 9 18

B：14

program main

implicit none

integer(2) A(7),B(5),Temp(7)

external Sub

data A /2,5,5,8,9,12,18/

data B /5,8,12,12,14/

integer m,n

Temp=A

print*,'删去相同的数后：'

print*,'A'

call Sub(Temp,B,7,5)

print*,'B'

call Sub(B,A,5,7)

end

subroutine Sub(C,D,m,n) integer(2) C(m),D(n) integer i,j,k,m,n,p p=m 10 do i=1,p do j=1,n if (C(i)==D(j)) then do k=i,p-1 C(k)=C(k+1) end do p=p-1 goto 10 end if end do end do print*,C(1:p)

end subroutine

7、 高斯消去法求解线性代数方程组：对于一般的n 阶方程组，

??????

?=+++=+++=+++n

n nn n n n n n n b x a x a x a b x a x a x a b x a x a x a 22112222212111212111.....................................................

高斯消去法步骤如下：第一步：若011≠a ，令n i a a l i i ,...,3,2,/1111==，用()1i l -乘第1个方程加到第i 个方程上()n i ,...,3,2=，得同解方程组

(1)(1)(1)(1)11112211(2)(2)(2)

22222(2)(2)(2)22...................................n n n n n nn n n a x a x a x b a x a x b a x a x b ?+++=?++=??

??++=?

其中 .,,3,2,,,,3,2,,,,,,3,2,1,11)2(11)2(1)1(11)1(1n i b l b b n j i a l a a b b n j a a i i i j i ij ij j j =-==-==== 第二步：若,0)2(22≠a 令n i a a l i i ,,4,3,/)

2(22)2(22 ==，用)(2i l -乘第2个方程加到第i 个方程上

),,4,3(n i =，则将),,4,3()2(2n i a i =消去。一般，设第k-1步后方程组化为如下的同解方程组

(1)(1)(1)(1)11112211(2)(2)(2)

22222

()()

()()()()

............................................................................n n n n k k k kk k kn n k k k k nk n nn n n a x a x a x b a x a x b a x a x b a x a x b ?+++=?++=????++=?

??

?++=?

则第k 步：若若,0)(≠k kk a 令n k k i a a l k kk k ik ik ,,2,1,/)()( ++==，用)(ik l -乘第k 个方程加到第i

个方程上(,,,2,1n k k i ++=)，得到如下的同解方程组

????

??

???????=++=++=++=++=++++++++++++++++)

1()1(1)1(1,)

1(1)1(,11)1(1,1)()

()()

2(2

)2(22)2(22)1(1)1(12)1(121)1(11............................................................................k n n k nn k k k n k k n k n k k k k k k k n k kn k k kk n n n n b x a x a b x a x a b x a x a b x a x a b x a x a x a 其中(1)()()(1)()(),.1,,,,1,,k k k ij ij ik kj k k k i

i ik k a a l a i j k n b b l b i k n

++?=-=+??=-=+?? 按上述做法，做完n-1步，原方程组化为同解的上三角形方程组

??

??

?

?

?????==++=++=+++)

()()()

()()

2(2

)2(22)2(22)

1(1)1(12)1(121)1(11....................................................................n n n n nn k k n k kn k k kk n n n n b x a b x a x a b x a x a b x a x a x a 最后，设0)(≠n nn

a ，逐步代回得原方程组的解()()

()()1/,()/(1,2,,2,1).n n n n nn n

k k k

k k kj j kk

j k x b a x b a x a k n n =+?=??=-??

?=--?

∑

注意：上述公式中)()

(,k i k ij b a 的上标k ，是用来区别消去过程中第k 步利用的量。在用编程求解

时，可把)(k ij a 存在ij a 位置，)(k i b 存在i b 位置。

解方程组（1）?????-=++-=-+=++311874351321321

3

21321x x x x x x x x x （2）???=+=+-232121021215

x x x x

program main implicit none

real,allocatable ::A(:,:) real,allocatable ::X(:) integer(2) i,j,k,n real(4) m

print*,'请输入方程组所含方程的个数：' read*, n

此方程病态，提示：列主元消去法，第k 步

消去过程选取第k 行绝对值最大元素a kq ，交换k 和q 列，然后继续消去过程

allocate(A(n,n+1))

allocate(X(n))

print*,'请输入方程组的系数矩阵：'

do i=1,n

read*,A(i,:)

end do

do k=1,n-1

do i=k+1,n

A(i,:)=A(k,:)*(A(i,k)/A(k,k)*(-1))+A(i,:)

end do

end do

X(n)=A(n,n+1)/A(n,n)

do k=n-1,1,-1

m=0

do i=k+1,n

m=m + A(k,i)*X(i)

end do

X(k)=(A(k,n+1)-m)/A(k,k)

end do

do i=1,n

print*,'X(',i,')=',X(i)

end do

deallocate(A)

deallocate(X)

End

选作：

魔方阵，古代又称“纵横图”，是指组成元素为自然数1、2…n的平方的n×n的方阵，其中每个元素值都不相等，且每行、每列以及主、副对角线上各n个元素之和都相等。

如3×3的魔方阵：

8 1 6

3 5 7

4 9 2

5×5的魔方阵：

17 24 1 8 15

23 5 7 14 16

4 6 13 20 22

10 12 19 21 3

11 18 25 2 9

魔方阵的排列规律如下：

(1)将1放在第一行中间一列；

(2)从2开始直到n×n止各数依次按下列规则存放；每一个数存放的行比前一个数的行数减1，列数加1（例如上面的三阶魔方阵，5在4的上一行后一列）；

(3)如果上一个数的行数为1，则下一个数的行数为n(指最下一行);例如1在第一行，则2应放在最下一行，列数同样加1；

(4)当上一个数的列数为n时，下一个数的列数应为1，行数减去1。例如2在第3行最后一列，则3应放在第二行第一列；

(5)如果按上面规则确定的位置上已有数，或上一个数是第一行第n列时，则把下一个数放在上一个数的下面。例如按上面的规定，4应该放在第1行第2列，但该位置已经被占据，所以4就放在3的下面；

试打印出奇数阶魔方阵。

program main

integer n,i,x,y

integer(2),allocatable::a(:,:)

print*,'请输入魔方的阶数（奇数）'

read*,n

allocate(a(n,n))

a=0;y=(n+1)/2;x=1;a(x,y)=1

do i=2,n*n

if(x==1.and.y==n) then

x=2

else

x=x-1

if(x==0) x=n

y=y+1

if(y==n+1) y=1

if(a(x,y)/=0) then

x=x+1;y=y-1

if(x==n+1) x=1

if(y==0) y=n

x=x+1

end if

end if

a(x,y)=i

end do

do i=1,n

print*,a(i,1:n)

end do

end program

大学FORTRAN考试试题样题

《高级语言程序设计(FORTRAN)》考试样题 一、选择题(每小题2分,共40分) 1. 下列标识符中不能作为FORTRAN的合法变量名的是_____D_____。 A) A_1 B) A1 C) A1A D) A+1 2. 在FORTRAN的变量类型说明语句中没有 B A) REAL语句B)CHAR语句C)INTEGER语句 D)LOGICAL语句 3．FORTRAN90 源程序的扩展名是 C 。 A. TXT B. F C. F90 D. FOR 4.下列关于符号常量的说法正确的是___C____。 A.符号常量的定义应该放在所有语句之前 B.符号常量不可以是字符型常量 C.在一个程序单位中定义一个符号常量后，不能再改变它的值 D.符号常量定义语句与赋值语句一样，是可执行语句 5、在下列写法中，__D______不是FORTRAN常量。 A. .TRUE. B.1E-2 C. (0.0,2.4) D.1,000,000 6．表达式15／4／2.0的值是 B 。 A．整数2 B．实数1．5 C．实数2．25 D．实数1. 3 7. 下列不是FORTRAN赋值语句的是 _____D______。 A) X=.TRUE. B) X=X*X C) X=X**2**2 D) X=X++ 8. 梯形的上底为A,下底为B,高为H,计算其面积的FORTRAN表达式中错误的是____A____。 A)1/2*(AB)*H B)(A+B)*H/2 C)(A+B)/2*H D)0.5*(A+B)*H 9. 要选拔身高T>1.7米且体重W<62.5公斤的人,FORTRAN的逻辑表达式是__C_____。 A)T.ce.1.7.AND.w.le.62.5 B)T.le.1.7.OR.w.ge.62.5 C)T.gt.1.7. AND.w.lt.62.5D)T.gt.1.7.OR.w.lt.62.5 10．FORTRAN 90规定，变量类型声明的优先顺序是 C 。 A.隐含约定（I-N规则）、IMPLICIT声明、类型声明 B.类型声明、隐含约定（I-N规则）、IMPLICIT声明 C.类型声明、IMPLICIT声明、隐含约定（I-N规则） D.IMPLICIT声明、类型声明、隐含约定（I-N规则） 11、假设CH是一个字符型变量，并且具有足够的长度，则对于CH(3:7),下面说法正确的是___D________。 A.它是一个数组片段 B.它是一个数组定义形式 C.它是CH的一个子串，并且长度是4 D.它是CH的一个子串，并且长度是5 12．数组声明语句为： INTEGER,DIMENSION(-5:-1,-3:3,11:15) ::A 数组共有 A 个元素。 A．175 B．150 C．120 D．17 13. 下列DATA语句中正确的是______C____。 A) DATA /A,B/-1.0, -1.0/ B) DATA A, B/2*(-1.0)/ C) DATA A, B/2*-1.0/ D) DATA A/-1.0, B/-1.0 14. 下列有关FORTRAN数组的说法正确的是______D_____。 A) 数组定义语句只能定义一维或二维数组 B) 数组元素下标不能出现负值 C) 在赋值语句中出现的A(3)肯定是数组元素 D) 程序中使用的数组必须进行说明

Fortran95习题答案

第四章 1.program main implicit none write(*,*) "Have a good time." write(*,*) "That's not bad." write(*,*) '"Mary" isn''t my name.' end program 2.program main real, parameter :: PI=3 implicit none.14159 real radius write(*,*) "请输入半径长" read(*,*) radius write(*,"(' 面积='f8. 3)") radius*radius*PI end program 3.program main implicit none real grades write(*,*) "请输入成绩" read(*,*) grades write(*,"(' 调整后成绩为'f8.3)") SQRT(grades)*10.0 end program 4.integer a,b real ra,rb a=2 b=3 ra=2.0 rb=3.0 write(*,*) b/a ! 输出1, 因为使用整数计算, 小数部分会无条件舍去write(*,*) rb/ra ! 输出1.5 5.program main implicit none type distance real meter, inch, cm end type type(distance) :: d write(*,*) "请输入长度:" read(*,*) d%meter d%cm = d%meter*100 d%inch = d%cm/2.54 write(*,"(f8.3'米='f8.3'厘米='f8.3'英寸')") d%meter, d%cm, d%inch end program 第五章 1.program main

Fortran复习题讲解

FORTRAN程序设计复习题 选择题 B (1)下列各FORTRAN表达式中合法的是 A) S+T*2P >= 0.0 B) .NOT. (A*B+C) C) A2+B2/(C+D) <= 0.0 D) (A+B).NOT．A*B．GT．0.0(.NOT.只跟一个表达式) C (2)数学式(3/5)ex+y的FORTRAN表达式是 A) 3*EXP(X+Y)/5 B) 3*E* *(X+Y)/5.0 C) (3/5)*EXP(X+Y) D) (3.0/5.0)EXP(X+Y) D (3)下列FORTRAN77表达式中不合法的是 A) A．GT．B．EQV．C．GT．D B) A．AND．B．AND．C．AND．D C) .NOT．(X．LE．0.0) D) A．LT．B．LT．C．LT．D D (4)下列叙述中不正确的是 A) FORTRAN子程序可以单独编译 B) 对一个FORTRAN源程序进行编译和连接无误后可生成可执行文件 C) 即使编译和连接都正确无误，FORTRAN程序运行时仍可能出错 D) FORTRAN连接的主要任务是把函数库中的函数翻译成机器指令(主要任务为连接目标文件) B (5)在下列FORTRAN77运算符中，优先级最高的是 A) .AND. B) .NOT. C) .OR. D) .EQ. B (6)FORTRAN表达式"6/5+9/2**3/2"的值为 A) 33 B) 1 C) 5 D) 3 A (7)下列FORTRAN77 表达式中，合法的是：A）相当于0.0.LE.X .AND.（.NOT.（X.GE.10.0 ））A) 0.0.LE.X .AND..NOT.X.GE.10.0 B) 0.0.GE.X.LT.10.0 C) 0.0.LE.X.GE.10. D) 10.0.LE.X.OR..AND.X.GT.0.0 D (8)关于编译一个FORTRAN源程序文件，下列说法中错误的是 A) 允许编译只有一个主程序而没有子程序的源文件 B) 允许编译有多个子程序的源文件 C) 允许编译只有一个子程序而没有主程序的源文件 D) 允许编译有多个主程序的源文件 C (9)在FORTRAN77源程序中，续行标志符必须放在 A) 第1列 B) 第1-6列 C) 第6列 D) 第5列 D (10)下列关于"SUBROUTIN E MAP(X，Y)"语句行的叙述中，不正确的是 A) 这是子程序的第一个语句 B) 字符串"MAP"是子程序名 C) 变量X是子程序的形参 D) 子程序执行后，MAP将返回整型数据（子程序无返回值，自定义函数才有） A (11)FORTRAN表达式"2/4+0.5"的值是 A) 0.5 B) 1 C) 1.0 D) 0 D (12)FORTRAN表达式"MOD(4.9，1.3)"的值是 A)1.3 B)0.0 C)4.9 D)1.0 A (13下列FORTRAN运算符中，优先级最低的是 A）逻辑运算符.AND. B）算术运算符* C)关系运算符 >= D)算术运算符＋ A (14下列语句函数的定义中正确的是 A）F(X，Y)＝(X＋Y)/(X*Y)＋7.0 B）FUNCTION FUN(I，J，K)＝3*I＋2*J＋0.5*K

fortran习题4答案

工程分析程序设计 上机作业（五） 模块化编程 上机目的：掌握内部例程、主程序、外部例程、模块等功能的使用方法。 1、 有一个六边形，求其面积。为求面积，作了3条辅助线。如图所示：（提示，三角形面积=()()()s s a s b s c ---，其中2 a b c s ++=，a 、b 、c 为三个边长）。要求用内部函数来计算每个三角形的面积。 l 1=10 l 2=30 l 3=16 l 4=13 l 5=21 l 6=14 l 7=20 l 8=36 l 9=28 program main implicit none real :: A(9)=(/10,30,16,13,21,14,20,36,28/) integer 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)) print*,'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=++++++++++- (1234)(345678)(3*4*5*6)(1*2*3) program main implicit none integer A(8) integer i,b,c,d,e,f,g,m,n,y b=0;d=0;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 print*,'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) integer s,t,N1,N2,i do i=N1,N2 s=s+A(i) t=t*A(i) end do end subroutine end 3、编写函数子程序GDC求两个数的最大公约数。求最大公约数的算法如下：把两个数中 大的那个数作为被除数，两数相除得到一个余数。把余数去除除数得到新一轮的余数。 不断重复这一过程直到余数为0，这时的除数就是两个数的最大公约数。

大学FORTRAN考试试题1

1．下列叙述中，正确的是(D ) A．语句标号的大小影响程度执行的顺序 B．程序完全按语句出现的先后顺序执行 C．不同程序单位不能有相同的语句标号 D．同一程序单位不能有相同的语句标号 2．下列标识符中，不能作为合法的FORTRAN90标识符的是( C) A．A3_B3 B．VOID C．_123 D．IF 3．下列哪一个为正确的常量( D ) A．123,000 B. 3.5E+2.5 C. TRUE D. “HELLO” 4．若A=2，B=2，I=3，则表达式A**B**I的值为( C ) A. 64 B.12 C.256 D16. 5．圆的直径存放在整型变量D之中，下列计算圆面积的表达式中正确的是( D ) A. 3.14159*（D/2）*（D/2） B. 3.14159*(D*D/4) C. 3.14159*(D/2)**2 D. 3.14159*D*D/4 6．下列运算符中，运算优先级最高的是(B ) A. 关系运算 B. 算术运算 C. 逻辑非运算 D. 逻辑与运算 7．下列FORTRAN的表达式中值为0.5的是( C) A．MOD(4.8，0.5) B．100/20/10 C．50.0/4/25 D．MOD（15，10）/10 8．下列是完整的FORTRAN程序，编译时出错的语句是( C ) A．PROGRAM PRINT B．PARAMETER (PI=3.1415926) C．PI=PI+1 D．WRITE(*，*)SIN(PI+0.5) END 9．变量的类型定义中，优先级由高到低的顺序为( A ) A．类型说明语句、IMPLICIT 说明语句、隐含约定 B．隐含约定、IMPLICIT 说明语句、类型说明语句 C．IMPLICIT 说明语句、类型说明语句、隐含约定 D．隐含约定、类型说明语句、IMPLICIT 说明语句 10．“两整数M与N不能同时小于零”的FORTRAN表达式( D ) A．M＞＝0．AND．Ｎ＞＝０ B．(M．GE．0．AND．N．LT．0)．OR．(M．LT．0．AND．N．GE．0）C．(M．LE．0．AND．N．LE．0)．OR．(M．LE．0．AND．N．GE．0）D．M．GE．0．OR．N．GE．0

FORTRAN程序设计复习题及答案

FORTRAN程序设计复习题 一、选择题 B (1)下列各FORTRAN表达式中合法的是 A) S+T*2P >= B) .NOT. (A*B+C) C) A2+B2/(C+D) <= D) (A+B).NOT．A*B．GT．(.NOT.只跟一个表达式) C (2)数学式(3/5)ex+y的FORTRAN表达式是 A) 3*EXP(X+Y)/5 B) 3*E* *(X+Y)/ C) (3/5)*EXP(X+Y)D) EXP(X+Y) D (3)下列FORTRAN77表达式中不合法的是 A) A．GT．B．EQV．C．GT．D B) A．AND．B．AND．C．AND．D C) .NOT．(X．LE．D) A．LT．B．LT．C．LT．D D(4)下列叙述中不正确的是 A) FORTRAN子程序可以单独编译 B) 对一个FORTRAN源程序进行编译和连接无误后可生成可执行文件 C) 即使编译和连接都正确无误，FORTRAN程序运行时仍可能出错 D) FORTRAN连接的主要任务是把函数库中的函数翻译成机器指令(正确描述：主要任务为连接目标文件) B (5)在下列FORTRAN77运算符中，优先级最高的是 A) .AND. B) .NOT. C) .OR. D) .EQ. B (6)FORTRAN表达式"6/5+9/2**3/2"的值为 A) 33 B) 1 C) 5 D) 3 A (7)下列FORTRAN77表达式中，合法的是： A) .AND.. B) 10.0 C) D) 提示：A）相当于 .AND.（.NOT.（）） D (8)关于编译一个FORTRAN源程序文件，下列说法中错误的是 A) 允许编译只有一个主程序而没有子程序的源文件 B) 允许编译有多个子程序的源文件 C) 允许编译只有一个子程序而没有主程序的源文件 D) 允许编译有多个主程序的源文件 C (9)在FORTRAN77源程序中，续行标志符必须放在 A) 第1列 B) 第1-6列C) 第6列D) 第5列 D (10)下列关于"SUBROUTIN E MAP(X，Y)"语句行的叙述中，不正确的是 A) 这是子程序的第一个语句 B) 字符串"MAP"是子程序名 C) 变量X是子程序的形参D) 子程序执行后，MAP将返回整型数据 提示：子程序无返回值，自定义函数才有） A (11)FORTRAN表达式"2/4+"的值是 A) B) 1 C) D) 0 提示：2/4默认等于整型，=》 D (12)FORTRAN表达式"MOD，"的值是 A) B)0.0 C) D) A (13下列FORTRAN运算符中，优先级最低的是 A）逻辑运算符.AND. B）算术运算符*

FORTRAN77试题B

贵州大学2006-2007学年第一学期考试试卷 B 高级语言程序设计 注意事项： 1. 请考生按要求在试卷装订线内填写姓名、学号和年级专业。 2. 请仔细阅读各种题目的回答要求，在规定的位置填写答案。 3. 不要在试卷上乱写乱画，不要在装订线内填写无关的内容。 4. 满分100分，考试时间为120分钟。 一、选择题（共50分,每小题2分，各题中A）、B）、C）、D） 四个选项中，只有一个选项是正确的，请将正确的编号填在括号内） 1. 使用超大规模集成电路制造的计算机应该归属于（） A) 第一代B) 第二代C) 第三代D) 第四代 2. 下列FORTRAN77运算符中，优先级最低的是（）A）逻辑运算符.AND. B）算术运算符* C）关系运算符.GE. D）算术运算符＋ 3．下列标识符中，不属于FORTRAN77常量的是（）A）.TRUE. B）FALSE C）ˊROOT＝ˊ D）ˊ123.0ˊ 4.下列叙述中错误的是（）A）PROGRAM语句可以是FORTRAN77主程序的第一个语句 B）每一个FORTRAN77程序块的最后一个语句必须是END语句 C）FORTRAN77的所有语句都是可执行语句 D）FORTRAN77的任何可执行语句都可以有标号 5.下列FORTRAN77表达式中，值为.FALSE.的是（）A）.NOT.(.FALSE. .AND. .NOT. .TRUE.) B）.NOT. .FALSE. .OR. .NOT. .TRUE. C）.NOT. .TRUE. .AND. .NOT. .FALSE. D）X .GT. Y .OR. X .LE. Y 6. 在计算机内部，一切信息存取、处理和传送的形式是（） A) ASCII码B) BCD码C) 二进制D) 十六进制

fortran课后习题答案

第一章 FORTRAN程序设计基础第15页 1、2 1．简述程序设计的步骤。 “程序设计”：反映了利用计算机解决问题的全过程，通常要经过以下四个基本步骤：(1)分析问题，确定数学模型或方法；(2)设计算法，画出流程图；(3)选择编程工具，编写程序；(4)调试程序，分析输出结果。 2. 什么是算法？它有何特征？如何描述算法？ 解决问题的方法和步骤称为算法。 算法的五个特征：(1) 有穷性。 (2) 确定性。 (3) 有效性。 (4) 要有数据输入。(5) 要有结果输出。 算法的描述有许多方法，常用的有：自然语言、一般流程图、N－S图等。 第二章顺序结构程序设计 第29页 1、2、3、4、5、6、7、8、9 1．简述符号常量与变量的区别？ 符号常量在程序运行过程中其值不能改变。变量在程序运行过程中其值可以改变。 2. 下列符号中为合法的FORTRAN 90标识符的有哪些? (1) A123B (2) M%10 (3) X_C2 (4) 5YZ (5) X+Y (6) F(X) (7) COS(X) (8) A.2 (9) ‘A’ONE (10) U.S.S.R. (11) min*2 (12) PRINT 3. 下列数据中哪一些是合法的FORTRAN常量？ (1) 9,87 (2) .0 (3) 25.82(4) -356231 (5) 3.57*E2 (6) 3.57E2.1 (7) 3.57E+2(8) 3,57E-2 4. 已知A=2,B=3,C=5(REAL)；且I=2,J=3(INTEGER),求下列表达式的值: (1) A*B+C 表达式的值: 11 (2) A*(B+C) 表达式的值: 16 (3) B/C*A 表达式的值: 1.2 (4) B/(C*A) 表达式的值: 0.3 (5) A/I/J 表达式的值: 0.33 (6) I/J/A 表达式的值: 0 (7) A*B**I/A**J*2 表达式的值: 4.5 (8) C+(B/A)**3/B*2. 表达式的值: 7.25 (9) A**B**I 表达式的值: 512 5. 将下列数学表达式写成相应的FORTRAN表达式： (1) 1E-2 （2）(-B+SQRT(B*B-4*A*C)/(2*A) (3) 1+X+X*X/2+X**3/2/3 (4) COS(ATAN((A**3+B**3)**(1.0/3)/(C*C+1))) (5) EXP(A*X**2+B*X+C) (6) COS(X*Y/SQRT(X*X+Y*Y))**3 6. 用FORTRAN语句完成下列操作： (1) 将变量I的值增加1。I=I+1 (2) I的立方加上J,并将结果保存到I中。 I=I**3+J (3) 将E和F中大者存储到G中。G=Max(E,F) (4) 将两位自然数N的个位与十位互换，得到一个新的数存储到M中(不考虑个位为0的情况) M=MOD(N,10)*10+N/10 第三章选择结构程序设计第43页 1、2、3、5、6、7、9 1．分析下列程序运行结果 (1) LOGICAL P INTEGER I,I1,I2,I3 P=.FALSE. READ*,I I1=MOD(I,10) I2=MOD(I/10,10) I3=I/100

fortran习题参考答案

工程分析程序设计上机作业（五） 模块化编程 上机目的：掌握内部例程、主程序、外部例程、模块等功能的使用方法。 1、有一个六边形，求其面积。为求面积，作了3条辅助线。如图所示：（提示，三角形面 积 2 a b c s ++ =，a、b、c为三个边长）。要求用内部函 数来计算每个三角形的面积。 =21 program main implicit none real :: A(9)=(/10,30,16,13,21,14,20,36,28/) integer 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)) print*,'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=++++++++++- (1234)(345678)(3*4*5*6)(1*2*3) program main implicit none integer A(8) integer i,b,c,d,e,f,g,m,n,y b=0;d=0;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 print*,'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) integer s,t,N1,N2,i do i=N1,N2 s=s+A(i) t=t*A(i) end do end subroutine end 3、编写函数子程序GDC求两个数的最大公约数。求最大公约数的算法如下：把两个数中 大的那个数作为被除数，两数相除得到一个余数。把余数去除除数得到新一轮的余数。 不断重复这一过程直到余数为0，这时的除数就是两个数的最大公约数。

计算机技术基础(Fortran)试题库

一、选择题 1、下列关于"SUBROUTINE MAP(X，Y)"语句行的叙述中，不正确的是( ) A) 这是子程序的第一个语句B) 字符串"MAP"是子程序名 C) 变量X是子程序的形参D) 子程序执行后，MAP将返回整型数据 2、FORTRAN表达式"2/4+"的值是( ) A) B) 1 C) D) 0 3、阅读下列FORTRAN程序： PI=3. WRITE(*,' ')PI END 程序运行后输出结果是( ) A) B) 3.1415 C) D) 4、圆的直径存放在整型变量K之中，下列计算圆面积的表达式中正确的是( ) A) *K*K/4 B) *(K*K/4) C) *(K/2)**2 D) *(K/2)*(K/2) 5、要判断"月收入M在2000元以上(含2000元)且5000元以下)不含5000元)"的职工，应该使用的逻辑表达式是( ) A) M .GE. 2000 .AND. M .LE. 5000 B) M .GE. 2000 .OR. M .LE. 5000 C) M .GE. 2000 .AND. M .LT. 5000 D) M .GE. 2000 .OR. M .LT. 5000 6、设FORTRAN子程序如下： SUBROUTINE SUB(I，X) B=I+1 X=X+B END 若主程序中的变量M，B已按I-N规则赋值，则调用上述子程序的语句中正确的是( ) A) CALL SUB(M，M) B) CALL SUB(M+4，B) C) CALL SUB(B，B) D) CALL SUB(M，B+4) 7、下列不是FORTRAN常量的是( ) A) , B) +00 C) 2/3 D) 'Very good!' 8、阅读下列FORTRAN程序： F= DO 20 K=3, 7, 3 DO 20 L=K-1, K

FORTRAN90练习题

一．选择题 1．语句OPEN（3，FILE=‘ABC’）打开的是文件。 A．有格式顺序存取B．无格式顺序存取 C．有格式直接存取D．无格式直接存取 2．调用子程序时，假设实参可以是常数、内在函数及表达式，其相应形参是。 A．变量名B．数组名C．过程名D．* 3．下列数据中，不符合FORTRAN90常量表示的是。 A．-25.6 B．3.2344D0 C．’BOOT’’=’D．TRUE 4．圆心在原点的两个同心圆,半径分别为2和4。描述点(X,Y)在小圆外，大圆内(包括在两圆周上)的表达式是。 A．ABS（X）<=4.0 .AND. ABS(Y)>=2.0 B．2.0>=SQRT(X*X+Y*Y)<=4.0 C．X*X+Y*Y<=16.0 .AND. X*X+Y*Y>=4.0 D．（X.OR.Y）>=2.0 .AND. (X .OR.Y)<=4.0 5.下列字符串中可用作FORTRAN90变量名的是。 A．ABC B. A-BCD C. 2ABC D. AB.C 6.下列数据中不符合FORTRAN90常数表示的为。 A. 21.50 B. +0.0D0 C. .TRUE. D. $1000 7.表达式：－B+SQRT(B*B-4.0*A*C).LE.3E-5值的类型是。 A. 逻辑型 B. 整型 C. 实型 D. 字符型 8.按隐含约定，设X=2.5,Y=1.0,Z=3.0, 执行赋值语句I=X+Y/Z后I的值为。 A. 2.0 B. 2 C. 3 D. 3.0 9.下列语句中那个语句是错误的是。 A. READ(*,*)X B. READ*,X C. PRINT(*,*)X D. PRINT *,X 10.说明语句REAL A(1:2,-1:10)说明数组A的元素个数为。 A. 11 B. 12 C. 22 D. 24 11.M能被N整除的FORTRAN表达式为。 A. (M/M)=0 B. (M MOD N).EQ.0 C. (M/N).EQ.0 D. MOD(M,N).EQ.0 12．以下关于FORTRAN90中有关虚参、实参的说法中，不正确的是。 A．实参和虚参个数一定要相等B．实参和虚参个数可以不相等 C．虚参可以用数组 D．实参可以用数组 13. 下列关于FORTRAN 90源程序编辑规则的叙述之中,正确的是。 A. 任何以符号“!”开头的内容均是注释行。 B. 主程序第一行必须是PROGRAM语句。 C. 要实现有格式输入输出，必须有FORMAT语句。 D. 每个子程序中必须有一个RETURN语句。 14. 阅读下列程序 F=1.0 DO J=3,7,2 DO K=J-1,J F =F*K ENDDO ENDDO WRITE(*,'(I3,F8.1)') J,F

Fortran 上机题汇总

Fortran 复习题汇总program exer101_1 implicit none realpi,r,t,s C 说明符号常量 parameter(pi=3.1415926) C 打开用于保存结果的文件 open(9,file='mydata.dat') write(*,*)"请输入r和t的值：" read*,r,t s=t/360.0*pi*r**2 C 在屏幕上显示结果 print*,'r=',r,'t=',t,'s=',s C 采用格式说明语句显示结果 write(*,100)'r=',r,'t=',t,'s=',s 100 format(1X,3(A2,F10.5,1X)) C 将结果保存至文件中 write(9,*)'r=',r,'t=',t,'s=',s end program exer101_2 implicit none realpi.r,t,s parameter(pi=3.1415926) !说明符号常量 open(9,file='mydata.dat') !打开用于保存结果的文件write(*,*)"请输入r和t的值：" read*,r,t s=t/360.0*pi*r**2 print*,'r=',r,'t=',t,'s=',s !在屏幕上显示结果 write(*,100)'r=',r,'t=',t,'s=',s !采用格式说明语句显示结果100 format(1X,3(A2,F10.5,1X)) write(9,*)'r=',r,'t=',t,'s=',s !将结果保存至文件中 end 1

2 program ex202 implicit none real F,C print*,"请输入一个华氏温度：" read*,F

Fortran语言编程练习题

平面驻点流动方程数值解 1问题描述 自编程序完成平面驻点流动方程：的数值解。 边界条件为 式中：。 2求解过程 由于上述方程是非线性方程，可采用MATLAB 7.0软件来求解，步骤如下： 第一步：将方程化为一阶常微分方程组。 边界条件为 第二步：建立ode.m和lbc.m两个M文件。 ode.m文件程序如下： function dfdx=ode(x, f) dfdx=[f(2);f(3);-f(1)*f(3)+f(2)^2-1]; lbc.m 文件程序如下： function res=lbc(f0,finf) res=[f0(1);f0(2);finf(2)-1]; 第三步：求解方程。 在MATLAB 7.0工作窗口输入程序： infinity=4; solinit=bvpinit(0:0.4:infinity,[0 0 0]); sol=bvp4c(@ode,@lbc,solinit); x=0:0.4:infinity f=deval(sol,x) plot(x,f(1,:),'ob',x,f(2,:),'rp',x,f(3,:),'b*') /*绘图命令*/ xlabel('轴\it \eta');ylabel('轴\it \phi')

legend('平面驻点流动\phi曲线','平面驻点流动d\phi/d\eta曲线','平面 驻点流动d^2\phi/d\eta^2曲线') title('平面驻点流动的数值解') 3结果分析 平面驻点流动数值解的计算结果见表1，图1是平面驻点流动数值解的散点图。从表和图中可以看出，从开始呈线性增长，随着的增加，偏离斜直线，当以后渐近于1。在左右，，即此时粘性流动的流速已接近 势流流速，只差百分之一。 00.40000.8000 1.2000 1.6000 2.000 2.4000 2.8000 3.2000 00.08800.31240.62200.9798 1.3620 1.7553 2.1530 2.5523 00.41450.68590.84670.93240.97320.99060.99710.9992 1.23250.84630.52510.29380.14730.06580.02600.00910.0028 表1 平面驻点流动数值解的计算结果

fortran考试试题(A)有答案版

中南林业科技大学课程考试试卷 课程名称： Fortran 语言程序设计 ；试卷编号：A 卷；考试时间：120分钟 题号 一 二 三 四 总分 应得分 30 10 20 40 实得分 一、填空题（30分，每题2分） 得分 评卷人 复查人 1. Fortran 语言中有：整数型 逻辑型 字符型 实数型 复数型五种基本数据类型。 2. Fortran 程序的书写格式有两种固定格式 自由格式。 3. 字符型常数’this ”s ’的字符长度为6 4. 下列程序运行结果为125。 Program main Implicit none Integer::floor=5,i Do i=1, floor If(i==3)cycle If(i==4)cycle Write(*,*)i End do end 5. 下列程序的运行结果为7 6.0,2.0。 Program main Implicit none Real :: a, b, m A=2.0 B=76.0 If(a

sum=0.0 do i=1,203,2 sum=sum+i end do write(*,*) i end 7.下列程序的运行结果是12000。Program test Implicit none Integer::i Integer::a(5)=(/ (i, i=1,5) /) Integer::b(5)=0 Where (a<3) b=a End where Write(*,”(5(i3,1x))”) b end 8.下列程序的运行结果是30。Program test Implicit none Integer,parameter:: limit=10 Integer counter Integer:: ans=0 Counter=2 do while(counter<=limit) ans=ans+counter counter=counter+2 end do Wrte(*,*)ans end 9.下列程序的运行结果是2,3。Program main implicit none integer ::a=1 integer ::b=2 call add ( a ) Call add ( b) Write(*,*) a , b Stop end Subroutine add (sum) Implicit none

Fortran练习题

● 1.输出a，b，c的值，其中a=28,b=34,c=60.5 integer a,b,c定义整型integer a=28;b=34;c=60.5‘=’为赋值符号，把右赋给左 print*,'a=',a,'b=',b,'c=',c输出格式：print*,’文字内容’,该输出变量 end结束必须以end作结 ● 2.输出x和y的值 real x,y定义实型real print*,'输入x和y:' read*,x,y读入格式：read*,读入变量 print*,x,y end ● 3.求三个数的和及平均值 real sum,ave read*,a,b,c sum=a+b+c ave=sum/3 print*,'和为',sum print*,'平均值为',ave end ● 4.现有长1500m宽1000m的地块要拍卖建造住宅小区，地块中央有一半径为100m的 圆形公共绿地不属于拍卖范围，每平方米地价为0.5万元，编写程序计算该地块拍卖总价，并输出。 real area,money read*,a,b,r area=a*b-3.14*r**2a b为a**b money=area*0.5 print*,'总价为',money,'万元' print*,'地长为',a print*,'地宽为',b print*,'绿地半径为',r end ● 5.求三个数中的最大数与最小数 real max_x,min_x read*,a,b,c min_x=a if(bmax_x)then max_x=b

fortran 编程练习题

Fortran 程序设计 一．循环篇 1.编程找出并输出100-150之间和400-450之间能被9整除的数。 program main implicit none integer a,b a=100 do 10 while(a<=150) if(mod(a,9)==0) write(*,*)a 10 a=a+1 b=400 do 20 while(b<=450) if(mod(b,9)==0) write(*,*)b 20 b=b+1 end 2.回文数是指正读和反读都一样的数。如：232，编程求100到999之间的回文数。 program main implicit none integer a,b,c,n do 5 n=100,999 a=mod(n,10) b=n/100 5 n=n+1 If(a==b)then write(*,10)n 10format(1x,"100-999之间的回文数:",i3) end if stop end 3.输出所有水仙花数（水仙花数是指一个三位正整数，其各位数字的立方和等于该数本身）。 program main implicit none integer a,b,c,d do 5 d=100,999 a=mod(d/100,10) b=mod(d/10,10)

c=mod(d,10) if (d==a**3+b**3+c**3) then write(*,15) 15format(1x,i3,"为a水仙花数" ) else write(*,25) 25format(1x,i3,"不为水仙花数" ) endif 5 d=d+1 stop end 4 .编程输出如下图形。 * *** ***** ******* *** * program main implicit none integer i,j i=0 do 5 while(i.le.3) j=0 do 10 while (j.le.2-i) write(*,1) 1format(1x," ") 10 j=j+1 j=0 do 20 while(j.le.2*i) write(*,2) 2format(1x,"*") 20 j=j+1 5 i=i+1 i=0 do 15 while(i.ne.2) j=0 do 30 while(j.le.i) write(*,3) 3format(1x," ") 30 j=j+1

Fortran90试题

Fortran 复习 题目类型： 1、选择题 2、填空题 3、简答题 4、编程题 样题 一、选择题 1．下列名称中，符合FORTRAN 90命名规则的是。 A. a-b-c B. $123 C. length_1 D. date.year 2．关于FORTRAN 90程序单元，以下四种说法中正确的是。 A.一个FORTRAN 90程序由多个程序单元组成，其中有一个且只能有一个主程序单元 B. 一个FORTRAN 90程序由多个程序单元组成，其中有一个且只能有一个模块单元 C. 一个FORTRAN 90程序允许有多个程序单元，但每类程序单元只能有一个 D. 一个FORTRAN 90程序允许有多个程序单元，但主程序单元和模块单元只能有一个3．关于FORTRAN 90的CASE结构，以下说法中正确的是。 A. CASE结构内必须取DEFAULT作为情况选择器的最终值 B. CASE结构的情况表达式不能是字符类型 C. CASE结构的情况表达式和情况选择器的取值类型必须相同 D. CASE结构的不同情况选择器取值可以重复 4．已知数组说明REAL A(20:22,0:9,-3:-1)，数组A的数组元素个数是。 A．164 B．360 C．450 D．90 5．子程序的形式参数（虚拟变元或哑元）可以是。 A．变量名、数组名、过程名 B．变量名、数组名、过程名和符号常量名 C．变量名、数组名、数组片段、过程名 D．任意符号名

6. 语句OPEN（3，FILE=’ABC’，ACCESS=’DIRECT’，RECL=40）打开的是文件。 A．有格式顺序存取B．无格式顺序存取 C．有格式直接存取D．无格式直接存取 7、下列关于"SUBROUTINE MAP(X，Y)"语句行的叙述中，不正确的是( ) A) 这是子程序的第一个语句 B) 字符串"MAP"是子程序名 C) 变量X是子程序的形参D)子程序执行后，MAP将返回整型数据 8、 FORTRAN表达式"2/4+0.5"的值是( ) A) 0.5 B) 1 C) 1.0 D) 0 9、阅读下列FORTRAN程序： PI=3.14159265 WRITE(*,'(F7.4) ')PI END 程序运行后输出结果是( ) A) 3.142 B) 3.1415 C) 0.31416 D) 3.1416 10、圆的直径存放在整型变量K之中，下列计算圆面积的表达式中正确的是( ) A) 3.1415926*K*K/4 B) 3.1415926*(K*K/4) C) 3.1415926*(K/2)**2 D) 3.1415926*(K/2)*(K/2) 二、填空题 1、下列FORTRAN函数子程序的功能是用以下公式计算一组数据Z1,Z2,…,Zn的标准差σ： σ2=(Z12 +Z22+…+Zn2)/n-[(Z1+Z2+…+Zn)/n]2 请在程序中的下划线处填入合适的内容。 【1】 SIGMA(Z，N) REAL Z(N) SUMX=0.0 SUMX2=0.0 DO K=1,N 【2】 SUMX=SUMX+X SUMX2=SUMX2+X*X ENDDO AV=SUMX/N AV2=SUMX2/N S2=AV2-AV*AV