An algorithm for evaluating Wigner coefficients of U(4) ⊃ SUS(2) ⊗ SUT(2) for the coupling [n14,n24,n34,n44]⊗[1,0,0,0]↓[m14,m24,m34,m44] with arbitrary irreducible representation [n14,n24,n34,n44] is provided based on the algebraic expressions with the expansion coefficients of the U(4)⊃SUS(2)⊗SUT(2) states in terms of those in the U(4) canonical basis and the Clebsch-Gordan coefficients of U(4) in the canonical basis. The state expansion coefficients are evaluated as the components of the null space vectors of the spin-isospin projection matrices. Program summaryProgram Title: Wigner-U4-SU2xSU2CPC Library link to program files:https://doi.org/10.17632/4b7gf2phh2.1Licensing provisions: GPLv3Programming language: Wolfram Mathematica version 9.0.1.0 and beyondNature of problem: Numerical evaluation of the elementary Wigner coefficients (EWCs) of U(4)⊃SUS(2)⊗SUT(2).Solution method: Mathematica implementation of the EWCs of U(4)⊃SUS(2)⊗SUT(2) [1] using the expansion coefficients of U(4)⊃SUS(2)⊗SUT(2) states in terms of those of the U(4) canonical basis determined as components of the null space vectors of the spin-isospin projection matrices [2].