基本矩阵运算

https://www.cnblogs.com/xxfss2/archive/2011/11/15/2249981.html

去博客设置页面，选择一款你喜欢的代码片高亮样式，下面展示同样高亮的 代码片.

// An highlighted block
program Console4
    implicit none
    parameter N=5 
    integer i,j
    real  array1(2,N),array2(N,2),array3(2,2)
    real array4(2,3),array5(2,3),array6(2,3)
    array1(1,1)=3.2;array1(1,2)=4.5;array1(1,3)=7.56;array1(1,4)=6.34;array1(1,5)=34.2;
    array1(2,1)=44.3;array1(2,2)=12.3;array1(2,3)=1.236;array1(2,4)=12.2;array1(2,5)=9.5; 
   !对array1 求转置矩阵 
    array2=transpose(array1)
   call PrintArray(array1) 
   call PrintArray(array2) 
   !对array1和array2求乘积 
    array3=matmul(array1,array2)
   call PrintArray(array3)  
   array4=3
   array5=2
   array6=array4+array5
   call PrintArray(array6)  
   array6=array4-array5
   call PrintArray(array6)  
   !array6= lincg(array5)
    read *,i

    end program Console4

subroutine PrintArray(array)
    real array(:,:)
   integer length,i
   length=size(array,1) 
   do i=1,length
        print *,array(i,:)
   end do
   print *,''
end subroutine