经检查是公式写错了，负号位置错了。
k=0.390803;mu=0.463506;sigma=0.219108;
PDF=@(x)(1/sigma).*exp(-(1+k*(x-mu)/sigma).^(-1/k)).*(1+k*(x-mu)/sigma).^(-1-(1/k));
S=quadl(PDF,0.3,2.9)
运行结果
S = 0.897133772014529
