Hi Ted, that is correct. If you input an NxP predictor matrix, you will get a beta vector that is P+1 in length, the extra term (the first term) is the beta0 term.
Note the default behavior is to exclude the column vector of 1s from the predictor matrix. glmfit adds that.