Port the ising pseudocode on page 349 to the program of your choice and show interesting plots.

8.2 The Ising Model of a Ferromagnet Monte Carlo simalation of a 2D Ising model using the Metropolis algorithm Width of square lattice Temperature in units of ek 349 program ising ize-10 T-2.5 initialize for iteration -1 to 100*size 2 do Main iteration loop 1-int (rand*size 1) Choose a random row number and a random column number Compute U of hypothetical flip If Hipping reduces the energy then flip it int (rand size1) deltaU,j,Ediff.) if Ediff0 then colorequare(i,j) else it randexpC-Edift/T) then otherwise the Boltzmann factor gives the probability of flipping colorsquare(i,j) end if end if next iteration Naw go back and start over end progran subroutine delta(i,j,Ediff) Compute AU of fipping a dipole (note periodic boundary conditions) if 1 then top - )lse top s(i-1,j) t size then bottom-(1,) olse bottom- s(1+1.) if j-1 then left s(1,size) else ieft-s(i,j-1) if j-size then right 8(1,1) else right " s(1,j+1) Edif - 2*8ij) (top+bottom left+right) and subroutine subroutine initialize Initialize to a random array fori 1 to size for j1 to size if rand