function S = Psi(K,m,n)
if n==0
    S = pl(K);
else
    S = sum(pl(K)) + Psi(getK1(K,m,n),m,n-1);
end
return