This paper presents a novel theoretical and numerical approach for an infinitesimal current source (ICS) located in a planar isotropic multilayer medium. Using the mixed-potential integral equation (MPIE) formulation for depicting the electromagnetic disturbance created by the ICS, a detailed definition of Green's functions of Lorenz potentials and fields is provided in this paper. The proposed Green's functions are valid for the considered multilayer isotropic medium, which can have arbitrary layer parameters. This paper also analyzes two commonly observed special cases of the multilayer medium -the multilayer soil including air and the multilayer lossless dielectric -and the proposed equations are modified to meet the requirements of the medium. Green's functions can be obtained from the systems of linear equations proposed in this study. In comparison to other approaches, the advantage of the proposed procedure is that the solutions of the equations are immediately obtained in any field layer of the multilayer medium. In addition, the proposed system of linear equations can be solved easily using well-known numerical computation methods. Furthermore, this paper offers an alternative way of obtaining Green's functions, which are closed-form expressions for the kernels of spectral-domain Green's functions.