n<-30

m1<-1000

p<-matrix(0,nrow=n,ncol=n)
for(i in 1:9){
p[i,i+1]<-1
}
p[4,15]<-1
p[7,28]<-1
p[9,30]<-1
for(i in 10:19){
p[i,i+1]<-1
}
p[12,25]<-1
p[17,29]<-1
p[19,30]<-1
for(i in 20:29){
p[i,i+1]<-1
}


duration<-matrix(0,nrow=m1,ncol=n)
for(i in 1:m1){
duration[i,]<-rnorm(n,50,10)
}

set.seed(1)
actual.duration<-rnorm(n,50,10)+rnorm(n,10,3)

time<-1524.737