function P = skel2p(n,a)
global SKEL
if 0
SKEL{1} = list_all_skel(1,'ab')
SKEL{2} = list_all_skel(2,'ab')
SKEL{3} = list_all_skel(3,'ab')
end

P = zeros(n);
alf =  SKEL{n}(1).ALF;
i = find(alf==a);
for ai = 1:length(SKEL{n})
  d =  SKEL{n}(ai).delta(:,i);
  for q=1:n
    P(q,d(q)) = P(q,d(q)) + 1;
  end
end
%P = pnormdim(P,1);
return