An integral equation approach with a new MoM solution method is used to obtain current on the surface of a printed dipole antenna on a dielectric finite substrate, which is required for computation of antenna characteristic. For this purpose, a coupled tensor-volume/surface integral equations(TVSIE) is employed. To solve the tensor-volume integral equation, a mixed-domain MoM expansion using modified Legendre polynomial basis function and cubic modeling are adopted. For current modeling on the printed structure, Rao-Wilton-Glisson (RWG) basis function is used. The accuracy of the proposed approach is verified through comparison with the MoM solution based on the spectral domain Green's function for infinitely large substrate, the results of FDTD method, and experimental data.