function p = inclexcl(xy)
if length(xy) == 1
  p = xy;
else
  x = xy(1);
  xy = xy(2:end);
  p = x + inclexcl(xy) - x*inclexcl(xy);
end
return