An analytical (strong form) solution to the boundary value problem of moderately thick shallow cylindrical panels, with arbitrary laminations is presented. A double Fourier series approach is used to solve five second-order highly coupled partial differential equations that arise from the shallow shell formulation based on popular Donnell-Mushtari-Vlasov shell theory, and the first-order shear deformation-based through-thickness theory. An admissible boundary condition is considered to obtain numerical results that constitute the study of convergence of displacements, and moments; and spatial variations of them presented in the form of contour plotting for Various parametric effects. These, hitherto unavailable, analytically obtained numerical results should serve as base-line solutions for future comparisons of popular approximate methods such as finite element, finite difference, Galerkin approach, Rayleigh-Ritz method, collocation method, least-squares method, and experimental results, for the case of arbitrary laminations.