mysubset <- function(small, large) {
# function p=mysubset(small,large)
# % MYSUBSET Is the small set of +ve integers a subset of the large set?
# % p = mysubset(small, large)

# % Surprisingly, this is not built-in.

# if isempty(small)
#   p = 1; % isempty(large);
# else
#   p = length(myintersect(small,large)) == length(small);
# end
}