ix(9)=[1,2,3,1,2,3,1,2,3]
iy(9)=[1,1,1,2,2,2,3,3,3]
G(i,j) = (ix(i)-ix(j))*ex + (iy(i)-iy(j))*ey

necrohan 发表于 2020-12-14 15:30
不知道你要的是不是这个意思

[mw_shl_code=fortran,true]program test

不是，我要的不是求Gij的值，我是想要每个位置的坐标

不知道你要的是不是这个意思

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

program test
implicit none
integer i,j,n
real ex,ey,g,tmp

n=9;
ex=1;ey=1

do i=1,n-1
  tmp=g(i,ex,ey)-g(i+1,ex,ey)
  write(*,*)i,tmp
enddo

end

real function g(i,ex,ey)
implicit none
integer i
integer row,col
real ex,ey
row=mod(i,3)
col=i/3+1
if(row==0)then
  row=3
  col=col-1
endif
!write(*,*)i,row,col
g=row*ex+col*ey
end

necrohan 发表于 2020-12-14 09:48
这是个3*3的矩阵吧

应该算是吧，应该是两个3*3的点坐标相减

这是个3*3的矩阵吧