n<-10

m1<-1000

p<-matrix(0,nrow=n,ncol=n)
p[1,2]<-1
p[2,3]<-1
p[3,5]<-1
p[5,9]<-1
p[9,10]<-1
p[4,6]<-1
p[4,7]<-1
p[7,8]<-1
p[8,10]<-1
p[3,7]<-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<-296.4193