function data = Gaussian2D(mu, variance, N)
    data = randn(N, 2);
    data(:, 1) = data(:, 1) * sqrt(variance(1)) + mu(1);
    data(:, 2) = data(:, 2) * sqrt(variance(2)) + mu(2);