from sage.all_cmdline import *
from phi4 import Simulation
from phase_plot import PhasePlot
plt = PhasePlot(32, 32, equilibrate=1512, iterations=1512, observable=Simulation.abs_phi)
plt.mass_squared = srange(-1.5, 0, 0.1)
plt.coupling_constant = srange(0.1, 1, 0.1)
gfx = \
plt.plot3d() + \
text3d(r'lambda', [0, 0.75, 1]) + \
text3d(r'mu squared', [-0.5, 0, 0.2]) + \
text3d(r'abs(phi)', [-1.5, 0, 2])
gfx.show()
