absk6(y) = abs(kappa(y,6))
function ica_factorization_alt(Y::Matrix, k::Integer=minimum(size(Y)),
        maxiters=1000;obj=absk4)
    m, n = size(Y)
    μy = mean(Y, 2)
    Ymean = μy * ones(1, size(Y, 2))
    Ytil = Y - Ymean

    UsV = svds(Ytil; nsv=k)[1]
    U = UsV[:U]
    s = UsV[:S]
    V = UsV[:V]
    S = Diagonal(s)

    Qica = optimQ(sqrt(size(V, 1))*V',k,maxiters;obj=obj)
    Vica = V*Qica ## TODO: Fill in formula for Vica'
    Wica = U*Diagonal(s)*Qica   ## TODO: Fill in formula (1) for Wica
    Xica = pinv(Wica)*Ymean + Vica'
    return Wica, Xica, Qica
end