n<-1000

m1<-1000

p<-matrix(0,nrow=n,ncol=n)
for(i in 1:99){
p[i,i+1]<-1
}
p[97,104]<-1
p[100,1000]<-1
for(i in 101:199){
p[i,i+1]<-1
}
p[197,204]<-1
p[200,1000]<-1
for(i in 201:299){
p[i,i+1]<-1
}
p[297,304]<-1
p[300,1000]<-1
for(i in 301:399){
p[i,i+1]<-1
}
p[397,404]<-1
p[400,1000]<-1
for(i in 401:499){
p[i,i+1]<-1
}
p[497,504]<-1
p[500,1000]<-1
for(i in 501:599){
p[i,i+1]<-1
}
p[597,604]<-1
p[600,1000]<-1
for(i in 601:699){
p[i,i+1]<-1
}
p[697,704]<-1
p[700,1000]<-1
for(i in 701:799){
p[i,i+1]<-1
}
p[797,804]<-1
p[800,1000]<-1
for(i in 801:899){
p[i,i+1]<-1
}
p[897,904]<-1
p[900,1000]<-1
for(i in 901:999){
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,7,3)

time<-47210.21