Fortran语言

Fortran语言

（Fortran77结构化程序设计）

谭浩强 田淑清 编著

普通高校计算机基础教学系列课程

说明

??

第一章 算法（略）

第二章 计算机和计算机程序（略） 第三章

Fortran语言程序设计初步

3.1 Fortran语言的发展概况

FORmula TRANslator 公式翻译器

Fortran语言是一种适合于科学和工程计算问题的程序设计语言；也是计算机高级语言发展史上使用最早的一种程序设计语言。 3.1 Fortran语言的发展概况 Fortran语言的发展过程：

?????????????

1951年由约翰·贝克斯等人开始研究Fortran语言； 1957年第一个Fortran程序在IBM704机上运行； 1958和1962年先后推出FortranⅡ和FortranⅣ； 1966和1978年先后颁布Fortran66和Fortran77； 1991年颁布了Fortran90。

现在已经有了Fortran 95, Fortran 2003 Fortran语言的特征 Fortran语言的最大特性是接近数学公式的自然描述，在计算机里具有很高的执行效率。 易学，语法严谨。

可以直接对矩阵和复数

!!"#$"%与&&’’"()*#*"+语言混合编程

陈建平!罗先启

!三峡大学土木水电学院"湖北宜昌!###,""!

摘要!对-./"分析了0.1与2233和456/6.7语言进行混合编程的常用方法进行了介绍#

其实现方式和各自的利弊#并用实例对-89文件实现方式进行了较详细的论述+关键词!-./"0.1$233$56/6.7$!2!4!混合编程

中图分类号##/:,&,!!!文献标识码#.!!!文章编号#&%(!$’#)9"!""#"%$";#($";

%!,-./0121455,6.78629221!"#$"%<,=8&&’’46/()*#*"+33#3:;

2<=>?@A>@>DE9@A>@!0BCF

"$$#2EGG==EH2@I@GJKMNEEO=N8>@>==N@>2<@>A/@I+S@T#,""!$2<@>ACLBCCCC#">?=147=<@QBA=N@

第三章

二、将数学式写成fortran表达式 (1)a**2+4*b**3)/(a-b)

(2)(-b+sqrt(b*b-4*a*c))/(2*a) (3)(6*sin((x+y)**2))/(2*a) (4)sin(y/(aqrt(x*x+y*y)))

(5)sin(atan(aqrt(x*x+y*y))/(abs(c))) 九、电路编程

Implicit real(a-z) R1=30 R2=60 R3=45 U=120

I=(u/r1+u/r2+u/r3) Print*,i End

结果：8.666667

十、六边形面积编程 ! main program

Implicit real(a-z) I1=10 I2=20 I3=16 I4=13 I5=21 I6=14 I7=30 I8=36 I9=28

