fortran课后习题答案.doc

1、1第一章 FORTRAN 程序设计基础 第 15 页 1、21简述程序设计的步骤。“程序设计”：反映了利用计算机解决问题的全过程，通常要经过以下四个基本步骤：(1)分析问题，确定数学模型或方法；(2)设计算法，画出流程图；(3)选择编程工具，编写程序；(4)调试程序，分析输出结果。2. 什么是算法？它有何特征？如何描述算法？解决问题的方法和步骤称为算法。算法的五个特征：(1) 有穷性。 (2) 确定性。 (3) 有效性。 (4) 要有数据输入。(5) 要有结果输出。算法的描述有许多方法，常用的有：自然语言、一般流程图、NS 图等。第二章 顺序结构程序设计 第 29 页 1、2、3、4、5、6、

2、7、8、91简述符号常量与变量的区别？符号常量在程序运行过程中其值不能改变。变量在程序运行过程中其值可以改变。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) AONE (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

3、+2 (8) 3,57E-24. 已知 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 表达式的值: 5125. 将下列数学表达式写成相应的 FORTRAN 表达式：(1) 1

4、E-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)*36. 用 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 的情况)

5、M=MOD(N,10)*10+N/10第三章 选择结构程序设计 第 43 页 1、2、3、5、6、7、91分析下列程序运行结果(1)LOGICAL PINTEGER I,I1,I2,I3P=.FALSE.READ*,II1=MOD(I,10)I2=MOD(I/10,10)I3=I/1002IF(I1+I3.EQ.2*I2)P=.TRUE.PRINT*,PEND输入 123 时，输出： T输入 132 时，输出： F(2)REAL X,Y,ZREAD*,X,YIF(X*Y)0.AND.(X.NE.1.0)THENZ=10.0ELSE IF(X*Y0.AND.x=2.5.AND.y0.AND.b0

6、.AND.c0.AND.(a+b)c.AND.(b+c)a.AND.(c+a)b(8) a,b 之一为 0 但不同时为 0FORTRAN 表达式: (a=0.OR.b=0).AND.(.NOT.(a=0.AND.b=0)(9) 如果 m 与 n 相等,则 k 为 1,否则 k 为 0FORTRAN 语句:IF(m=n)THENk=1ELSEk=0ENDIF(10) a,b 是否相同号FORTRAN 表达式: a*b03.如果 A=2.5,B=7.5,C=5.0,D=6.0,L=.TRUE.,M=.FALSE.,请求出下列逻辑表达式的值。(1) (A+B).LT.(C+D).AND.A.EQ.3

7、.5 逻辑表达式的值: .FALSE.(2) A+B/2.0.NE.C-D.OR.C.NE.D 逻辑表达式的值: .TRUE.(3) .NOT.L.OR.C.EQ.D.AND.M 逻辑表达式的值: .FALSE.(4) C/2.0+D.LT.A.AND.NOT.TRUE.OR.C.EQ.D 逻辑表达式的值: .FALSE.(5) (C.GT.D).OR.C.NOT.(A+B.LT.D) 该表达式错误(6) (A.LT.B).AND.(B.LT.A) 逻辑表达式的值: .FALSE.5、某运输公司在计算运费时，按照运输距离 S 对运费打一定的折扣 D，其标准如下：S-15.AND.X=0.AND

8、.X15.AND.X1E-5)Y0=Y1Y1=2*Y0/3+X/(3*Y0*2)ENDDOPRINT*,Y1END6P.66 习题 12 Real a,b,m,xf(x)=x*3-x*2-1Read*,a,bDo while(Abs(a-b)1e-6)m=(a+b)/2If(f(a)*f(m)0)thena=mElseb=mEndifEnddoPrint*,(a+b)/2END第五章 FORTRAN 90 数据类型 第 88 页 4、104有一个三角形，顶点为 X,Y,Z，其坐标分别为(1.5,2.0),(4.5,4.5),(18.0,10.5)。求三角形的面积和重心（提示：重心坐标=（ X+

9、Y+Z）/3 利用复数计算） 。COMPLEX X,Y,Z,WREAL A,B,C,S,PREAD(*,*)X,Y,ZA=ABS(X-Y)B=ABS(Y-Z)C=ABS(Z-X)P=(A+B+C)/2S=SQRT(P*(P-A)*(P-B)*(P-C)W=(X+Y+Z)/3WRITE(*,*)SWRITE(*,*)WEND10设有 4 个圆，圆心分别在(2,2),(-2,-2),(2,-2),(-2,2),半径为 1。当 X,Y 坐标点落在圆上和圆内时， Z=1；落在圆外时,Z=1;落在圆外时，Z=0。编出相应的的程序（要求利用逻辑变量） 。IMPLICIT NONELOGICAL P,Q,R

10、,SREAL X,YINTEGER ZREAD*,X,YZ=0P= (X-2)*2+(Y-2)*2X(MA)MA=IENDDOT=X(1)X(1)=X(MA)X(MA)=TWRITE(*,10)(X(I),I=1,10)10 FORMAT(1X,10I4)END87. 求 34 矩阵中的最大元素和最小元素。REAL A(3,4),AMAX,AMININTEGER I,JDO I=1,3READ(*,*)(A(I,J),J=1,4)ENDDOAMAX=A(1,1)AMIN=A(1,1)DO I=1,3DO J=1,4IF(A(I,J)AMAX)AMAX=A(I,J)IF(A(I,J)1E8 时停

11、止计算。计算公式是：.!1.!21neFUNCTION JC(N)JC=1DO I=1,NJC=JC*IENDDOENDPROGRAM EX6REAL*8 E,T10E=0T=1.0I=0DO WHILE (T1E-8)E=E+TI=I+1T=1.0/JC(I)END DOPRINT*,EEND8对任意自然数 n，设计一个求 n 的各数位立方和的函数子程序 F(n)(如 F(121)=1*3+2*3+1*2=1+8+1=10)，并调用 F(n)求 1000 以内的所有水仙花数(水仙花数的概念请参看例 8.5)。!子程序FUNCTION F(N)INTEGER F,N,S,TS=0T=NDO WHILE(T0)K=MOD(T,10)S=S+K*3T=T/10ENDDOF=SEND!主程序PROGRAM EX8INTEGER F,IDO I=100,999IF(F(I)=I)PRINT*,IENDDOENDDO12设计一个子例行程序，计算一个任意数组 A 中正数、负数和 0 的个数及其和。！子程序SUBROUTINE TJ(A,M,P,N,Z,S)INTEGER A(M),M,P,N,Z,S,IP=0N=0Z=0S=0DO I=1,MIF(A(I)0)THENP=P+1ELSEIF(A(I)0)THENN=N+1