A computer routine for personal computer(PC/AT class) is developed to analysize the mode characteristics of silica based optical waveguides whose cross sections are of semi circular and other typical shapes. The basic algorithm used in the routine is to convert the wave equation into a matrix equation by expanding the wave function in terms of simple harmonic functions. The matrix elements are a set of overlap integrals of sinusoidal funtions with appropriate weight given by the distribution of refractive index over the waveguide cross section. The eigenvectors and eigenvalues of the matrix is then computed via diagonalization. We explain some practical problems that arises when implementing the algorithm into the routine. By using this routine we analyze the mode characteristics of silica based optical waveguides of semi circular and some other typical cross sections.