Area=x(i2,i3,i7)+x(i1,i7,i8)+X(i4,i8,i9)+x(i5,i6

imsl7.0中用use linear_operators这句话是会出错的，当时intel论坛上也有人问，后来intel给出了个X64 imsl的补丁，但是32位的没有。

要使用包含在linear_operators这个库中的函数时，要用use+原函数 例如：上面那个例子把use linear_operators改成use operation_xt即可 另外，imsl7.0引用函数和以前版本不一样的 补充一点：

linear_operators这个文件是这样子的（一看就知道原因了）： modulelinear_operators

usecond_int usedet_int usediag_int

usediagonals_int usefft_int useifft_int useeye_int

uselin_eig_self_int uselin_sol_self_int usenorm_int useoperation_i useoperation_ix useoperation_t useoperation_h useoperation_tx useoperation_hx useoperation_x useoperation_

imsl7.0中用use linear_operators这句话是会出错的，当时intel论坛上也有人问，后来intel给出了个X64 imsl的补丁，但是32位的没有。

要使用包含在linear_operators这个库中的函数时，要用use+原函数 例如：上面那个例子把use linear_operators改成use operation_xt即可 另外，imsl7.0引用函数和以前版本不一样的 补充一点：

linear_operators这个文件是这样子的（一看就知道原因了）： modulelinear_operators

usecond_int usedet_int usediag_int

usediagonals_int usefft_int useifft_int useeye_int

uselin_eig_self_int uselin_sol_self_int usenorm_int useoperation_i useoperation_ix useoperation_t useoperation_h useoperation_tx useoperation_hx useoperation_x useoperation_

FORTRAN 90标准函数（一）

(2012-07-03 17:14:57) 转载fortran 函数 教育

▼

标签： 分类： 学习

符号约定：

? I代表整型;R代表实型;C代表复型;CH代表字符型;S代表字符串;L代表逻辑型;A代表数组;P代表指针;T代表派生类型;AT为任意类型。

? s:P表示s类型为P类型(任意kind值)。s:P(k)表示s类型为P类型(kind值=k)。 ? […]表示可选参数。 ? *表示常用函数。

表1 数值和类型转换函数 函数名 ABS(x)* AIMAG(x) AINT(x[,kind])* AMAX0(x1,x2,x3,…)* AMIN0(x1,x2,x3,…)* ANINT(x[,kind])* CEILING(x)* CMPLX(x[,y][,kind])) CONJG(x) DBLE(x)* DCMPLX(x[,y]) DFLOAT(x) DIM(x,y)* DPROD(x,y) FLOAT(x)* FLOOR(x)* IFIX(x)* IMAG(x) INT(x[,kind])* LOGICAL(x[,kind])* MAX(x1,x2,x3,…)* MAX1(x

中南大学

本科生课程设计(实践)设计报告

(程序设计基础)

题 目 学生学号 学生姓名 指导教师 学 院 专业班级

微分方程求解问题

1208100101 马晶 蔡 旭 晖 土木建筑学院 土建类10-01班

计算机基础教学实验中心

2011 年 6 月 30 日

Fortran语言程序设计

摘要：《FORTRAN程序设计》课程设计是在学生修完“FORTRAN程序设计基础”之后的一门实践性课程，要求学生完成一个大型程序设计任务。通过本课程，培养学生利用程序设计技术进行计算机与专业知识的交叉应用、综合解决实际问题的能力。让学生通过自己分析问题、设计算法、编写调试程序，以解决具体的实际问题，要求掌握结构化程序设计与调试方法，得到分析问题、解决问题的初步训练，提高综合应用能力。

关键字：微分方程求解，QuickWin绘图。

题目：意大利生物学家Ancona曾致力于鱼类种群相互制约关系的研究，他从第一次世界大战期间，地中海各港口捕获的几种鱼类捕获量百分比的资料中，发现鲨鱼的比例有明显增加。为了反映在没有人工捕获的自然环境中食饵与捕食者之间的制约关系，建立一个食饵—捕食者系统的数学模型如下：

?dx1?x1(r1??1x2)?

习题二

一、问答题

1. 给出下面变量名称,哪些是合法变量？哪些是非法变量？说明原因。

Count、num_2、x&y、4x+5y、china-suzhou、$us、AbCdE、Mr.bai、t5、 _another、school_class_25、#125、2002Y、π、β、A01/02、alpha、date(1) 1. 判定下面整数,指出哪些是合法整数,哪些是非法整数？说明原因。

-0、+ 215、$125、3,245,895、5.3245、5#384、-524_3、#5DFE、 23-345、16#1A2B、38#ABCD、8#275_2、+327890、4 #321 2. 判定下面实数,指出哪些是合法实数,哪些是非法实数？说明原因。

-0E2、45.2345E3.5、-5489E25_8、-.2345E-35、$185.45E、+ 2.753425E24_3、 58D85、+0.E-0、-00000.001E5、5,443,223.44、-12 34E+2、+ 18.5E 18、2.5E42

习题三

一、选择题

1．下面是Visual Fortran中正确的表达式是 。

（A）A*COS(X)+∣B

三、选择结构程序设计

1.开机；进入Fortran PowerStation4.0集成工作环境，了解工作环境及其菜单的使用，输入如下FORTRAN90源程序，并以A1.F90存盘。

PROGRAM EX1 IMPLICIT NONE PRINT *，'2+3=',2+3

PRINT *，'This is a program！' END PROGRAM EX1

编译、链接、运行后，结果是：2+3=5 This is a program!

2.输入并调试以下程序，查看运行结果。 PROGRAM EX2 IMPLICIT NONE INTEGER::A,B,C READ *,A,B C=A+B

PRINT *,'A+B=',C END PROGRAM EX2

运行时，输入2，3↙，则结果是：2+3=5 This is a program! 3.程序设计

2⑴ 已知一个球的半径，求它的表面积（4?R）、体积（432直径(2R)、周长（2?r）和面积（?R）。

?R3）及经过球心的横切面的

⑵求方程X2+3X-1=0的根

3.(1). (2

Fortran 库参考

Sun Studio 11

Sun Microsystems,

文件号码 819-4757-102005 年 11 月，修订版 A

请将有关本文档的意见和建议提交至：/hwdocs/feedback

版权所有 © 2005 Sun Microsystems, Inc., 4150 Network Circle, Santa Clara, California 95054, U.S.A. 保留所有权利。

美国政府权利 － 商业用途。政府用户应遵循 Sun Microsystems, Inc. 的标准许可协议，以及 FAR（Federal Acquisition Regulations，即“联邦政府采购法规”）的适用条款及其补充条款。必须依据许可证条款使用。本发行版可能包含由第三方开发的内容。

本产品的某些部分可能是从 Berkeley BSD 系统衍生出来的，并获得了加利福尼亚大学的许可。UNIX 是 X/Open Company, Ltd. 在美国和其他国家/地区独家许可的注册商标。

Sun、Sun Microsystems、Sun 徽标、Java 和 JavaHelp 是 Sun Microsystems, Inc.