The solitary wave, so-called because it often occurs as a single entity and is localized, was first observed by J. Scott Russell on the Edinburgh-Glasgow Canal in 1834. It is known that many nonlinear evolution equations have soliton solutions, such as the Korteweg de Vries equation, the Sin-Gordon equation, the nonlinear Schrödinger equation, the Kadomtsev-Petviashvili equation, the Davey-Stewartson equation, etc. In order to study the property of nonlinear evolution equations, methods are developed to derive solitary wave solution or soliton solution to nonlinear evolution equations. Some of the most important methods are the inverse scattering transformation (IST)  method, the bilinear method  -  , symmetry reduction method  , the Bäcklund or Darboux transformation method  and so on. Having soliton solutions is one of the basic integrable properties of nonlinear evolution equations.
In this paper, we are interested in the general expression of N-soliton solution to the (2 + 1)-dimensional KdV equation,
which was first derived by Boiti et al. by using the idea of the weak Lax pair  . This system can also be obtained from the inner parameter-dependent symmetry constraint of the KP equation  . Recently, the dromion solutions and some exact solutions are studied by Lou and Wazwaz respectively  -  . While as for the uniformed ex- pression of its N-soliton solution is unknown yet.
In this article, we’ll study the N-soliton solution to the (2 + 1)-dimensional KdV system (1). A compact form of the N-soliton solution to Equation (1) is obtained by means of pfaffian technique, which is given in section 2. Conclusion and further discussions are given in section 3.
2. N-Soliton Solution to the (2 + 1)-Dimensional KdV Equation
Given a nonlinear evolution equation, if it has 3-soliton solution, then this equation is of great possibility of having N-soliton () solution. Pfaffian technique is one of the methods that can help us to determine whether the evolution equation has multi- soliton solutions or not. In this section, we first review some properties of pfaffian.
Pfaffians are antisymmetric functions with respect to its independent variables
An n-th order pfaffian can be defined inductively by the expansion rule 
where denotes the absence of letter j. For example, when n = 2, we have
There are various kinds of pfaffian identities. In this article, we just introduce the so-called Plüker relation for pfaffians 
which we are going to use. Hereafter, we let denote pfaffian
2.2. N-Soliton Solutions
The Hirota form of the (2 + 1)-dimensional KdV system (1) is
which is obtained by the dependent variable transformations
Here the Hirota bilinear operator is defined by
with n and m are arbitrary nonnegative integers.
In  , the 3-soliton solution to the (2 + 1)-dimensional KdV system (3) is obtained
via the perturbation method. It claims that the N-soliton solutions for can also be obtained by using perturbation method, but the explicit expression of the multi- soliton solution is not given.
In this article, we’ll study the the multi-soliton solution to Equation (3) using the pfaffian technique   . A compact form of the N-soliton solution is given in terms of N-th order pfaffian.
Proposition 1. If the t-function f of the (2 + 1)-dimensional KdV system (3) is given by the pfaffian function
whose entries, for, are defined by
then this particular pfaffian function (6) gives an N-soliton solutions to the (2 + 1)-di- mensional KdV system (3).
Proof. In the following, we will prove that the pfaffian function (6) satisfies the (2 + 1)-dimensional KdV Equation (3). By defining “differential operators” ()
we obtain the following differential formulae
In order to find the pfaffian expression for the differential functions with derivative of variable y, we need to define another letter 
Then we have
Substituting formulae (9) and (11) into the right hand side of Equation (3), we obtain nothing but the Plücker relation for pfaffians (2)
where denotes. Therefore the pfaffian function (6) solves the (2 + 1)-dimensional KdV system (3).
Note that in order to derive the differential formulae of the pfaffian function (6), we have to define another extra letter besides the “differential operators”. The multi-soliton solution to the nonlinear (2 + 1)-dimensional KdV system (1) can be obtained by substituting pfaffian function (6) into the dependent variable transfor- mation (4) directly.
In this article, a compact form of the multi-soliton solution to the (2 + 1)-dimensional KdV system is given via the pfaffian technique. As one can see, the key point of the proof is to derive suitable expressions of the differential formulae of pfaffian t-function f. It is worth pointing out that the method used in this article is different as the one for the proof of the BKP equation, which the differential formulae of the pfaffian t-fun- ction depend only on the “differential operators”.
We thank the editor and the referee for their comments. This research work is supported by the National Science Foundation of China (Grant Numbers 11271362, 11271266 and 11501510) and a President’s Grant from the University of Chinese Academy of Sciences. These supports are greatly appreciated.