The mathematical model of light propagation in a planar gradient optical waveguide consists of the Maxwell’s equations supplemented by the matter equations and boundary conditions. In the coordinates adapted to the waveguide geometry, the Maxwell’s equations are separated into two independent sets for the TE and TM polarizations. Each of the systems can be transformed to a second order ordinary differential equation. The boundary conditions forMaxwell’s equa- tions are reduced to two pairs of boundary conditions for the obtained equations. Thus, the problemof describing a complete set of modes in a regular planar waveg- uide is formulated in terms of an eigenvalue problem. For each polarization there are three types of waveguide modes: guided modes, substrate radiation modes, and cover radiation modes. In this work we implement the numerical-analytical calculation of all types of waveguide modes.