function e = f(x, y)
  e = x^2 + exp(-1/(y+1));
endfunction

% f([2,3], [1,4]) gives and error, but
arrayfun(@f, [2, 3], [1,4])
% works