data {
     int<lower=0> N;
     vector[2] y[N];
}
parameters {
         vector[2] mu;
         cov_matrix[2] sigma;
}
model {
      y ~ multi_normal(mu,sigma);
}