— Loadable Function: y = fprod (x,dim)

Product of elements along dimension dim. If dim is omitted, it defaults to 1 (column-wise products).