#
read cyl0;
#
#Bounded polynomials of the form fni*cos(n*theta)
#or fni*sin(n*theta)
#
F11:=-diff(F2,r);
F12:=-diff(F3,r):
F13:=-diff(F4,r):
F14:=-diff(F5,r):
F15:=-diff(F6,r):
F16:=-diff(F7,r):
F17:=-diff(F8,r):
F18:=-diff(F9,r):
#
F22:=-expand(simplify(diff(F13,r)-F13/r));
F23:=-expand(simplify(diff(F14,r)-F14/r)):
F24:=-expand(simplify(diff(F15,r)-F15/r)):
F25:=-expand(simplify(diff(F16,r)-F16/r)):
F26:=-expand(simplify(diff(F17,r)-F17/r)):
F27:=-expand(simplify(diff(F18,r)-F18/r)):
#
F33:=-expand(diff(F24,r)-2*F24/r);
F34:=-expand(diff(F25,r)-2*F25/r):
F35:=-expand(diff(F26,r)-2*F26/r):
F36:=-expand(diff(F27,r)-2*F27/r):
#