Polynomial interpolation solutions for the non-linear multi-dimensional integral equation for the evolution of plasma wave amplitude