Please have a look at this post: Please read: make it easier to help you . It explains how to ask a question, to make it easier for people to understand and answer.
What is the problem you have with this code? Anyway just trying to execute it, it fails on one of the first lines: pdp = [LinRange(box[1],box[2],ninit),box[3].+1e-4.rand(ninit,1)] because 1e-4.rand(ninit,1) doesn’t make sense. Do you mean 1e-4.*rand(ninit)?