function v=percentile(x,p) p=p/100; xord=sort(x); n=length(x); if ceil(n*p)==n*p v=1/2*(xord(n*p)+xord(n*p+1)); else v=xord(ceil(n*p)); end