function??para?=?wbl3corrfit（x）
%?f（x）?=?b*a^（-b）*（x-c）^（b-1）*exp（-（（x-c）/a）^b）
%?a?------??尺度參數
%?b?------??形狀參數
%?c?------??位置參數
%?para?=?[abc]

if?~isvector（x）
???error（‘wbl3corrfit函數的輸入必須為向量‘）;
end

if?find（length（x）==size（x））==1
????x?=?x‘;
end

step?=?0.01;
xSort?=?sort（x）;
R_xSort?=?1-（1:length（xSort））/（length（xSort）+1）;
miu=0:s