Imaging formulas of Gaussian beams in nonlinear conical GRIN rod lenses are derived from the variational method and the ABCD law. It is shown that the imaging properties and the ray paths of such lenses depend on the optical power. The imaging formulas reduce to the known results when n2=0 ( n2: the nonlinear coefficient) and under the geometrical optics approximation.