Notes. The epoch is T₀ = 2 456 224.724705. The values below correspond to that of best-fit solutions with χ² = 2578, with individual contributions: , , and . The masses of components in M_⊙ units were , , and . Component B is on its lower limit . The mass of component C is very poorly constrained, it is more of a distant test mass compared to the others. The 3σ uncertainties of the elements were determined by a simplified one-dimensional χ² mapping, assuming a relative increase of χ² by a factor of 1.13, i.e. suitable for the number of degrees of freedom we have (ν = 908). The uncertainty values were verified using the bootstrap method with 100 random selections of datasets and corresponding simplex optimisation, but realistic uncertainties are likely to be larger than that because there are a number of local minima with statistically equivalent χ² values. We do not report a full correlation matrix of our solution here. Its non-diagonal terms indicate higher values of uncertainties for those elements that are strongly correlated or anti-correlated with others (e.g. r_m^Aa,ω2 = 0.74, r_a1,Ω2 = −0.77, r_a1,ω2 = −0.80, r_a3,i3 = −0.79).