An isoparametric element method for the analysis of arbitrary dielectric optical wavefuides is described. This technique is based on the rigorous Ee-He variational formulation. The eight-node curve-sided quadrilateral element has been used. This type of isoparametric element is epecially suited to the study of the complex optical waveguides and optical waveguide devices with curved boundaries. As numerical examples, the step-index fiber, the graded-index fiber, the isotropic fiber coupler and anisotropic fiber coupler have been studied in thiS paper.