华中科技大学 FORTRAN课件 第11章 极力推荐入门教程!!!!!

华中科技大学土木工程与力学学院 《Fortran语言》
9
11.1 函数子程序
11.1.1 函数子程序的定义:
函数体中的执行语句完成求解函数值的运算.在 程序的执行流程返回调用单位之前,必须把所 求得的函数值赋给函数名。 在子程序中必定有如下语句存在: 函数名=表达式
如: MYABS= -1*X 和 MYABS=X 而不能: MYABS(X)= -1* X 和 MYABS(X)=X
华中科技大学土木工程与力学学院 《Fortran语言》
13
11.1 函数子程序
11.1.2 函数子程序的调用:
函数子程序的调用与系统函数（内部函数）的 调用形式完全相同。
变量名= 含有（函数名（实参） 或 函数名（））的表达式
注意以下四点：
1）只能在表达式中对函数进行调用； 2）对于主程序要调用外部函数来说，应在说 明部分对函数名的数据类型加以说明。如果没 有说明，则按隐含类型处理。必须在调用程序 中使外部函数名与函数子程序定义中函数名的 类型完全相同。
e e sinh( x) 2
x
x
10
x x x e 1 x ... 2! 3! n 0 n!
华中科技大学土木工程与力学学院 《Fortran语言》
3

n
20
华中科技大学土木工程与力学学院 《Fortran语言》
21
Program main1 Real:: MyExp,y,x Write(*,*) 'Enter x:' Read(*,*)x Y=(MyExp(x)-MyExp(-x)) /2.0 Write(*,*) Y END PROGRAM Function MyExp(x) result (Exp_res) Real::Exp_res,x,s,T integer::N N=1; T=1; S=1 do T=T*X/N IF(ABS(T)<=1E-5)EXIT S=S+T;N=N+1 end do x Exp_res=S END FUNCTION
。
华中科技大学土木工程与力学学院 《Fortran语言》
25
GO…
华中科技大学土木工程与力学学院 《Fortran语言》
26
11.2 子例行程序
Fortran77除了可以使用函数子程序 (FUNCTION)之 外,还可以使用子例行程序(SUBROUTINE).它们都属 于子程序的范围。 函数子程序与子例行程序区别：
华中科技大学土木工程与力学学院 《Fortran语言》
6
11.1 函数子程序
用子程序(程序段)来定义一个函数—称之
为“外部函数”； 通常是在内部函数没有这个功能(未提 供),而语句函数又不能完成所需的功能时, 采用函数子程序; 函数子程序是一个单独的程序单位—它 有明显标志（语句函数属于所在的程序 单位，只能被本程序单位中的表达式引 用），被其它程序调用。
。。。。
。。。。 。。。。
end
。。。。
。。。。 。。。。
。。。。
end
华中科技大学土木工程与力学学院 《Fortran语言》
5
本章要点
11.1 11.2 11.3 11.4 11.5 11.6 函数子程序; 子例行程序; 实参与虚参之间的数据传递(输) 利用子程序实现程序的模块化设计; 子程序中的SAVE语句和DATA语句; 程序举例
华中科技大学土木工程与力学学院 《Fortran语言》
17
PROGRAM MAIN N=3 Y=(S(N)+S(N+1)+S(N+2))/(S(N+3)+S(N+4)) WRITE(*,100) ‘Y=‘,Y 100 FORMAT(1X,A,E13.6) END FUNCTION S(X) INTEGER X (1 2 3) (1 2 3 4) (1 2 3 4 5) S=0 Y (1 2 3 4 5 6) (1 2 3 4 5 6 7) DO 10 I=1,X S=S+I 10 CONTINUE END
27
11.2 子例行程序
11.2.1 子例行程序定义
调用方式不同：函数子程序的调用总出现 在表达式中，子例行程序的调用是由一条 独立的语句CALL来完成的。 使用场合：如果只需要得到一个计算值， 采用函数子程序。在设计函数子程序时， 不主张在函数体内改变形参的值，从而保 证不使对应的实参发生变化。如果想通过 程序运算后得到一批数据（矩阵运算）或 一些特定的操作，就使用子例行程序。
第十一章 子 程 序
华中科技大学土木工程与力学学院 《Fortran语言》
1
概
念
一个Fortran应用程序往往是由一个主程 序(必要条件)和多个子程序组成，对于 一个完整的应用系统几乎都是这样。 Fortran子程序分为: 1)函数子程序; 2)子例行程序; 3)数据块子程序(后章介绍)
华中科技大学土木工程与力学学院 《Fortran语言》
计算y=f（x）+g（x）值。 f（x）=x+sin30° x+1，f(x)>0 g （ x ）= 0， f(x)=0 x-1， f(x)<0 要求用外部函数求g（x）的值
FUNCTION GF(X,F) REAL X, F, GF IF (F.GT.0) THEN GF=X+1 ELSEIF (F.LT.0) THEN GF=X-1 ELSE GF=0 ENDIF END
华中科技大学土木工程与力学学院 《Fortran语言》
14
11.1 函数子程序
11.1.2 函数子程序的调用:
如上例中：由于外部函数MYABS是实型， 主程序中就应该说明了MYABS为实型， 不然就造成类型不一致的错误。 3）在调用外部函数时，实参表中的实参 个数必须与对应函数的形参个数完全相 同，且实参与形参之间是建立的一一对 应关系，类型也要求一致。
华中科技大学土木工程与力学学院 《Fortran语言》
10

