

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

<< cyl0.m

(* This file generates functions Fni such that the potentials 
    Fni*Cos[n*\[Theta]] and Fni*Sin[n*\[Theta]] 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]]
