This peper reports an flexible, interactive and atuomatic system, based on the vector finite element method, and implemented within IMSL'S PDE/Protran problem solving environment. Numerical results for all Propegsting modes in an anisotropic dielectric rectangular waveguide are presented and verified through compeison with an earlier study (where, in fact, fewer modes were found). Furthermore, in comparison with other studies, the evaluation of accuracy of the method is demonstrated in terms of both propagation constants as well as modal field profiles. These results also clearly illustrate the speed and the ease of use of the modified WAVEGIDE program.