j**u
发帖数: 6059 | 1 ☆─────────────────────────────────────☆
lzmaths (飞刀) 于 (Sat Mar 17 11:14:24 2007) 提到:
我想用子函数，但是子函数返回是 数组，要用Interface，
但是我写的程序老有错误，不知道能不能指点一下，很简单的一个程序，是求单位矩阵
的子程序。
谢谢！
program test
implicit none
interface
function eye(n)
implicit none
integer :: n
real :: eye(n,n)
end function
end interface
integer, parameter :: n=5
real :: A(n,n)
A=eye(n)
write(*,*) A
end
real function eye(n)
implicit none
integer :: n
real :: B(n,n)
integer :: i
B=0
forall(i=1:n) B(i,i)=1
eye=B
return
end
错误是：Err |
|
