矩阵算法合集问题续II

PotsyYZhou

按提示， 变换了一下格式如下：

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

subroutine eulerian_inverse ( n, a )
  do j = 1, n
    do i = 1, n
      if ( i .lt. j ) then
        a(i,j) = 0.0D+00
      else if ( i .eq. j ) then
        a(i,j) = 1.0D+00
      else if ( j .lt. i ) then
        dot = 0.0D+00
        do k = j, i - 1
          dot = dot + b(i,k) * a(k,j)
        end do
        a(i,j) = - dot
      end if
    end do
  end do
  return
end

      还是不行，报错如下：
main.f90:14269:8:14269 |         a(i, j) = 0.0D+00      |        1Error: Unclassifiable statement at (1)
main.f90:14281:8:14281 |         a(i, j) = -dot      |        1Error: The function result on the lhs of the assignment at (1) must have the pointer attribute.
？？？！

li913

a需要定义为数组。