powered by
sumising(niter = 10^3, numb, beta)
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