x <- cbind(1, 3, 2, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 3, 1, 0, 1,0, 0, 0, 2, 0, 0,1, 0, 0, 0, 1, 0,1, 0, 0, 0, 0, 0,1, 1, 1, 0, 1, 0,0, 0, 0, 0, 1, 0,0, 0, 0, 0, 0, 1,0, 0, 1, 0, 0, 0,0, 1, 2, 0, 2, 1,1, 1, 0, 0, 0, 2,1, 1, 1, 1, 1, 0)
x <- matrix(x, ncol=6)
svd <- svd(x)
// löschen der singulärwerte, zB
svd$d[3] <- 0
svd$d[4] <- 0
svd$d[5] <- 0
svd$d[6] <- 0
// neue term2doc matrix in xx speichern:
xx <- svd$u %*% diag(svd$d) %*% t(svd$v)