  This paper discusses the algorithms and implementations of three Mathematica
packages for the study of integrability and the computation of closed-form
solutions of nonlinear polynomial PDEs.
  The first package, PainleveTest.m, symbolically performs the Painleve
integrability test. The second package, PDESpecialSolutions.m, computes exact
solutions expressible in hyperbolic or elliptic functions. The third package,
PDERecursionOperator.m, generates and tests recursion operators.
