工程分析程序设计 上机作业（一）
Fortran 程序设计基础
上机目的：掌握程序书写、字符集和标识符、数据类型、声明、算术表达式、表控输入输出语句等编程基本概念。
1、 若有实型数A=1.0 B=3.5 T=10.0 X=5.0 整型数I=-5 J=7 K=3，求出下面表达式的值。 -(A+T) (B+(X/T))/(4.0*A) (I*J)/K (I/K)*J+T/X
-(K+1)/5+I*A-B SQRT(REAL(ABS(K)+1)) MAX(J,MOD(J,K)) J+INT(T/B)/2
2、 从键盘输入三个角度值a, b, c ,计算如下式子的值：
22
34a b a b +- c b a b a tan |)||cos(||
|sin +++ 26ln()1403b c a ++
3、 已知一圆柱底面半径R 为6.4，高为5，编程求圆柱底面圆周长、圆面积，圆柱的体积和表面积；以
R=6.4为半径的圆球表面积和体积、（要求在变量声明的同时赋初值）
4、 编写一个程序，要求完成以下要求：
(1) 提示用户输入任意的3个小数；
(2）显示这三个小数；
(3）将这三个小数相加，并显示其结果；
(4）将结果按四舍五入方法转换成整数并显示。
5、 有一个六边形，求其面积。为求面积，作了3条辅助线。如图所示：（提示，三角形面积
2
a b c s ++=，a 、b 、c 为三个边长）。
=21
6、 我国现在人口13.0亿人，假定每年增长率为1.5%，求n 年以后的人口数。n 从键盘输入。
7、 请写一个程序，可以让使用者输入初始速度V 0（单位为m/s ）与角度θ（单位°），利用下列公式计算
飞行时间与距离：
g V t /sin 20θ=
t V D ?=θcos 0
《高级语言程序设计(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) 程序中使用的数组必须进行说明
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）算术运算符*
工程分析程序设计 上机作业（三） 模块化编程（2） 上机目的：进一步掌握内部例程、外部例程、接口块、模块等功能的使用方法。联系例程重载、例程递归的使用方法。 1、 编写一个子例程子程序，SUM(S ，T ，N1，N2)。把整型数N1到N2进行求和，并把求和的结果放置到S ，把整型数N1到N2进行求积，并把求积的结果放置到T 。并用这个子程序来计算 (1234)(345678)(3*4*5*6)(1*2*3)y =++++++++++- 2、 编写一个模块程序，提供以下服务：定义出常量π、e 。定义出子程序，实现求和21n i n =∑、 求阶乘n !。并在主程序中计算如下结果：从键盘上输入整数n 、实型数A 、R 、R 0，求21! n i n n =?? ???∑（实型）和0202n n R R An R e R R π??- ????? ???（实型）。 3、 编写函数子程序GDC 求两个数的最大公约数。求最大公约数的算法如下：把两个数中大的那个数作为被除数，两数相除得到一个余数。把余数去除除数得到新一轮的余数。不断重复这一过程直到余数为0，这时的除数就是两个数的最大公约数。 调用此函数，求1260，198，72三个数的最大公约数。 选作： 编写一个函数子程序求sinh(x)，以下分三个步骤完成：（1）用递归算法，求出1!1!n n x x x n n n -=?-；（2）自行编程求出234012!3!4!!n x n x x x x e x n ∞=≈+++++???=∑，要求计算精度是610!n x n -<[不能用Fortran 的内部函数EXP(x)]；（3）计算sinh()2x x e e x --=
FORTRAN笔记 2014.10.07 目录 第七讲_FORTRAN的基本知识.ppt (2) FORTRAN语言程序设计初步 (2) FORTRAN源程序的书写格式（以77为例） (2) 变量 (2) 变量类型 (2) 算术运算符和运算优先级 (3) 赋值语句 (3) 参数语句（PARAMETER语句） (3) END语句 (3) PAUSE语句 (3) 逻辑运算和选择结构 (4) 关系表达式 (4) FORTRAN中数组的定义及使用 (4) 其他 (5) 1. fortran语言定义CHARACTER*6 TTL(14,3),CNAM(400)是什么意思? (5) 2. fortran里character*10 是什么意思 (5) 3. Fortran中kind是什么函数？ (5)
第七讲_FORTRAN的基本知识.ppt FORTRAN语言程序设计初步 FORTRAN是Formula Translation的缩写，意为“公式翻译”，它是为科学、工程问题或企事业管理中的那些能够用 数学公式表达的问题而设计的，其数值计算的功能较强。 常用的是FORTRAN77和FORTRAN90两种标准。 1、一个程序由若干个程序单位组成。主程序和每一个子程序分别是一个独立的程序单位。 2、每一个程序单位都是以“END”结束的。 3、一个程序单位包括若干行。 1）语句行。由一个FORTRAN语句组成。 2）非语句行，即注释行。 4、FORTRAN程序中的语句可以没有标号，也可以有标号，根据需要而定。标号的作用是标志一个语句以便被其 他语句引用。 5、一个程序单位中各类语句的位置是有一定规定的。 6、FORTRAN源程序必须按一定的格式书写。 FORTRAN源程序的书写格式（以77为例） 每一行有80列，分别如下： 1、第1-5列为标号区。一行中第一列为“C”或“*”，该行即被认为是注释行。 2、第6列为“续行标志区”，如果在一行的第6列上写一个非空格和非零的字符，则该行作为其上一行的续行。 3、第7-72列为语句区。 4、第73-80列，注释区。 变量 变量名：一个变量需要用一个名字（变量名）来识别。在同一个程序单位中不能用同一个变量名代表两个不同的变 量。 FORTRAN的变量名按以下规则选定： 1）第一个字符必须是字母，即变量名必须以字母开头； 2）在一个字母后面可以跟1-5为数字或字母。 如果选定的变量名超过6个字符，则只有前面6个字符有效。 注：在变量名中大写与小写字母是等价的。 变量类型 整型变量Integer、实型变量Real、双精度变量Double Precision、复型变量Complex、逻辑型变量Logical和字符型变量Character。 1、隐含约定（I-N规则） FORTRAN规定：在程序中的变量名，凡以字母I，J，K，L，M，N六个字母开头的，即认为该变量为整型变量。 在程序中，凡是变量名以字母I,J,K,L,M,N,i,j,k,l,m,n开头的变量被默认为整型变量，以其他字母开头的变量被 默认为实型变量。 2、用类型说明语句确定变量类型 1）INTEGER语句（整型说明语句） 2）REAL语句（实型说明语句） 3）DOUBLE PRECISION语句（双精度说明语句） 4）COMPLEX语句（复型说明语句） 5）LOGICAL语句（逻辑型说明语句）
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
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 复习题汇总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
第四章 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.p rogram 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 implicit none integer money real tax write(*,*) "请输入月收入" read(*,*) money if ( money<1000 ) then tax = 0.03 else if ( money<5000) then tax = 0.1 else tax = 0.15 end if write(*,"(' 税金为'I8)") nint(money*tax) end program 2.program main implicit none integer day character(len=20) :: tv write(*,*) "请输入星期几" read(*,*) day select case(day) case(1,4) tv = "新闻" case(2,5) tv = "电视剧" case(3,6) tv = "卡通" case(7) tv = "电影" case default write(*,*) "错误的输入" stop end select write(*,*) tv end program 3.program main implicit none integer age, money real tax write(*,*) "请输入年龄" read(*,*) age write(*,*) "请输入月收入" read(*,*) money if ( age<50 ) then if ( money<1000 ) then tax = 0.03 else if ( money<5000 )then tax = 0.10 else tax = 0.15 end if else if ( money<1000 ) then tax = 0.5 else if ( money<5000 )then tax = 0.7 else tax = 0.10 end if end if write(*,"(' 税金为'I8)") nint(money*tax) end program 4.program main implicit none integer year, days logical mod_4, mod_100, mod_400 write(*,*) "请输入年份" read(*,*) year mod_4 = ( MOD(year,4) == 0 ) mod_100 = ( MOD(year,100) == 0 ) mod_400 = ( MOD(year,400) == 0 ) if ( (mod_4 .NEQV. mod_100) .or. mod_400 ) then days = 366 else days = 365 end if write(*,"('这一年有'I3'天')") days stop end program 第六章 1.program main implicit none integer i do i=1,5 write(*,*) "Fortran" end do stop end program 2.program main implicit none integer i,sum sum = 0 do i=1,99,2 sum = sum+i end do write(*,*) sum stop end program 3.program main implicit none integer, parameter :: answer = 45 integer, parameter :: max = 5 integer weight, i do i=1,max write(*,*) "请输入体重" read(*,*) weight if ( weight==answer ) exit end do if ( i<=max ) then write(*,*) "猜对了" else write(*,*) "猜错了" end if stop end program 4.program main implicit none integer, parameter :: max=10 integer i real item real ans ans = 1.0 item = 1.0 do i=2,max item = item/real(i) ans = ans+item
贵州大学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) 十六进制
工程分析程序设计上机作业（五） 模块化编程 上机目的：掌握内部例程、主程序、外部例程、模块等功能的使用方法。 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语言) (总分：36.00，做题时间：90分钟) 一、 (总题数：36，分数：36.00) 1.编译程序能将高级语言编写的源程序转换成( )。 A．解释程序 B．汇编程序 C．映象程序 D．目标程序 （分数：1.00） A. B. C. D. √ 解析： 2.一个完整的FORTRAN源程序( )。 A．至少包括一个主程序 B．至少包括一个主程序和一个子程序 C．由一个主程序与一个以上的子程序组成 D．由一个主程序与一个子程序组成 （分数：1.00） A. √ B. C. D. 解析： 3.语句函数定义语句在程序内合法的位置是( )。 A．在程序块开头语句之后，END语句之前 B．在程序块开头语句之后，可执行语句之前 C．在说明语句之后，END语句之前 D．在说明语句之后，可执行语句之前 （分数：1.00） A. B. C. D. √ 解析： 4.下列关于函数子程序虚实参数的错误说法是( )。 A．可以没有形参数 B．虚实结合的数组长度可以不同 C．实参表与虚参表类型可以不同 D．函数名可以作为虚参
（分数：1.00） A. B. C. √ D. 解析： 5.下列叙述中正确的是( )。 A．FORTRAN程序块中，无名公用语句只能有一个B．FORTRAN子程序中，至少应有一个RETURN语句C．FORTRAN程序块中，最后一行必须是END语句D．FORTRAN程序块中，必须有变量说明语句 （分数：1.00） A. B. C. √ D. 解析： 6.运行下面的程序时得不到所需的结果，其主要原因是( )。INTEGER X(11) DATA X/9，8，7，6，5，4，3，2，1，0，-1/ DO 10 1=1，X(1) ，-1 WRITE(*，*)1.0/SQRT(25.0-REAL(X(1) )* * 2) 10 CONTINUE END A．没有给X(11) 赋初值 B．发生除以零的情况 C．发生负数开平方的情况 D．循环参数设置错误 （分数：1.00） A. B. C. D. √ 解析： 7.下列数据中，不符合FORTRAN常量表示法的是( )。 A．-25.6 B．2.758D3 C．'FOOT"=' D．TRUE （分数：1.00） A. B. C. D. √ 解析：
《复合材料力学》课程上机指导书（力学101，力学C101-2） 河北工业大学机械学院力学系 2013年9月
目录 作业1单向板刚度及柔度的计算(2学时) (1) 作业2单向板的应力、应变计算(2学时) (2) 作业3绘制表观工程常数随 的变化规律(3学时) (3) 作业4绘制强度准则的理论曲线(包络线)(3学时) (4) 作业5层合板的刚度计算(3学时) (5) *作业6层合板的强度计算(4学时) (6) 附录作业提交说明……………………………………………. . 7 注：带“*”的题目可作为自愿选做题。
作业1 单向板刚度及柔度的计算 一、要 求 （1）选用FORTRAN 、VB 、MAPLE 或MATLAB 编程计算下列各题； （2）上机报告内容：源程序、题目内容及计算结果； （3）材料工程常数的数值参考教材自己选择； （4）上机学时：2学时。 二、题 目 1、已知单层板材料工程常数1E ，2E ，12G ，计算柔度矩阵[S ]和刚度矩阵[Q ]。（玻璃/环氧树脂单层板材料的MPa 1090.341?=E ，MPa 1030.142?=E ， MPa 1042.0412?=G ，25.021=μ， MPa 1001=σ，MPa 302-=σ，MPa 1012=τ） 2、已知单层板材料工程常数1E ，2E ，12G ，21μ及θ，计算柔度矩阵][S 和刚度矩阵][Q 。（M P a 1090.341?=E ，MPa 1030.142?=E ，MPa 1042.0412?=G ，25.021=μ，?=30θ）
作业2 单向板的应力、应变计算 一、要 求 1、选用FORTRAN 、VB 、MAPLE 或MATLAB 编程计算下列各题； 2、上机报告内容：源程序、题目内容及计算结果； 3、材料工程常数的数值参考教材自己选择； 4、上机学时：2学时。 二、题 目 1、已知单向板的应力x σ、y σ、xy τ，工程常数1E ，2E ，12G ，21μ及θ，求x ε、 y ε、xy γ；1σ、2σ、12τ；1ε、2ε、12γ。 （知?=30θ，应力MPa 160=x σ，MPa 60=y σ，MPa 20=xy τ，工程常数MPa 1090.341?=E ，MPa 1030.142?=E ，MPa 1042.0412?=G ，25.021=μ，?=30θ） 2、已知1σ、2σ、12τ，工程常数1E ，2E ，12G ，21μ及θ，求1ε、2ε、12γ；x ε、y ε、 xy γ；x σ、y σ、xy τ。 （知MPa 1001=σ，MPa 302-=σ，MPa 1012=τ，MPa 1090.341?=E ，MPa 1030.142?=E ，MPa 1042.0412?=G ，25.021=μ，?=30θ）
一、选择题 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
FORTRAN第二次实习作业 2014301610307 -D-马莉1 、根据最新个人所得税计算方法，如下面的个税税率表。工资、薪金所得，适用七级超额累进税率，税率为百分之三（3% ）至百分之四十五（45% ）。3500 元为起征点，即只需要缴纳高于3500 元部分的个税。请根据下表编程，通过给定个人收入，计算其所需邀纳的个税。 program jiaoshui implicit none integer m real n character i 100 print*,'应缴纳所得税是否含税，请回答Y或N：' read*,i print*,'请输入工资所得：' read *,m if(i=='Y')then select case(m) case(0:1500) n=m*3/100-0 case(1501:4500) n=m*10/100-105 case(4501:9000) n=m*20/100-555 case(9001:35000) n=m*25/100-1005 case(35001:55000) n=m*30/100-2775 case(55001:80000) n=m*35/100-5505 case default n=m*45/100-13505 end select elseif(i==’N’)then select case(m)
case(0:1455) n=m*3/100-0 case(1456:4155) n=m*10/100-105 case(4156:7755) n=m*20/100-555 case(7756:27255) n=m*25/100-1005 case(27256:41255) n=m*30/100-2775 case(41256:57505) n=m*35/100-5505 case default n=m*45/100-13505 end select else print*,"输入格式不正确，请重新输入：" goto 100 endif print*,n End 2、编写程序求某天在当年中是第多少天，例如2010 年1 月10 号是2010 年的第10 天。要求分别使用IF 、CASE结构来编写。program tianshu2 implicit none integer M1(12),i,y,m,d,sum,n data M1/31,28,31,30,31,30,31,31,30,31,30,31/ sum=0 100 print *,"请输入日期格式例如：2010 10 23" read 10,y,m,d 10 format(i4,1x,i2,1x,i2) if(m>12.or.d>31)then
工程分析程序设计上机作业（五） 数组（1） 上机目的：练习数组的声明、存储、操作，以及数组参数、动态数组、数组函数的使用。 1、打印杨辉三角形（格式不限）。 1 1 1 1 2 1 1 3 3 1 1 4 6 4 1 program yanghui_triangle implicit none real::T(5,5) !第一个数代表列，第二个代表行 integer i,j Do j=1,5 T(1,j)=1 T(j,j)=1 end Do DO j=3,5 DO i=2,j T(i,j)=T(i-1,j-1)+T(i,j-1) end Do end Do print*,T end program 2、输入两个矩阵，并用矩阵作为子程序的参数，用子程序完成任意两个矩阵的乘法。（如 果可能，用数组函数子程序来完成这一功能） program array integer,parameter::m=2,n=2 !只需要改变常量的大小即可改变数组形状 real A(m,n),B(m,n),C(m,n) read*,A,B call ar(A,B,m,n,C) print*,C end program
subroutine ar(A,B,m,n,C) integer m,n real A(m,n),B(m,n),C(m,n) C=A*B end subroutine 3、用“冒泡算法”对一个数列A(I)进行从小到大排序，步骤如下：（1）若A(2)
中南林业科技大学课程考试试卷 课程名称： 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
相关主题
文本预览