library(jack)
mu <- c(2, 1)
nu <- c(3, 2, 1)
LR <- LRmult(mu, nu, output = "list")
LRcoeffs <- LR$coeff
LRparts <- LR$lambda
LRterms <- lapply(1:length(LRcoeffs), function(i) {
LRcoeffs[i] * SchurPol(3, LRparts[[i]])
})
smu_times_snu <- Reduce(`+`, LRterms)
smu_times_snu == SchurPol(3, mu) * SchurPol(3, nu)
Run the code above in your browser using DataLab