In this paper,numerical simulations and optimizations for the bilateral-pumping P-doped Raman fiber laser(RFL) are performed using the shooting method based on Newton-Raphson method.