We present a largely analytic series-based approach for steady saturated-unsaturated flow problems in irregular shaped porous domains. We consider a porous medium with irregulat shaped boundary, both at the surface and at the base. Over the ground surface is steady rainfall with any excess allowed to freely run off without ponding. We provide a variational formulation which has explicit dependence on the position of the saturated-unsaturated interface (phreatic surface) and derive a series solution for the integrand of the penalty functional. A simple direct numerical method may then be applied to minimize the functional. This approach has several advantages. It affords a realistic description of the water distribution in both zones and an accurate location of the phreatic surface. The exact separated solutions as basis functions remove much of the arbitrariness that applies to spectral methods. No spatial discretisations are necessary and the global solution errors are readily estimated from maximum principles. We apply the solution to a canonical hillslope geometry, viz. an inclined permeable region whose cross-section is a long, thin parallelogram. We find that the height of the phreatic surface is relatively independent of the sorptive number. We prove that the phreatic surface is tangential to the soil surface at the point where the phreatic surface meets the soil surface.