Integrable variable coefficient equations describe the real world in many fields of physical and engineering sciences. Many researchers are devoted to discussing these equations by utilizing different methods ref.  - . In ref.  , Dai and Jeffrey extended the dressing method to a generalized version for solving nonlinear evolution equations associated with matrix spectral problems and variable coefficient cases, in which a key is that variable coefficient dressing operators are transformed to different variable coefficient ones. By using the generalization, we have studied integrable variable coefficient coupled Hirota equation in ref. . In ref.  , integrable variable coefficient Manakov model and cylindrical NLS equation are discussed in detailed, respectively. In ref. , we developed the generalized dressing method to the discrete system and an integrable variable coefficient Toda equation is researched. Recently, the dressing method is extended to a matrix Lax pair for Camassa-Holm equation in ref. , in which interactions between soliton and cuspon solutions of the system are studied. The dressing method as nonlinear superposition in Sigma models has been researched by Dimitrios Katsinis et al. in ref. . Multi-lump solutions of KP equation with integrable boundary are discussed in ref.  by using the generalized dressing method. Nabelek etal. in ref.  studied Kaup-Broer system and derived its solutions.
In the present paper, we extend the generalized dressing method to discrete operators similar to ref. . Through direct calculations, we derive a new integrable variable coefficient Toda equation
where, the coefficient is related to n, . Equation (1.1) is an extension of the well known two-dimensional Toda equation. We will construct one soliton solution of (1.1).
The present paper is organized as follows. In Section 2, we obtain a new integrable variable coefficient Toda equation based on the generalized dressing method. In Section 3, as an application, we derive one soliton solution of (1.1) by utilizing the separation of variables.
2. Integrable Variable Coefficient Toda Equation
In this section, we first summarize the variable coefficient version of the dressing method. We extend the generalized version of the dressing method to discrete systems and derive different integrable cylindrical Toda lattice equations by choosing different operators.
First, we consider three linear differential difference operators ref. 
Similar to the generalized dressing method application to continuous system, we introduce the triangular factorization about the operator “ ”
where is the identity operator, for and for . It is assumed that
for all . For convenience, we denote , . The discrete Gelfand-Levitan-Marchenko (GLM) equation can be obtained from (2.2), which reads in ref. 
We introduce two differential-difference operators and defined by
where is the shift operator of the discrete variable n, defined by , , t and y are continuous variables.
The dressing operators and can be derived from the relations
Similar to a theorem ref.  for continuous systems, it can be proved that and are differential-difference operators. For sake of simplicity, we denote .
We write the dressing operators
Acting on function on (2.6) and with aid of (2.8), which is reduced to
from which, comparing coefficient of , we have
Letting , with aid of (2.7) and (2.9), we have
from which, comparing coefficient of , we have
and we derive
The following theorem in ref.  is an extension of original dressing method, which can yield a wide range of integrable variable-coefficient nonlinear evolution equations.
Theorem: If the operators and satisfy a relation
where are arbitrary functions of , then their corresponding dressing operators will satisfy the relation
Proof: According to (2.6), (2.7) and (2.13), we can give simple proof as follows through simple calculation. In fact,
Actually, variable-coefficient Toda equations are obtained from (2.14). From (2.14), we derived
then the above Equations (2.15) and (2.16) are reduced to
According to (2.19), we assume that
then (2.18) is reduced to a new integrable variable coefficient Toda equation
Let , , then the above equation is reduced to a new 2 + 1 dimensional Toda lattice equation
The above equations are new and different to classical Toda lattice equation in ref.     . Because the coefficient of equation is related to n, this is an important physical meaning.
3. Explicit Solution of Integrable Variable Coefficient Toda Equation
In this section, we shall use the generalized dressing method to construct explicit solutions of the variable coefficient Toda Equation (2.21). Using the relation , we have
Assume that (3.1) and (3.2) have N-soliton solutions in the form of separation of variables
moveover, we suppose that
Substituting (3.3) and (3.4) into the GLM (2.3) yields that
where L is defined by
and is Kronecker’s delta.
In what follows, we will obtain one soliton solution of (2.21). First, we give separation of variables solutions for in (3.3) and (3.4),
From (3.5), we derive
with , using (2.17), we have
Under transformation , we derive one soliton solution of (2.21)
The authors thank the authors of the references. The work described in this paper is supported by National Natural Science Foundation of China.
The data used to support the findings of this study are available from the corresponding author upon request.
The work described in this paper is supported by National Natural Science Foundation of China (Grant No.11301149).
Yanan Huang and Ting Su do derivation and calculations. Junhong Yao mainly draw soliton solution picture.
 Hou, S. and Liu, X.D. (2005) A Numerical Method for Solving Variable Coefficient Elliptic Equation with Interfaces. Journal of Computational Physics, 202, 411-445.
 Sun, Z.Y., Gao, Y.T., Liu, Y., et al. (2011) Soliton Management for a Variable-Coefficient Modified Korteweg—De Vries Equation. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics, 84, Article ID: 026606.
 Jeffrey, A. and Dai, H.H. (1990) On the Application of a Generalizzed Version of the Dressing Method to the Integration of Variable-Coefficient KdV Equation. Rendiconti di Matematica e delle sue Applicazioni, 10, 439-455.
 Su, T., Dai, H.H. and Geng, X.G. (2009) On the Application of a Generalized Dressing Method to the Integration of Variable-Coefficient Coupled Hirota Equations. Journal of Mathematical Physics, 50, 226.
 Su, T., Ding, G.H. and Fang, J.Y. (2014) Integrable Variable-Coefficient Coupled Cylindrical NLS Equations and Their Explicit Solutions. Acta Mathematicae Applicatae Sinica, 30, 1017-1024.
 Dai, H.H. and Su, T. (2010) The Generalized Dressing Method with Applications to the Integration of Variable-Coefficient Toda Equations. Proceedings of the Estonian Academy of Sciences, 59, 293-298.
 Nabelek, P.V. and Zakharov, V.E. (2020) Solutions to the Kaup? CBroer System and Its (2 + 1) Dimensional Integrable Generalization via the Dressing Method. Physica D Nonlinear Phenomena, 409, Article ID: 132478.
 Date, E. and Tanaka, S. (1976) Analogue of Inverse Scattering Theory for the Discrete Hill’s Equation and Exact Solutions for the Periodic Toda Lattice. Progress of Theoretical and Experimental Physics, 56, 457-465.