function visualize(X, y)
    %% Binary classificatio problem
    y = logical(y);
    %% Compute pairwise Euclidean Distances
    D = pdist(X, 'euclidean');
    %% Transform using 2D MDS
    [Yfull, eigvals] = cmdscale(D);
    retained = sum(eigvals(1:2))/sum(eigvals);
    fprintf('%5.2f percent information retained\n', retained*100);

    Y = Yfull(:, 1:2);

    %% Scatter Plot
    scatter(Y(~y, 1), Y(~y, 2), 100, 'black', 'o', 'LineWidth', 2);
    axis(reshape([min(Y)*1.2; max(Y)*1.2], 1, []));
    hold on;

    scatter(Y(y, 1), Y(y, 2), 100, 'black', '*', 'LineWidth', 2);
    hold off;
    l = legend('Control', 'Abnormal', 'FontSize', 15);
    grid;
end
