Problem with bsxfun and custom function

I am trying to apply elementwise a custom 'division' function.

Here is an example

A = rand(4)
B = randn(5)

bsxfun(na_div,A,B)     

with

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
    

I get the following error message: 'Not enough input arguments'.

The function works properly when applied to two numbers. Any idea why? Or any other solution to get the same result?

Thanks

Answer 1

[Star Strider]

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.