  The infinite many symmetries of Davey-Stewartson (DS) system are closely
connected to the integrable deformations of surfaces in a four-dimensional
space. In this paper, we give a direct algorithm to construct the expression of
the DS hierarchy by two scalar pseudo-differential operators involving partial
derivatives.
