Pascal constants for diamagnetic susceptibility

Several sources are available for the Pascal constants [4-8]. One of them [4] is used below. 

C The number of terms was increased by 1 because of the DO cycle X2-ZERO DO 20 KK N1,N2 

C Number of all terms for non-negative factorials T1=DMAX1(ZERO,A+AL-C-GA,B+BE-C-GA) 

---