第九章函数与子例行程序.pdfVIP

第九章 函数与子例行程序 9.1 概述 9.2 语句函数 9.3 函数子程序 9.4 子例行程序 9.5 形参与实参的结合 9.6 外部子程序 9.7 递归子程序 9.1 概述 Fortran 程序单元大体分为主程序、子程序两种。当程序中多次 出现相同或相似的运算时，或者需要将复杂问题分解为若干个 小问题分别处理时，可以用子程序编程。 两种 Fortran 子程序： (1) 函数子程序 包括标准函数 ( 例如 sqrt, exp, int) 和用户自定义函数。 (2) 子例行程序 包括标准子例行程序 ( 例如 random(x), date(string) ) 和自定义子 例行程序。 9.2 语句函数 9.2.1 语句函数的定义 一般形式：语句函数名 (x1,x2,…,xn) ＝表达式 例如： REAL d, k, x, y, A, T d(x,y)=sqrt(x*x + y*y) k(T)=A*exp(-Ea/R/T) (1) 语句函数的类型可用类型声明语句定义或遵守 I-N 规则。 语句函数名不能与变量重名。 (2) 自变量x1,x2,…,xn 只能是普通变量，表达式必须引用全部 自变量。 (3) 语句函数的定义是非执行语句，只能写程序说明部分。 9.2.2 语句函数的调用 语句函数的使用方法与标准函数相同。 例 : 计算空间一点 P(x,y,z) 的方向余弦 x cosα = 2 2 2  x y z y cosβ = 2 2 2  x y z z cosγ = 2 2 2  x y z α 、 β 、γ 为点 P 到原点的直线与坐标轴 x 、y 、z 的夹角。 REAL F,DA,DB,DG,a,b,c,r,x,y,z F(a,b,c,r) = r/SQRT(a*a+b*b+c*c) READ(*,*) x,y,z DA = F(x,y,z,x) DB = F(x,y,z,y) DG = F(x,y,z,z) WRITE(*,*) DA,DB,DG END 调用时用实在参数（有实际值的常量、变量、表达式）代替 定义时的形式参数（自变量 x1,x2,…,xn ）。实参的个数和类 型必须与形参完全一致。 9.3 函数子程序 语句函数只能定义在一个表达式内。对于复杂的函数，必须 使用函数子程序。 函数子程序的定义： 类型说明 FUNCTION 函数名 ( 形式参数表 ) 说明语句 可执行语句 END FUNCTION 例：求 n! 的函数子程序 REAL FUNCTION fac(n) INTEGER i,n; REAL :: f=1.0 DO i=1,n f=f*i ENDDO fac=f END FUNCTION (1) 函数体的说明部分对形式参数和内部变量说明类型。内部 变量只在本函数内有效，与其他程序单元中的同名变量无关。 (2) 函数名在函数体中至少要赋值一次：函数名= 表达式，因 为函数名负责把函数值带回调用程序。 (3) 函数体中可设置零条、一条或多条RETURN 语句，执行到 此语句时返回调用程序。 (4) 调用时实在参数与形式参数的个数和类型必须完全相同。 例：求 1!+3!+5!+…+9! INTEGER i; REAL :: sum=0 DO i=1,9 sum=sum+fac(i) ENDDO PRINT *, sum CONTAINS REAL FUNCTION fac(n)