in dz_from_phi, included alternative code that uses fzero for root finding (result: slow)
in find_regions, included the old code using bwconncomp as alternative code
removed computation and use of n2_ns in optimize_surface_exact.m
removed some unnecessary lines in the construction of lsqr matrix