The majority of massive stars are located in binary or multiple star systems. This poses additional challenges to quantitative analyses that use model atmospheres because little information is currently available on the chemical composition of such systems. The members of the quadruple star system HD 37061, which excites the H ii region Messier 43 in Orion, are fully characterised. Accurate and precise abundances for all elements with lines traceable in the optical spectrum are derived for the first time. A hybrid non-local thermodynamic equilibrium (non-LTE) approach, using line-blanketed hydrostatic model atmospheres computed with the Atlas12 code in combination with non-LTE line-formation calculations with Detail and Surface was employed. A high-resolution composite spectrum was analysed for the atmospheric parameters and elemental abundances of the individual stars. Fundamental stellar parameters were derived based on stellar evolution tracks, and the interstellar reddening was characterised. We determined the fundamental parameters and chemical abundances for three stars in the HD\,37061 system. The fourth and faintest star in the system shows no distinct spectral features, as a result of its fast rotation. However, this star has noticeable effects on the continuum. The derived element abundances and determined ages of the individual stars are consistent with each other, and the abundances coincide with the cosmic abundance standard. We find an excellent agreement between our spectroscopic distance and the Gaia Data Release 3 parallax distance.