Posterior distributions for the radial M/L_F814W (left) and M/L_F160W (right) variation profiles obtained after the burn-in phase of the ADAMET MCMC optimization process for the JAM applied to the stellar kinematics of the nucleus of NGC 4736, measured using the NIRSpec G235H/F170LP observations, are presented. In each triangle plot, the analysis involves four parameters (i, M_BH, $M/L_{2\stackrel{″}{.}0}$, and β_r), shown as scatter plots for their 2D projected PDF with white indicates maximum likelihood and a CL at 1σ, while black corresponds to a confidence level at 3σ. The top histograms display the 1D projected distributions of the best-fitting model parameters, including the best-fit values (indicated by vertical black lines, corresponding to the highest likelihood within the PDF) and 3σ uncertainty (indicated by either side black dash lines). Insert plots in the top-right of each triangle plot illustrate the agreement by directly comparing the V_rms data with that generated from the best-fit JAM using the radially varying M/L_F814W profile (left) and the constant M/L_F160W value (right) in the same velocity range. It is evident that the constrained model shows excellent agreement with the data across the NIRSpec FoV.