fn initial_j(j1: i32, j2: i32, m1: i32, m2: i32) -> f64
The closed form of the symbol at the lower end j3 = j1 - j2 of a row: Eq. (14) of Xu (1998), treams _initforwardj.
j3 = j1 - j2
_initforwardj