double precision function t(j1,j2,j3) implicit none include 'constants.f' include 'sprods_com.f' integer j1,j2,j3 t=s(j1,j2)+s(j2,j3)+s(j1,j3) end