% FEX33  A function used in Example 3.33.
%
function y=fex33(x,om)
n=size(x,1);
y=zeros(n,1);
for k=1:n
  y(k)=gamma(1+x(k))/(x(k)+om);
end
