global GRID

% pre-processing
problem_setup

make_grid

[IX,M] = make_geometry(ax,bx,xc,yc,R,GEOM);

separate_boundaries

set_rhs

% ejiim solver
run_ejiim

% post-processing
visualise_solution