dmult extension

[John W. Eaton]

On 10-Mar-2008, Jaroslav Hajek wrote:

| please consider the attached patch. It rewrites dmult to allow scaling
| columns of matrix (like A*diag(B)), and also along arbitrary dimension
| (like dmult(ones(3,3,3,3), [1 2 3], 3)).

Will you please edit this to conform more to the Octave coding
conventions (lower case variable names and whitespace) and add
some tests?

Thanks,

jwe