Linear diffusion problems defined within irregular multidimensional regions are analytically solved through integral transforms, requiring numerical routines only for integration purposes, when a general functional boundary representation is considered. Auxiliary one-dimensional eigenvalue problems mapping the irregular region are applied with an integral transformation procedure so that the original differential SturmLiouville system gives place to an algebraic eigenvalue problem. The exact analytical inversion formula is then employed to yield the desired potential, explicitly, at any point within the domain. To allow for improved flexibility and further applicability, the related integration is simplified through an approximate boundary representation using lines connecting user provided points instead of the former exact representation of the irregular bounds, which is particularly advantageous when a functional description of the boundaries is not available. A cylindrical region test case with known exact solution is considered, and treated as an irregular region in the Cartesian coordinates system. Convergence behavior and error analysis are carefully undertaken and illustrated.