Statistics of The effect of heterogeneity on the agreement between analytical and numerical solutions describing the recession flow in unconfined hillslope aquifers