The Luk-Ulo Melange Complex (LMC) is characterized by a chaotic assemblage of mixed rocks with a block-in-matrix fabric. The exposed blocks consist of various scattered rock types, trending in an ENE-WSW direction. In the case of Mt. Parang, the origin of the diabase remains uncertain, with ongoing debate as to whether it is associated with in situ volcanic activity or represents an exotic block within the melange deposit. Subsurface data obtained through geophysical investigation can aid in modeling the geometry of intrusive bodies using inverse modeling techniques. In this study, we conducted a gravity survey and performed 3D inverse modeling to investigate the subsurface beneath Karangsambung. A total of 818 gravity data points and 28 rock density measurements were integrated with existing geological data to construct an a priori 3D geological model. To ensure the results align with geological concepts, the 3D inversion utilized a stochastic approach, allowing for the incorporation of multiple geological constraints over fifty million iterative procedures. Ultimately, the inversion successfully reduced the misfit between observed and calculated data from 2.71 to 0.55 mGal. Based on the inverted 3D model, the diabase rock in Mt. Parang is identified as having an intrusive origin. The intrusion model exhibited minimal changes in density, volume, and shape during the inversion process. Additionally, the model suggests the presence of a solidified magma reservoir at a depth of approximately 3 km, potentially related to Dakah volcanism. The inverted model also reveals the block-in-matrix structure of the Luk-Ulo Melange Complex in the northern area.