Part_II_数值分析程序设计.docxVIP

PAGE PAGE 10 数值分析程序设计 Part II Fortran 程序设计 线性代数 矩阵的加法和乘法 矩阵作为数值分析的重要工具，在数值计算中具有非常重要的作用，在Fortran 语言中利用数组来表示矩阵。在 Fortran 90 中可以直接对数组进行计算， 一个命令就可以完成矩阵的加法、减法。 Real a(m,n), b(m,n),c(m,n) C=a+b ！矩阵加法 C=a-b ！矩阵减法 在 Fortran 77 中必须利用循环完成矩阵的加减法： do r=1,m do c=1,n c(r,c)=a(r,c)+b(r,c) end do end do 矩阵的乘法，在 Fortran 90 中不能直接用乘号来做，必须调用库函数。 C=a*b ！这个命令是做 c(I,j)=a(I,j)*b(I,j) C=matmul(a,b) ！库函数 matmul 可以做矩阵相乘 在 Fortran 77 中则要自己写矩阵乘法程序代码进行计算，下面是 Fortran 77 固定格式编写的矩阵相乘程序代码： C C 矩阵乘法范例 C By Perng 1997/9/17 PROGRAM MATMUL_DEMO IMPLICIT NONE INTEGER N 数值分析程序设计 PARAMETER(N=3) INTEGER A(N,N) ! MATRIX A INTEGE