mx=116; % mx and my should be multiples of 58
my=116;
epack=.1; % small for uniform grid, ~1 for more points at edge
alpha=0.0; % specify linear lambda profile,
beta=0.0; % coefficient of the second derivative term of lambda(psi) (NOT p/B^2)
rms_tol=1e-9;
max_iter=300; % total iterations before quitting
restart=0; % set to 1 if no grid parameters have been changed from previous
           % case, else set to zero.
remeshed=0;