  The rules to write out any one of the linearly independent functions
belonging to the infinite set of those in polynomial form that satisfy the
sourceless Grad-Shafranov equation as stated in the toroidal-polar coordinate
system are established. It is found that a polynomial solution even in the
poloidal angle is given by the product of an integral power of the radial
coordinate variable by a complete polynomial of equal degree in this same
variable with angular-dependent coefficient functions that are linear
combinations of a finite number of Chebyshev polynomials in the cosine of the
poloidal angle, the numerical coefficients of these being expressed in terms of
the binomial numbers of Pascal's arithmetic triangle. Tables of the ten
polynomial solutions of the lowest degrees are provided in variables of the
toroidal-polar and of the cylindrical coordinate systems.
