FORTRAN复习课.ppt

9.8 外部子程序/示例2/外部子例行程序3 9.8 外部子程序 [9.23]设计和编写矩阵相加的外部子例行程序。 解：已知矩阵Am×n、Bm×n,求Cm×n=Am×n+Bm×n。 主程序单元为： 输入矩阵数据的外部子例行程序为： 输出矩阵数据的外部子例行程序为： 计算矩阵之和的外部子例行程序为： 输入两个3×4矩阵数据(整数)： 12 34 41 55 10 10 10 10 A：10 21 33 18 B：20 20 20 20 50 42 31 14 30 30 30 30 输出结果矩阵为： 22 44 51 65 30 41 53 38 80 72 61 44 !外部子例行程序 !计算矩阵之和,A和B为已知矩阵，C为结果矩阵 SUBROUTINE matrix_add(A,B,C,n1,n2) INTEGER n1,n2,i,j !声明3个n1×n2数组 INTEGER,DIMENSION(n1,n2) :: A,B,C DO i=1,n1 DO j=1,n2 C(i,j)=A(i,j)+B(i,j) ENDDO ENDDO END SUBROUTINE 概述 定义 调用 EXTERNAL 示例1 示例2 第九章 函数与子例行程序 * aaaaaaaaaaa * List-directed 表控 * INTEGER :: sum,n=2,I,x REAL av sum=0; DO I=1,n READ *,x sum=sum+x ENDDO av=sum/n WRITE(*,*) 平均成绩为：,av * INTEGER sum,n,x REAL av sum=0;n=0 READ *,x DO WHILE (x .ne. -1) n=n+1; sum=sum+x READ *,x ENDDO av=sum/n WRITE(*,*) 平均成绩为：,av * INTEGER sum,n,x REAL av sum=0;n=0 READ *,x DO n=n+1;sum=sum+x;READ *,x IF (x==-1) EXIT ENDDO av=sum/n WRITE(*,*) 平均成绩为：,av M=1;N=5;L=1 DO I=M,N,L PRINT*,I,M,N,L M=M+1;N=N*I;L=2*L+I ENDDO PRINT*,I DO I=1,100,2 PRINT*,I IF (I.GE.9) EXIT ENDDO I=1 DO WHILE (I.LE.100) PRINT *,I IF (I.GE.9) EXIT I=I+2 ENDDO I=1 DO PRINT *,I IF (I.GE.9) EXIT I=I+2 ENDDO END DO I=2,5 PRINT (1X,I3\)“,I IF (I3) CYCLE PRINT (1X,I3\),I END DO PRINT *,‘结束！’ END * 矩阵乘积：matmul; dot_product; 8.7 动态数组/分配 integer,allocatable::a(:) read*,n allocate(a(n)) !allocate(a(n),stat=ierr) do i=1,n a(i)=i enddo sum=0 do i=1,n sum=sum+a(i) enddo av=sum/n print*,av deallocate(a) end 8.7 动态数组 概述 声明 分配 释放 示例 第八章 数组 1. 概述 2. 标准子程序 3. 语句函数 4. 内部子程序 5. 子程序形参与实参之间的数据传递关系 6. 递归子程序 7. 外部子程序 第九章 函数与子例行程序 《FORTRAN90程序设计》课多媒体课件 9.1 概述/说明 9.1 概述 相关基本概念： 子程序（过程 subprogram）： 函数 (function) 子例行程序(subroutine)。 函数：只产生一个运算结果并与函数名相关 联，通过操作数调用执行的子程序。 子例行程序：完成一个处理过程，产生多个处 理结果，通过CALL语句调用执行 的子程序。 FORTRAN 90子程序分类如图9-3所示。 概述 情况1 情况2 情况3 说明 第九章 函数与子例行程序 9.3 语