The numerical methods for finding the optimal parameters in 640 Gbps (16 channels $ imes$ 40 Gbps) WDM system with optical phase conjugator (OPC) are proposed, which effectively compensate the distorted overall WDM channels. The considered optimal parameters are the OPC position and the dispersion coefficient of fibers. The numerical approaches are accomplished through two different procedures. One of these procedures is that the optimal OPC position is previously searched and then the optimal dispersion coefficient is searched at the obtained optimal OPC position. The other is the reverse of the above procedure. From the numerical results, it is confirmed that two optimal parameters depend on each other, but less related with the searching procedure. The methods proposed in this research will be expected to alternate with the method of making a symmetrical distribution of power and local dispersion in real optical link which is a serious problem of applying the OPC into multi-channels WDM system.