“…We employ the expansion–collocation method as described by Farina and Martin 16 . To tackle the weaker singularity we have used the Korobov–Conroy number theoretic method during the computations of the integrals as proposed by Farina et al 20 . In cylindrical polar coordinates and and the disc is given as Let, , then Therefore, Equation () can be written as …”