program identity implicit none integer ,dimension(4,4) ::ra integer i ra=0 do i=1,4 ra(i,i)=1 end do print *,ra end program identity