This runs without error:
A = rand(4)
B = randn(5)
bsxfun(@na_div,A,B)
function[out] = na_div(x,y);
if isnan(x) & ~isnan(y)
out = 1/y;
elseif ~isnan(x) & isnan(y)
out = x;
else
out = x/y;
end
end
except for throwing:
Error using bsxfun
Non-singleton dimensions of the two input arrays must match each other.
which is true for bsxfun generally, so not specifically with your function.
It runs without error of both matrices are the same sizes, although different arguments will likely crash it in some situations.
