function sc3(X,Z,j)

scatter3(X(j(1),:),X(j(2),:),X(j(3),:))
hold on
scatter3(X(j(1),Z==1),X(j(2),Z==1),X(j(3),Z==1),'r')

xlabel(['Metagene ',int2str(j(1))])
ylabel(['Metagene ',int2str(j(2))])
zlabel(['Metagene ',int2str(j(3))])
hold off