求互相关的fortran小程序

do i=1,4 sumx=0; do ii=1,n read(i,*)ux(i,ii) sumx=sumx+ux(i,ii) end do ax(i)=sumx/n; sumx=0 do ii=1,n sumx=sumx+(ux(i,ii)-ax(i))**2 end do vx(i)=sqrt(sumx/(n-1)) end do do itao=-id,id ……