REAL FUNCTION MYABS(X) REAL X IF(X .LT. 0) THEN MYABS= -1*X RETURN ELSE MYABS=X ENDIF END PROGRAM MAIN REAL MYABS, X, A READ(*,*) A X=MYABS(A) WRITE(*,*) A, X END
华中科技大学土木工程与力学学院 《Fortran语言》
28
11.2 子例行程序
11.2.1 子例行程序的定义
子例行程序必须以Subroutine语句开头,以END语 句结束.其基本形式如下: Subroutine 子例行程序名 (虚参表) 其中子例行程序名的命名规则与变量名相同 .但 子例行程序名只是标识一个子例行程序，不代表 任何值. (虚参表)是子例行程序与调用单位之间进行数据 传输的主要渠道。没有虚参表时,其括号可以省 略;多个虚参时,以逗号相隔。 有关其它规则与函数子程序基本相同。
!
函数子程序--函数子程序的名字代表一个值（在 函数子程序中求出的函数值并存放在函数名中， 因此函数名有数据类型）； 子例行程序—在子例行程序中求出的值不由该程 序名带回调用单位，是通过实参与形参的联系而 传递的；所以子例行程序名不属于某个类型，也 不代表某个值，也不需要对该名进行说明。
华中科技大学土木工程与力学学院 《Fortran语言》
华中科技大学土木工程与力学学院 《Fortran语言》
18
如果没有对X进行整型说明就 会产生编译错误．
华中科技大学土木工程与力学学院 《Fortran语言》
19
PROGRAM MAIN REAL MYEXP WRITE(*,*) ‘ 输入X的值：’ READ（*，*） X Y=（MYEXP（X）-MYEXP（-X））/2.0 WRITE（*，*）Y END FUNCTION MYEXP（X） REAL MYEXP N=0 T=1 S=1 DO 10 ，WHILE（ABS（T）.GT.1E-5） N=N+1 2 T=T*X/N x S=S+T CONTINUE MYEXP=S END
华中科技大学土木工程与力学学院 《Fortran语言》
一定要注意主、子 程序中变量类型说 明和函数名的类型 说明。
11
华中科技大学土木工程与力学学院 《Fortran语言》
12
11.1 函数子程序
11.1.1 函数子程序的定义:
函数在执行中遇到Return语句或End语句的时 候,立即返回到调用单位并由函数名把函数值带 回. Return语句:只能在子程序中,可以有多个; End语句:在主、子程序只能有一个,不能省略. 注意:函数是一个独立的编译单位,除函数名和 虚参外与其它程序单位没有任何关系或联系; 因此在函数子程序中不必考虑标号、名字等是 否和其它程序单位相重(具有很大的灵活性)。
华中科技大学土木工程与力学学院 《Fortran语言》
15
பைடு நூலகம்
如果没有对 Myabs说明 的话，就会 出现错误。
华中科技大学土木工程与力学学院 《Fortran语言》
16
11.1 函数子程序
11.1.2 函数子程序的调用:
在函数调用时，实参的值传递给对应位 置上的形参，程序就接着转去执行函数 体内的执行语句。在遇到End或Return 语句返回调用单位，并把函数值带回,然 后继续参与主程序的表达式的运算。 4）如果函数没有形参，则在调用时也没 有实参，但括号不能丢。
华中科技大学土木工程与力学学院 《Fortran语言》
7
11.1 函数子程序
11.1.1 函数子程序的定义: