A boundary correction method is proposed for the numerical solution of the degenerate four-wave mixing in photorefractive crystals. The present method can easily deal with the problem of split boundary conditions. Results from the present method shows that the incidence angles of and ratio of pump beam intensity in the PCR lead to a significant improvement; The line absorption of the interacting beams by the crystals is nonlinear in the PCR.