| sumising | R Documentation | 
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)
| niter | number of iterations | 
| numb | size of the square grid for the Ising model | 
| beta | Ising model parameter | 
returns a vector of 21 values for Z(\beta) corresponding to a regular sequence
of \beta's between 0 and 2
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.