data {
  int<lower=0> n;
  vector[n] y;
  vector[n] x;
}

parameters {
  real beta0;
  real beta1;
  real<lower=0> sigma;
}

model {
  y ~ normal(beta0 + beta1 * x, sigma);
  sigma ~ uniform(1/10.0, 100);
  beta0 ~ normal(0, 100);
  beta1 ~ normal(0, 100);
}