Hi I am having a problem with the output of my proc file. I wrote a proc as follows: Energy:=proc(n) local x,y,..; . . . . . plots[listdensityplot](H, labels=[Energy, Magnetisation]); W:=T->add(add(N[a,b]*b^2*exp(-a/(k*T)),a=Emin..Emax),b=Mmin..Mmax)/add(add(N[a,b]*exp(-a/(k*T)), a=Emin..Emax), b=Mmin..Mmax) plot(W(T)/n^4,T=0..5, labels=[Temperature, "