The determination of the geoid requires knowledge of the actual topographical density distribution, which is often approximated using the constant density of 2.67 g/cm³ due to the challenges in determining the actual density. This approximation introduces inaccuracies in reduced gravity, impacting the geoid. In this paper, geoidal heights are computed for a region within 36° – 39° latitude and 252° – 255° longitude using Stokes-Helmert's method, a gravimetric solution developed at the University of New Brunswick (UNB). The computation includes using different density models, including anomalous density. Direct density effect (DDE) and Primary indirect density effect (PIDE) are calculated using anomalous density derived from the UNB_Topo density model. Two geoidal heights are computed: one including both the DDE and PIDE, and one without including them. The effects of DDE and PIDE on geoidal height range from 10 cm to 27.6 cm. The geoidal heights obtained are evaluated by comparing them with the GSVS17 data. Both models agree with the GSVS17 height anomalies, after which the models behave similarly to each other, but differ from the leveling by 2 cm to 33 cm. The results indicate that the effect of topographical density is significant and must be accounted for when determining the Geoid.