function [ha, pos] = tight_subplot(Nh, Nw, gap, marg_h, marg_w) ...