VisualFortran90第六章1.pptVIP

第6章 程序单元和过程 复习：子程序编写思路 程序单位定义 虚元表 子程序所有输入和除结果名以外的输出 结果名 不能和函数名同名 说明语句 虚元表中的所有虚元、函数结果名；其它中间变量 不能说明函数名 执行语句 完成子程序单位的所有功能 至少有一次对结果名赋值 程序单位结束语句 END Function 主调程序： PROGRAM MAIN_SUM IMPLICIT NONE INTEGER::N,K,ISUM REAL::X,S READ*,N,X S=0.0 DO K=1,N S=S+X**K/ISUM(K) ENDDO PRINT*,S=,S END PROGRAM MAIN_SUM 算法2：把求和的工作全部由函数实现，这时需要给函数设计两个变元（n 和x），如下所示： FUNCTION SUM2(N,X) RESULT(SUM2_RESULT) IMPLICIT NONE INTEGER::N,I,J,T=0 REAL::X,SUM2_RESULT SUM2_RESULT=0 DO I=1,N T=T+I SUM2_RESULT=SUM2_RESULT+X**I/T ENDDO END FUNCTION SUM2 主调程序为： PROGRAM MAIN_SUM2 IMPLICIT NONE INTEGER::N REAL::SUM2,X,S READ*,N,X S=SUM2(N,X) PRINT*,S=,S END PROGRAM MAIN_SUM2 调用子程序的执行流程 主调程序： PROGRAM MAIN IMPLICIT NONE INTEGER::N,PRIME DO N=20,50 IF(PRIME(N)==1) PRINT*,N ENDDO END PROGRAM MAIN 改写PRIME函数（使用逻辑变量）： FUNCTION PRIME(M) RESULT(PRIME_RESULT) IMPLICIT NONE INTEGER::M,I LOGICAL::PRIME_RESULT PRIME_RESULT=.FALSE. DO I=2,M/2 IF(MOD(M,I)==0) EXIT ENDDO IF(IM/2) PRIME_RESULT=.TRUE. END FUNCTION PRIME 首先将该多边形划分成若干个三角形，每个三角形的边长分别用L1、L2、L3表示。 然后用下面的公式（海伦公式）计算每个三角形的面积。最后将这几个三角形的面积相加： 以上问题要经过三次完全相同的运算，但参加运算的L1、L2、L3量值不同。我们可以编写一段程序实现这一功能，并为该程序段命名。在计算每个三角形面积时调用该程序，并代入不同的参数完成运算，从而简化设计，便于结构化的程序设计。 练习题： 一、选择题 1．下面的函数过程说明语句中，错误的是 ： A）FUNCTION MYFUN(X )RESULT( ) B）FUNCTION MYFUN(X )RESULT(MYFUN_RESULT) C）FUNCTION MYFUN( )RESULT(RESULT_MYFUN) D）FUNCTION MYFUN( )RESULT(MYFUN_RESULT) 练习题： 2．下面的子例行程序过程说明语句中，正确的是 ： A）SUBROUTINE MYSUB( ) RESULT(RE) B）SUBROUTINE MYSUB C）SUBROUTINE MYSUB(1,2) D）SUBROUTINE MYSUB(AREA) RESULT(RE) FUNCTION AREA(A,B,C) RESULT(AREA_RESULT) REAL,INTENT(IN)::A,B,C REAL::P,AREA_RESULT P=(A+B+C)/2. AREA_RESULT=SQRT(P*(P-A)*(P-B)*(P-C)) END FUNCTION AREA PROGRAM MAIN IMPLICIT NONE REAL::L1,L2,L3,L4,L5,L6,L7,S1,S2,S3,S REAL::AREA READ*,L1,L2,L3,L4,L5,L6,L7 S1=AREA(L1,L2,L3) S2=AREA(L3,L4,L5) S3=AREA(L5,L6,L7) S=S1+S2+S3 PRINT*,S=,S END PROGRAM MAIN 方法一：使用函数实现 表达式方式调用 方法2：使用子例行子程序实现： SUBROU