model {
	for (i in 1:N) {
	    y[i] ~ dnorm(mu[i],tau)
	    mu[i] <- alpha + beta*X[i]
	}
#	alpha ~ dnorm(0,0.00001)
#	beta ~ dnorm(0,0.00001)
	alpha ~ dunif(-k,k)
	beta ~ dunif(-k,k)
	lsigma ~ dunif(-k,k)
	sigma <- exp(lsigma)
	tau <- pow(sigma,-2)
}