>> arrayfun(@primalityByTrialDivision,(1:14)) ans = 0 1 1 0 1 0 1 0 0 0 1 0 1 0