Integral equations method to solve Lapace equation numerically with applications to torsion problems for multiply connected and compositive cylinders