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