The method combing the function transformation and second order homogeneous linear ordinary differential equation with constant coefficients is proposed. With the help of symbolic calculation system Mathematica, the new soliton-like solutions of (3+1) dimensional Burgers equation with variable coefficients are constructed, which consists of the exponential function, trigonometric function and rational function.