improved test harness, point evaluation test stengthend

This the start of a change to the recursive Horner
evaluation, we will require that the first dimension
of x should range over the variables of p (rather than
the last).  To start we modify the point-evaluation
test to use a column vector x (and it now fails)
polyval_cube: polynomial variables in first dimension of x
polyval multi-argument form
other functions/tests use polyval variable-in-dim-1
updated maxmodnd to first-variable-in-dim-1
fix for matlab reshape prissiness
version bump