Various changes

• Yk (closest value) now includes a third output with indices of selected values
• Z^ (Cartesian power) now allows inputs in reverse order if the second is an array of at least two elements
• m (ismember) and u (unique) now accept numeric content in cell array input
• Z+ (2D convolution; see also Y+) now produces output the same size as the input
• Bug corrected in e (reshape) for logical input
• Modular indexing has been added, only for scalar indices (thanks to @dennis for the suggestion)
• Bug corrected in ZA. In addition, now it allows cell array input, and ignores non-recognized digits
• o with cell input now allows specifying fill value and padding side
• Z} (split array) now by default splits along the first non-singleton dimension
• Added W function: 2 raised to input, element-wise
• Updated the Fibonacci example in the documentation (1 byte fewer using y), and included the new features
• YS (circshift) now allows second input to be a vector with third input specifying dimension when first input is a 2D array
• Xd (diag) now includes spdiags functionality
• format function has been integrated with ZG
• h (horizontal concatenation) and v (vertical concatenation) now linearize their inputs if sizes don't match