This paper presents an efficient numerical method for calculating spatial power profiles of both signal and pump with significant Interchannel Stimulated Raman Scattering (ISRS) and backward Raman amplification in multiband systems. This method was evaluated in the optimization of a C+L+S/C+L+S+E 1000km link, employing three backward Raman pumps, by means of a closed-form EGN model (CFM6). The results show a 100x computational speed increase, enabling deep optimization which made it possible to obtain very good overall system performance and flat GSNR.