Pade approximants are used to derive the Gamow-Siegert wave function by back rotation of its complex rotated forms. This provides an unambiguous way to normalize a resonance wave function. 2 X (r)= 1l n m. r sin L The complex rotated resonance wave function is written as P", (r exp(i8))= g C"(8)y"(r) .(2) Stability of the energy versus 8 guarantees that the righthand side of Eq. (2) is a function of r exp(i8) [13,19]. The Gamow-Siegert definition [1, 2] of a resonance allows for the calculation of a quantized complex energy and a complex wave function: The Schrodinger equation is solved with outgoing boundary conditions in all channels open for dissociation of the system. Since this leads to complex wave numbers, a component P(r) of the wave function behaving as A exp[ikr] with k =ko ik & -(ko k]positive) becomes divergent as r, the dissociative coordinate, becomes infinite. Much attention has been paid to this circumstance [3 -10] with the aims of (1) offering a normalization method for the wave function (essentially, how to choose the amplitude A ) and (2) allowing the use of Gamow-Siegert wave functions as a basis set, just as is done for normalizable bound-state wave functions. Most of these problems are bypassed by the now very popular complex rotation method [11,12]: Changing r into p=r exp[i8] with 8)tan '(k&lko) leads to a localized wave function that can be normalized to unity [with P (r), however, instead of~P(r)~in the normalization integral [13]]. Since the complex rotated wave function is defined as the analytically continued Gamow-Siegert resonance wave function upon changing r into a complex coordinate, the reverse problem has also been given some attention. Changing p into p exp( i 8) leads back to the original coordinate r. This operation has been called back rotation [14,15]. One use of it is to ease the calculation of matrix elements in basis set approaches to resonance determination [16). When operating on a complex rotated wave function, back rotation should produce the resonance wave function, with the advantage that a unique function is obtained in this way. Several attempts to obtain this function have been described in the literature [15,17,18]. The simplest procedure for a singlechannel case consists in using a basis of integrable functions, such as free-particle wave functions in a box of length L:Back rotation is then obtained by writing P", (r) = g C"(8)y"(rexp( i 8-)) .( 3) However, the sine functions of Eq. (1) become combinations of hyperbolic functions with the result that it is difficult for large r to get satisfactory interference between these very oscillatory basis functions. The use of other basis sets meets with similar difficulties [16,17].We describe in this Brief Report a different approach which is made possible by the accurate determination of complex rotated resonance wave functions through numerical solution of coupled-channel differential equations with matching of functions obeying shortand long-range boundary conditions. This is based on the Fox-Goodwin procedure [...