I want to solve equation (1) using the finite difference method employing the initial conditions (eqs. 2-3, 4). The initial condition at z = -2400 nm for every value of r, gives the P(r,z) and will serve as a source for the differential equation that needs to be solved. Any suggestions will be appreciated.