An approximated analytical solution of the Navier-Stokes equations describing the laminar flow through a three-dimensional rectangular channel with wall suction or injection is presented. The solution is developed in terms of a stream function whose variation is described by a separation approach in combination with the polynomial approximation of the spatial velocity profiles. The approximated profiles agree well with simulated results obtained from the commercial code ANSYS CFX 10.