They consider the solution of the pde
$$ Lu(x,y)=f(x,y). $$
So, $CF$ stands for solution of $Lu=0$ and
$PI$ the particular solution for $Lu=f$, then adding the
Solutions to get the general solution,
since $L$ is a linear operator. See here for a related technique.
Added: I think what they mean by $D$ is $\frac{d}{dx}$, $D'$ is $\frac{d}{dy}$, $\frac{1}{D^2}$ or ($D^{-2}$) is the second integral w.r.to $x$. Now,
$$ \frac{1}{D^2}\left(1-3\frac{D'}{D}-...\right)(x+y)= \frac{1}{D^2}\left( (x+y)-\frac{D'}{D}(x+y)-... \right)$$
$$=\frac{1}{D^2}\left( (x+y)-\frac{1}{D}(1)-0-0-... \right). $$
$$ \implies\frac{1}{D^2}\left(x+y-\frac{3}{D}(1)\right) = \frac{1}{D^2}(x+y-3x)=\frac{1}{D^2}(y-2x)=y\frac{x^2}{2}-\frac{x^3}{3} $$