www.gusucode.com > 地统计学中的克里格插值法matlab源码程序 > devico.m

    function [Cs]=devico(key)
keybar=mean(key);
s=0;
for i=1:length(key)
    s=s+(key(i)-keybar)^3*length(find(key==key(i)));
end
Cs=s/(length(key)*std(key)^3*keybar);