Posterior distributions HD 25723 system orbital parameters. Each panel contains ~50 000 samples that are tested for 100 kyr dynamical stability using the MVS integrator. Stable samples are plotted in black, while non-stable samples are plotted in red. The upper panels show the probability density distribution of each orbital parameter of the stable samples (black) and of the non-stable samples (red). Contours indicate the 68.3%, 95%, and 99.7% confidence interval levels (i.e., 1σ, 2σ, and 3σ). The blue and green cross indicate the dynamical best-fit solution and the stable solution, respectively.