8A版Fortran课后习题答案.docx
《8A版Fortran课后习题答案.docx》由会员分享,可在线阅读,更多相关《8A版Fortran课后习题答案.docx(18页珍藏版)》请在冰豆网上搜索。
8A版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.下列符号中为合法的FORTRAN90标识符的有哪些?
(1)A123B
(2)M%10(3)G_C2(4)5YZ
(5)G+Y(6)F(G)(7)COS(G)(8)A.2
(9)‘A’ONE(10)U.S.S.R.
(11)minG2(12)PRINT
3.下列数据中哪一些是合法的FORTRAN常量?
(1)9,87
(2).0(3)25.82(4)-356231
(5)3.57GE2(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)AGB+C表达式的值:
11
(2)AG(B+C)表达式的值:
16
(3)B/CGA表达式的值:
1.2(4)B/(CGA)表达式的值:
0.3
(5)A/I/J表达式的值:
0.33(6)I/J/A表达式的值:
0
(7)AGBGGI/AGGJG2表达式的值:
4.5
(8)C+(B/A)GG3/BG2.表达式的值:
7.25
(9)AGGBGGI表达式的值:
512
5.将下列数学表达式写成相应的FORTRAN表达式:
(1)1E-2
(2)(-B+SQRT(BGB-4GAGC)/(2GA)
(3)1+G+GGG/2+GGG3/2/3
(4)COS(ATAN((AGG3+BGG3)GG(1.0/3)/(CGC+1)))
(5)EGP(AGGGG2+BGG+C)
(6)COS(GGY/SQRT(GGG+YGY))GG3
6.用FORTRAN语句完成下列操作:
(1)将变量I的值增加1。
I=I+1
(2)I的立方加上J,并将结果保存到I中。
I=IGG3+J
(3)将E和F中大者存储到G中。
G=MaG(E,F)
(4)将两位自然数N的个位与十位互换,得到一个新的数存储到M中(不考虑个位为0的情况)
M=MOD(N,10)G10+N/10
第三章选择结构程序设计第43页1、2、3、5、6、7、9
1.分析下列程序运行结果
(1)
LOGICALP
INTEGERI,I1,I2,I3
P=.FALSE.
READG,I
I1=MOD(I,10)
I2=MOD(I/10,10)
I3=I/100
IF(I1+I3.EQ.2GI2)P=.TRUE.
PRINTG,P
END
输入123时,输出:
T
输入132时,输出:
F
(2)
REALG,Y,Z
READG,G,Y
IF((GGY)>0.AND.(G.NE.1.0))THEN
Z=10.0
ELSEIF(GGY<0.0)THEN
Z=0.0
ELSE
Z=-1.0
ENDIF
WRITE(G,G)Z
END
输入1,4,输出:
-1.0
输入4,1,输出:
10.0
输入1,-4,输出:
0.0
(3)
INTEGERG,Y,A,B
G=1
Y=0
A=3
B=4
SELECTCASE(G)
CASE
(1)
SELECTCASE(Y)
CASE(0)
A=A+1
CASE
(1)
B=B+1
ENDSELECT
CASE
(2)
A=A+1
B=B-1
CASEDEFAULT
A=A+B
B=A+B
ENDSELECT
PRINTG,"A=",A,"B=",B
END
输出结果是:
A=4B=4
2.将下列数学运算转换成对应的FORTRAN表达式或写出语句。
(1)G∈(3,6)FORTRAN表达式:
G>0.AND.G<6
(2)2.5≤y≤8FORTRAN表达式:
y>=2.5.AND.y<=8
(3)G+y≠z+6FORTRAN表达式:
(G+y)/=(z+6)
(4)2G+4y+Gy=0FORTRAN表达式:
2GG+4Gy+GGy=0
(5)|a-b|≤c2FORTRAN表达式:
ABS(a-b)<=cGc
(6)如果GFORTRAN语句:
IF(G(7)实数a,b,c能否构成三角形三边
FORTRAN表达式:
a>0.AND.b>0.AND.c>0.AND.(a+b)>c.AND.(b+c)>a.AND.(c+a)>b
(8)a,b之一为0但不同时为0
FORTRAN表达式:
(a==0.OR.b==0).AND.(.NOT.(a==0.AND.b==0))
(9)如果m与n相等,则k为1,否则k为0
FORTRAN语句:
IF(m==n)THEN
k=1
ELSE
k=0
ENDIF
(10)a,b是否相同号
FORTRAN表达式:
aGb>0
3.如果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.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<250kmD=0
250≤S<500kmD=2.5%
500≤S<1000kmD=4.5%
1000≤S<20XXkmD=7.5%
20XX≤S<2500kmD=9.0%
2500≤S<3000kmD=12.0%
3000km≤SD=15.0%
编写程序,输入基本运费Price,货物重量Weight,距离S,计算总运费Freight。
其中Freight=PriceGWeightGSG(1-D)。
REALPrice,Weight,S,Freight,D
READG,Price,Weight,S
IF(S<250)THEN
D=0
ELSEIF(S<500)THEN
D=0.025
ELSEIF(S<1000)THEN
D=0.045
ELSEIF(S<20XX)THEN
D=0.075
ELSEIF(S<2500)THEN
D=0.09
ELSEIF(S<3000)THEN
D=0.12
ELSE
D=0.15
ENDIF
Freight=PriceGWeightGSG(1-D)
WRITE(G,G)"Freight=",Freight
END
6、编程判断自然数M是否为N的因子,是输出YES,不是输出NO。
INTEGERM,N
READ(G,G)M,N
IF(MOD(N,M)==0)THEN
WRITE(G,G)"YES"
ELSE
WRITE(G,G)"NO"
ENDIF
END
7、编程判断两位整数M是否为守形数。
所谓守形数是指该数本身等于自身平方的低位数,如25是守形数,因为252=625,而625的低两位为25。
是输出YES,不是输出NO。
INTEGERM
READ(G,G)M
IF(MOD(MGM,100)==M)THEN
WRITE(G,G)"YES"
ELSE
WRITE(G,G)"NO"
ENDIF
END
9已知
从键盘输入G,求出y的值。
REALG,Y
READ(G,G)G
IF(G>-15.AND.G<0)THEN
Y=COS(G+1)
ELSEIF(G>=0.AND.G<10)THEN
Y=LOG(GGG+1)
ELSEIF(G>15.AND.G<20)THEN
Y=GGG(1.0/3)
ELSE
Y=GGG
ENDIF
WRITE(G,G)"Y=",Y
END
第四章循环结构程序设计第64页1、2、5、11、12
1.写出下列程序的执行结果
(1)1.0
(2)811
(3)k=6
(4)145
(5)5167
P.652.利用下式计算π的近似值。
IMPLICITNONE
REALPI
INTEGERI
PI=0
DOI=1,N
PI=PI+1.0/(4GI-3)-1.0/(4GI-1)
ENDDO
PRINTG,PIG4
END
P66习题5
IMPLICITNONE
INTEGERN,S,M,I
LOGICALFLAG
DON=3,999,2
FLAG=.TRUE.
M=SQRT(REAL(N))
DOI=2,M
IF(MOD(N,I)==0)THEN
FLAG=.FALSE.
EGIT
ENDIF
ENDDO
IF(FLAG)THEN
S=N/100+MOD(N/10,10)+MOD(N,10)
IF(MOD(S,2)/=0)PRINTG,N
ENDIF
ENDDO
END
P.66习题11
REALY0,Y1
READ(G,G)G
Y0=G
Y1=2GY0/3+G/(3GY0GG2)
DOWHILE(ABS(Y1-Y0)>1E-5)
Y0=Y1
Y1=2GY0/3+G/(3GY0GG2)
ENDDO
PRINTG,Y1
END
P.66习题12
Reala,b,m,G
f(G)=GGG3-GGG2-1
ReadG,a,b
Dowhile(Abs(a-b)>1e-6)
m=(a+b)/2
If(f(a)Gf(m)>0)then
a=m
Else
b=m
Endif
Enddo
PrintG,(a+b)/2
END
第五章FORTRAN90数据类型第88页4、10
4.有一个三角形,顶点为G,Y,Z,其坐标分别为(1.5,2.0),(4.5,4.5),(18.0,10.5)。
求三角形的面积和重心(提示:
重心坐标=(G+Y+Z)/3利用复数计算)。
COMPLEGG,Y,Z,W
REALA,B,C,S,P
READ(G,G)G,Y,Z
A=ABS(G-Y)
B=ABS(Y-Z)
C=ABS(Z-G)
P=(A+B+C)/2
S=SQRT(PG(P-A)G(P-B)G(P-C))
W=(G+Y+Z)/3
WRITE(G,G)S
WRITE(G,G)W
END
10.设有4个圆,圆心分别在(2,2),(-2,-2),(2,-2),(-2,2),半径为1。
当G,Y坐标点落在圆上和圆内时,Z=1;落在圆外时,Z=1;落在圆外时,Z=0。
编出相应的的程序(要求利用逻辑变量)。
IMPLICITNONE
LOGICALP,Q,R,S
REALG,Y
INTEGERZ
READG,G,Y
Z=0
P=(G-2)GG2+(Y-2)GG2<=1
Q=(G+2)GG2+(Y+2)GG2<=1
R=(G-2)GG2+(Y+2)GG2<=1
S=(G+2)GG2+(Y-2)GG2<=1
IF(P.OR.Q.OR.R.OR.S)Z=1
PRINTG,Z
END
第七章数组第124页1、3、4、7、9、12
1.写出下列程序的执行结果
(1)
5.0
5.0
5.0
5.0
4.0
(2)K=5
(3)
123
456
789
101112
14710
25811
36912
(4)S=15.00
3.将一个数组的元素按逆序重新存放,例如,原来存放顺序为:
8,6,5,4,1。
要求改为:
1,4,5,6,8。
IMPLICITNONE
INTEGER,PARAMETER:
:
N=5
INTEGERA(N),I,T
DATAA/8,6,5,4,1/
WRITE(G,10)(A(I),I=1,5)
DOI=1,N/2
T=A(I)
A(I)=A(N-I+1)
A(N-I+1)=T
ENDDO
WRITE(G,10)(A(I),I=1,5)
10FORMAT(1G,5I4)
END
4.求数组中的最大和最小元素,并将最大和最小元素分别与数组的第一和第二元素互换。
IMPLICITNONE
INTEGERG(10),I,T,MA
!
变量MA存放最大数下标序号
G=(/10,55,25,70,45,15,25,85,45,35/)
WRITE(G,10)(G(I),I=1,10)
MA=1
DOI=2,10
IF(G(I)>G(MA))MA=I
ENDDO
T=G
(1)
G
(1)=G(MA)
G(MA)=T
WRITE(G,10)(G(I),I=1,10)
10FORMAT(1G,10I4)
END
7.求3×4矩阵中的最大元素和最小元素。
REALA(3,4),AMAG,AMIN
INTEGERI,J
DOI=1,3
READ(G,G)(A(I,J),J=1,4)
ENDDO
AMAG=A(1,1)
AMIN=A(1,1)
DOI=1,3
DOJ=1,4
IF(A(I,J)>AMAG)AMAG=A(I,J)
IF(A(I,J)ENDDO
ENDDO
WRITE(G,G)"最大值:
",AMAG
WRITE(G,G)"最小值:
",AMIN
END
9.输入5×5矩阵a,完成下列要求:
(1)输出矩阵a。
(2)将第2行和第5行元素对调后,输出新的矩阵a。
(3)用对角线上的各元素分别去除各元素所在行,输出新的矩阵a。
PARAMETER(N=5)
INTEGERA(N,N),I,J,T
DOI=1,N
READ(G,G)(A(I,J),J=1,5)
ENDDO
DOI=1,N
WRITE(G,10)(A(I,J),J=1,N)
ENDDO
DOI=1,N
T=A(2,I)
A(2,I)=A(5,J)
A(5,J)=T
ENDDO
DOI=1,N
WRITE(G,10)(A(I,J),J=1,N)
ENDDO
DOI=1,N
DOJ=1,N
A(I,J)=A(I,J)/A(I,I)
ENDDO
ENDDO
DOI=1,N
WRITE(G,10)(A(I,J),J=1,N)
ENDDO
10FORMAT(1G,5I5)
END
12.按以下格式打印出杨辉三角形的前N行。
1
11
121
1331
14641
15101051
PARAMETER(N=10)
INTEGERA(N,N),I,J
DOI=1,N
A(I,1)=1
A(I,I)=1
ENDDO
DOI=3,N
DOJ=2,I-1
A(I,J)=A(I-1,J)+A(I-1,J-1)
ENDDO
ENDDO
DOI=1,N
WRITE(G,10)(A(I,J),J=1,I)
ENDDO
10FORMAT(1G,10I5)
END
第八章子程序第173页1、6、8、12、13
1.指出下列错误的语句函数定义:
(1)F(G,Y)=G+Y+2+A-B正确
(2)SUM(G
(2),Y,Z)=3GY+(G
(2)+Z)GZ错误
(3)F(G,Y,G)=GGG2-Y+G+C-SIN(A)错误
(4)G2(Z,Y(I)=EGP(Z+1)-AGY(I)错误
(5)LAN(A,B,C)=AGG+(B-C)GG2-GGY正确
6.设计一个计算n!
的函数子程序,并调用该子程序计算数e的近似值。
当n!
>1E8时停止计算。
计算公式是:
FUNCTIONJC(N)
JC=1
DOI=1,N
JC=JCGI
ENDDO
END
PROGRAMEG6
REALG8E,T
E=0
T=1.0
I=0
DOWHILE(T>1E-8)
E=E+T
I=I+1
T=1.0/JC(I)
ENDDO
PRINTG,E
END
8.对任意自然数n,设计一个求n的各数位立方和的函数子程序F(n)(如F(121)=1GG3+2GG3+1GG2=1+8+1=10),并调用F(n)求1000以内的所有水仙花数(水仙花数的概念请参看例8.5)。
!
子程序
FUNCTIONF(N)
INTEGERF,N,S,T
S=0
T=N
DOWHILE(T>0)
K=MOD(T,10)
S=S+KGG3
T=T/10
ENDDO
F=S
END
!
主程序
PROGRAMEG8
INTEGERF,I
DOI=100,999
IF(F(I)==I)PRINTG,I
ENDDO
ENDDO
12.设计一个子例行程序,计算一个任意数组A中正数、负数和0的个数及其和。
!
子程序
SUBROUTINETJ(A,M,P,N,Z,S)
INTEGERA(M),M,P,N,Z,S,I
P=0
N=0
Z=0
S=0
DOI=1,M
IF(A(I)>0)THEN
P=P+1
ELSEIF(A(I)<0)THEN
N=N+1
ELSE
Z=Z+1
ENDIF
S=S+A(I)
ENDDO
END
!
主程序
INTEGERG(10),P,N,S,Z,I
DOI=1,10
READ(G,G)G(I)
ENDDO
CALLTJ(G,10,P,N,Z,S)
WRITE(G,G)"正数个数",P
WRITE(G,G)"负数个数",N
WRITE(G,G)"零数个数",Z
WRITE(G,G)"所有数的和",S
END
13.设计一个子例行程序SORT(A,N,K),其中A是一个一维数组,N是A的元素个数,SORT的功能是:
当K=1时,将数组A按升序排列;当K=0时,将数组A按降序排列;当K为其它数值时,数组A保持原序。
!
子程序
SUBROUTINESORT(A,N,K)
INTEGERA(N),K,P,T
DOI=1,N-1
P=I
DOJ=I+1,N
IF(K==1)THEN
IF(A(J)ELSE
IF(A(J)>A(P))P=J
ENDIF
ENDDO
T=A(I)
A(I)=A(P)
A(P)=T
ENDDO
END
!
主程序
PARAMETER(N=10)
INTEGERG(10),I
READ(G,G)G
WRITE(G,10)(G(I),I=1,N)
CALLSORT(G,N,1)!
将数组A按升序排列
WRITE(G,10)(G(I),I=1,N)
CALLSORT(G,N,0)!
将数组A按降序排列
WRITE(G,10)(G(I),I=1,N)
10FORMAT(1G,10I5)
END