function [muX, SXX] = marginalize_gaussian(mu, Sigma, X, Y, ns)% MARGINALIZE_GAUSSIAN Compute Pr(X) from Pr(X,Y) where X and Y are jointly Gaussian.% [muX, SXX] = marginalize_gaussian(mu, Sigma, X, Y, ns)[muX, muY, SXX, SXY, SYX, SYY] = partition_matrix_vec(mu, Sigma, X, Y, ns);