BSXFUN is very useful and fast.
In R2015b we can implement singleton expansion without BSXFUN, see the test suite for details. Have fun :-)
Input and Output:
yf = yourFcn(1) % Output a function handle supported singleton expansion.
yf = yourFcn(0) % Output a function handle unsupported singleton expansion.