# NOT RUN {
a <- sapply(1:500, function(i) rnorm(n=50, mean=0))
b <- sapply(1:500, function(i) rnorm(n=50, mean=0))
x <- t(cbind(a,b))
l <- c(rep(0,500), rep(1,500))
h <- hpb(D=x, L=l, r=30,t=50)
# }
Run the code above in your browser using DataLab