"MarkovKette" <-
function (uemat, N=4, start=1) 
{
# Simuliert eine Markovkette mit der Übergangsmatrix uemat
# für die Zustände 0,1,...,ncol(uemat)-1
# Startwert:
i<- start
mk<-i-1
for(n in 1:N){
  u<-runif(1)
  for(l in 1:ncol(uemat)){
    if(i == l){
       psum<-cumsum(uemat[,l])
       for(k in 1:nrow(uemat)){
           if(k==1){ 
                if(u<=psum[1])
                {j<-1}
           }
           else{
                if(psum[k-1]<u & u<=psum[k])
                {j<-k}
           }
       }
    }       
  }
  i<-j
  mk<-c(mk,i-1)
}
list(Markovkette=mk)
}