(*   BOUNDED NON-AXISYMMETRIC POTENTIALS IN CYLINFRICAL COORDINATES   *)

SetDirectory[NotebookDirectory[]];
<< cyl0.m

(* This file generates functions Fni such that the potentials 
Fni*Cos[n*θ] and Fni*Sin[n*θ] are harmonic.  *)

F11=-D[F2,r]
F12=-D[F3,r]
F13=-D[F4,r]
F14=-D[F5,r]
F15=-D[F6,r]
F16=-D[F7,r]
F17=-D[F8,r]
F18=-D[F9,r]

F22=-Expand[Simplify[D[F13,r]-F13/r]]
F23=-Expand[Simplify[D[F14,r]-F14/r]]
F24=-Expand[Simplify[D[F15,r]-F15/r]]
F25=-Expand[Simplify[D[F16,r]-F16/r]]
F26=-Expand[Simplify[D[F17,r]-F17/r]]
F27=-Expand[Simplify[D[F18,r]-F18/r]]

F33=-Expand[Simplify[D[F24,r]-F24/r]]
F34=-Expand[Simplify[D[F25,r]-F25/r]]
F35=-Expand[Simplify[D[F26,r]-F26/r]]
F36=-Expand[Simplify[D[F27,r]-F27/r]]