Question

为了对 Fortran 程序有一个初步了解，下面先介绍几个简单的 Fortran 源程序。

a) 基本语句

[例 1.1] 输入两个数，求算数平均和几何平均值。[ e_121_01.f ][ e_121_01.f90 ]

[计算例]

1.0 2.0 ←键盘输入(a,b) 值

1.500000 1.414214 ←计算结果输出至屏幕

程序説明：程序中第 1 行是注释行，对程序起说明作用。F77 注释行是以“C”或“*”作为该行第一个字符的，F90 可在任意一行末以“!”开始作为注释符。第 2 行是主程序名，第 3 行是变量类型定义，第 4 行是输入语句，第 5－6 行是赋值部分，第 7 行是打印输出语句，最后是程序结束。注意 F77 的固定书写格式和 F90 的自由格式。F90 中用“；”将两行并为一行。

PROGRAM 语句： 宣布程序开始，其后跟程序名。可省略。

REAL 语句： 定义 a,b,...等变量为实数型数据。如为整数型，则用 INTEGER 定义。

变量名： 可用 a,b 等无具体意义的文字，或用 average_value 等英文缩写。使用的字符和文字长度有具体规定。

输入输出语句： 可以用以下任一种默认格式，注意星号(*) 和逗号(,)。星号意指默认的输入输出硬件（键盘和屏幕) 及格式。

READ *, 变量名

READ(*,*) 变量名

PRINT *, 变量或常数，关系式

WRITE(*,*) 变量或常数，关系式

赋值语句： 将等式右边的变量取值赋予左边的变量。

变量 = 变量、常数、关系式

例：sum = sum + x

算术运算：

２项运算：加(+)，减(-)，乘(*)，除(/)，乘方(**)。

运算的优先顺序：加，减 ＜ 乘，除 ＜ 乘方，括号中优先(　)。

单项运算：（例　-a）

数据类型： 按定义有整数、实数、复数、双精度数等。

a) 输出字符

[例 1.2] 输入圆錐底面半径 r 和高 h，求体积和表面积。[ e_121_02.f90 ]

[计算例]

Input radius r and height h ? ← 提示待输入数据的物理含义

3.0 5.2

Volume = 49.00885 ← 打印计算值

Area = 84.85442

数据的输入：READ * 语句执行时进入等待数据输入的状态。数个数据输入时以英文逗号、空格或换行符作区别，单个数据中间不能有空格。

字符的输出：在执行 READ 语句时，计算机已经进入等待数据输入的状态，但它不会给出任何提示。除了程序员以外，谁也不会知道需要输入什么数据，即使是程序员自己可能也会忘记。另外，PRINT 语句打印的数据到到底是什么物理含义也需指明。因此，需要输出文字内容。

PRINT *, 字符常量

WRITE(*,*) 字符常量

字符常量： 用'　'或"　"括起来的文字字符。

[例 1.3] 已知放射性元素的半衰期，求给定时间后的衰减量。[ e_121_03.f90 ]

[例 1.4] 给定一整数，求其自平方至５次乗方的各次乘方。[ e_121_04.f90 ]

[例 1.5] 函数的计算。[ e_121_05.f90 ]

[例 1.6] 分别求半径 R=1，3，12.5 时的圆周长。此源程序由二部分组成：主程序和子程序。[ e_121_06.f90 ]