x=cbind(c(1,1)) c=x A=cbind(c(.5,.75),c(-.6,1.1)) while(length(c)<100) { x<-A%*%x c<-cbind(c,x) } plot(c[1,],c[2,])