dyn.load("dynamic.so")
is.loaded('gibbs')

dynamic<-function(n=10000,thin=500,seed=trunc(runif(1)*1e6))
{
  tmp=.C("gibbs",as.integer(n),as.integer(thin),
         as.integer(seed),x=as.double(1:n),
         y=as.double(1:n))
  mat=cbind(1:n,tmp$x,tmp$y) 
  colnames(mat)=c("Iter","x","y")
  mat
}

results<-dynamic()

dim(results)

plot(results[,2],type="l");
plot(results[,3],type="l")
plot(results[,2:3],pch=18,cex=0.6)