Derivative of Schur function In his answer to https://mathoverflow.net/questions/129854, R. Stanley says that the partial derivative (over the relevant x[i]) of the Schur function of a partition lambda of n equals the sum the Schur functions over the trimmed partions of lambda:
"Let $s_λ$ denote a Schur function and $p1=s1=x1+x2+\dots$. Then, $\frac{\partial\lambda}{\partial p_1}=\sum_\mu s_\mu$, where the $μ$'s are obtained by removing a single box from $λ$."
If I work this out for $\lambda = \{3,1\}$, I find for the derivative 
$12 X[3]+ 108 X[2,1]+ 60 X[1,1,1]$ (substituting $X[i,j,k]$ for terms $x[1]^i*x[2]^j*x[3]^k$) while the trimmed partitions $\{2,1\}$ and $\{3\}$ produce 
$4 X[3]+ 24 X[2,1] + 12 X[1,1,1]$. 
What am I missing? Can anyone provide an example of what is intended?
 A: Stanley uses a cryptic description whose interpretation is easy for the initiated, but not for others. Remember that the ring of symmetric functions is equipped with an inner product for which the Schur functions form an orthonormal basis (you can take that as its definition, although there are also other descriptions). Now the operation written $\frac\partial{\partial p_1}$ is just the adjoint (in the sense of linear mappings) of the operation of multiplication by $p_1=s_1=\sum_i x_i$. Since that multiplication applied to $s_\lambda$ gives the sum of all $s_\mu$ where the diagram of$~\mu$ can be obtained by adding a square to that of$~\lambda$, its adjoint applied to $s_\lambda$ gives the sum of all $s_\mu$ where the diagram of$~\mu$ can be obtained by removing a square from that of$~\lambda$.
A way to view this operation as applied to any symmetric function $f$ is to single out one of the variables in which $f$ is  a symmetric function (since it is symmetric it doesn't matter which one), call that variable $y$ and the others $x_1,x_2,\ldots$; observing that $f$ is certainly symmetric in the $x_1,x_2,\ldots$, write it as a polynomial in $y$ with as coefficients symmetric functions in the $x_1,x_2,\ldots$, and then take the coefficient of the first power $y^1$ of$~y$  as outcome (we've used up one variable, but given the infinite supply this makes no difference; this is the Hilbert hotel in operation). In practice one doesn't need infinitely many variables: one can write $f$ as symmetric polynomial in variables $x_1,\ldots,x_n$ where $n$ is at least $\deg f$, and choose $y=x_n$ as "ultimus inter pares" which avoids renumbering (the degree conveniently drops by $1$, so $x_n$ is expendable).
If one takes the description of $s_\lambda(x_1,\ldots,x_n)$ in terms of semistandard tableaux, one can partition the set of its terms according to the set of squares containing an entry$~n$ in the tableau. The parts of this partition of terms each define a product $s_\mu(x_1,\ldots,x_{n-1})x_n^d$, where $d$ is the number of squares in the fixed set containing an entry$~n$, and the diagram of $\mu$ consists of the remaining squares. It follows that the coefficient of $x_n^d$ is the sum of all such $s_\mu(x_1,\ldots,x_{n-1})$, where the diagram of$~\mu$ is obtained from that of$~\lambda$ by removing a horizontal strip of size$~d$. In particular, our operator maps $s_\lambda$ to the sum of all $s_\mu$ whose diagram is obtained from that of $\lambda$ by removing one corner.
To see that the operator coincides with formal differentiation with respect to $p_1$ when applied to a polynomial in $p_1,p_2,\ldots$, interpreted as power sum symmetric functions, one can either prove that it is a derviation which vanishes on $p_2,p_3,\ldots$ and sends $p_1\mapsto1$, or one can directly apply our description to a product $p_1^mf$ where $f$ is a factor on which the operator vanishes (a monomial in $p_2,p_3,\ldots$), using $p_1(x_1,\ldots,x_n)=p_1(x_1,\ldots,x_{n-1})+x_n$ and expansion of $(p_1(x_1,\ldots,x_{n-1})+x_n)^m$, giving as result of the operation $mp_1^{m-1}f$.
