In the present study, a compact analytical model is developed to determine the pressure drop of fully-developed, incompressible, and constant properties slip-flow through arbitrary cross-section microchannels. An averaged first-order Maxwell slip boundary condition is considered. Introducing a relative velocity, the difference between the bulk flow and the boundary velocities, the axial momentum reduces to the Poisson's equation with homogeneous boundary condition. Square root of area is selected as the characteristic length scale. Bahrami et al.'s model, which was developed no-slip boundary condition, is extended to cover the slip-flow regime in this study. The proposed model is a function of geometrical parameters of the channel: cross-sectional area, perimeter, polar moment of inertia and the Knudsen number. The model is successfully validated against existing numerical and experimental data from different sources in the literature for several shapes, including: circular, rectangular, trapezoidal, and double-trapezoidal cross-sections and a variety of gases such as: nitrogen, argon, and helium.