model {for (t in 2:500) {
y[t] ~ dnorm(0,P[t])
h[t] <- gam+alph*pow(y[t-1],2)
P[t] <- 1/h[t]}
gam ~ dgamma(1,1)
alph ~ dunif(0,1)}
