Part_II_数值分析程序设计.docVIP

数值分析程序设计 Part II Fortran程序设计 1 线性代数 1.1 矩阵的加法和乘法 矩阵作为数值分析的重要工具，在数值计算中具有非常重要的作用，在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) INTE