The problem of optimization of fed-batch fermentations using the substrate feed rate as the control variable is singular in nature. Previous approaches, including the boundary condition iteration method and transformation to a nonsingular problem using a different control variable, do not work well for solving optimization of systems governed by more than four differential equations. The applicability of a first-order conjugate gradient algorithm for optimizing fed-batch fermentions was tested for systems of varing complexity. This approach does not need any variable transformation of a priori knowledge of the control arc sequence. Constraints on the feed rate are handled in a simple and direct manner. The algorithm worked very well for three, four, and five-dimensional singular systems. The correctness of the optimal profile was judged by observing the variation in the sign of the gradient of the Hamiltonian. The gradient was found to be zero during the singular period and had the appropriate sign on the boundary arcs. The optimization method based on conjugated gradient approach can be complementary to the boundary condition iteration method for determination of the exact optimum profile.