We study the dilute and ultracold unitary Bose gas, characterized by a universal equation of state due to the diverging s-wave scattering length, under a transverse harmonic confinement. From the hydrodynamic equations of superfluids we derive an effective one-dimensional nonpolynomial Schrödinger equation (1D NPSE) for the axial wavefunction which, however, also takes into account the transverse wavefunction. By solving the 1D NPSE we obtain meaningful analytical formulas for the dark (gray and black) solitons of the bosonic system.