A QUMOND galactic N-body code, I: Poisson solver and rotation curve fitting