求上三角矩阵问题

麦田守望者 · 发表于 2014-4-28 20:56:07

本帖最后由麦田守望者于 2014-4-28 21:40 编辑

[Fortran] 纯文本查看 复制代码

subroutine Upper(matrix)
implicit none
real::matrix(:,:)
integer::m,n
integer::i,j
real::e
m=size(matrix,1)
n=size(matrix,2)
do i=1,n-1
   do j=i+1,m
   e=matrix(j,i)/matrix(i,i)
   matrix(j,i:m)=matrix(j,i:m)-matrix(i,i:m)*e
   end do
end do
return
end subroutine Upper

n是列号，m是行号吧，这个matrix(j,i:m)=matrix(j,i:m)-matrix(i,i:m)*e感觉不对啊，为什么列从i:m呢？这句整个都没看懂，要是我写就写成matrix(j,i:n)=matrix(j,i:n)-matrix(i,i:n)*e，哪位大哥帮忙看看啊