Shapes and lay-up configurations of laminated composite shallow shells are optimized simultaneously to maximize the fundamental frequency by a simple genetic algorithm method. The surface shape is defined by a cubic polynomial and this makes it possible to express a variety of surfaces with inconstant curvature radii by varying the values of the coefficients for each polynomial term. The coefficients and the lay-up configurations of the laminated shells are directly employed as design variables, and constraints are imposed on the coefficients and curvature radii to maintain the shallowness of the shells. The frequencies are calculated by using the Ritz method due to its flexibility with surface shapes. The results of the present analysis agreed well with experimental and finite element analysis results in terms of frequencies and mode shapes. The obtained optimum solutions resulted in higher fundamental frequencies than for shells with commonly emplyed surface shapes and lay-up configurations.