findChildren(parents)

function children = findChildren(parents)

n = numel(parents);
children = cell(1,n);

for i=1:n
    if parents(i)
        children{parents(i)} = [children{parents(i)} i];
    end
end