wolffd@0: function p=mysubset(small,large)
wolffd@0: % MYSUBSET Is the small set of +ve integers a subset of the large set?
wolffd@0: % p = mysubset(small, large)
wolffd@0: 
wolffd@0: % Surprisingly, this is not built-in.
wolffd@0: 
wolffd@0: if isempty(small)
wolffd@0:   p = 1; % isempty(large);
wolffd@0: else
wolffd@0:   p = length(myintersect(small,large)) == length(small);
wolffd@0: end