powered by
This function implements a path sampling approximation of the normalising constant of an Ising model with a four neighbour relation.
sumising(niter = 10^3, numb, beta)
returns a vector of 21 values for \(Z(\beta)\) corresponding to a regular sequence of \(\beta\)'s between 0 and 2
number of iterations
size of the square grid for the Ising model
Ising model parameter
isingibbs,isinghm
isingibbs
isinghm
Z=seq(0,2,length=21) for (i in 1:21) Z[i]=sumising(5,numb=24,beta=Z[i]) lrcst=approxfun(seq(0,2,length=21),Z) plot(seq(0,2,length=21),Z,xlab="",ylab="") curve(lrcst,0,2,add=TRUE)
Run the code above in your browser using DataLab