To solve the equilibrium, save all files in some directory, and type "main" on matlab command windos

Current parameters compute equilibrium with a coarse grid and small number of iterations

To get a more accurate solution, increase parameter "iternum" and decrease grid parameters "dxinit" and "dsinit"


