Algorithm for reconstructing complex coefficients of Laguerre–Gaussian modes from the intensity distribution of their coherent superposition