The optimal design code of an axial flow pump has been developed to determine geometric and fluid dynamic variables under hydrodynamic as well as mechanical design constraints. The design code includes the optimization of the complete radial distribution of the geometry by determining the coefficients of 2$^{nd}$ order polynomials to represent the three-dimensional geometry. The optimization problem has been formulated with a nonlinear multivariable objective function, maximizing the efficiency and stall margin, while minimizing the net positive suction head required. Calculation of the objective function is based on the mean streamline analysis and through-flow analysis using the present state-of-the-art model. The optimal solution is calculated using the penalty function method in which the genetic optimizer is employed. The optimized efficiency and design variables are presented in this paper as a function of non-dimensional specific speed in the range, 2$leq$ $n_{s}$ $leq$10. The results can be used in preliminary design of axial flow pumps.