data {
  int n;
  real y[n];
}
parameters {
  real<lower = 0> sigma;
  real mu;
}
model {
  y ~ normal(mu, sigma);
  mu ~ normal(0, 1000);
  sigma ~ uniform(1.0/100.0, 1000);
}