##TODO: Your unmixing code here
image1 = "mixed1.png";
M1 = Float64.(Gray.(load(image1)))
m,n = size(M1)
image2 = "mixed2.png";
M2 = Float64.(Gray.(load(image2)))

Y = [vec(M1) vec(M2)]';
Wica_k6, Sica_k6, Qica = ica_factorization_alt(Y; obj=absk6)
Wica_k4, Sica_k4, Qica = ica_factorization_alt(Y; obj=absk4)

Sica_k6 = abs.(Sica_k6)
unmixed1 = reshape(Sica_k6[1,:]/maximum(Sica_k6[1, :]), size(M1))
unmixed2 = reshape(Sica_k6[2,:]/maximum(Sica_k6[2, :]), size(M2))
Sica_k4 = 0.0
Sica_k6 = 0.0
gc()