Semi-analytic numerical method for unsteady flow in leaky confined aquifers