Some new solutions of the first kind of elliptic equation and formula of nonlinear superposition of the solutions are given. Combing the solutions with function transformation to construct Jacobi-like elliptic function exact solutions, a number of soliton-like solutions and triangular solutions of (3+1)-dimensional Zakharov-Kuznetsov equation with variable coefficients are given with the help of symbolic computation system Mathematica.