matlab中一组数取百分位数_一段计算分位数的程序，哪位帮忙看看。

function qth=quantle(p)

A=-5;

M2=100;

qth=0;

for i=1:length(k)

for j=1:M2

B=A+j*(0-A)/M2;

if simprl(gst(B),A,B)-P<=0.00001

qth(i)=B;

break

end

end

qth;

function s=simprl(f,a,b)

M1=10000

h=(b-1)/(2*M1);

s1=0;

s2=0;

for k=1:M1

x=a+h*(2*k-1);

s1=s1+feval(f,x);

end

for k=1:(M1-1)

x=a+h*2*k;

s2=s2+feval(f,x);

end

s=h.*(feval(f,a)+feval(f,b)+4.*s1+2.*s2)./3;

function f=gst(x)

global k n lambda sigma

slambda=(1+3.*lambda.^2-4.*lambda.^2.*beta(2./k,(n-1)./k).^2.*beta(3./k,(n-2)./k).^(-1)).^1/2;

c=0.5.*k.*beta(1./k,n./k).^(-3/2).*beta(3./k,(n-2)./k).^1/2.*slambda