I am stuck with the solution of coupled nonlinear equation for beam propagation inside a non linear crystal in Second Harmonic Process.
these equations can be either solved by symmetric split step beam propagation method or Fourth order Runge - kutta algorithm solver .
can anyone assist me on how to go with these method to solve the equations given in the pdf attached