This study presents a new form of velocity distribution in laminar liquid flow in rectangular microchannels using the eigenfunction expansion technique. Darcy friction factor and Poiseuille number are also obtained analytically. Due to the symmetry of the solutions, the effects of changing the aspect ratio from 0 to ∞ are also discussed. Using finite element method (FEM), the obtained analytical results are further compared with the 3D numerical simulations for the rectangular microchannels with different range of aspect ratio and pressure gradient, and excellent agreements were found. These findings provide additional insights in interpreting the results of the pressure-driven flows in finite aspect ratio microchannels, in which very precise comparison with the macroscale theory is crucial.