Computer simulation of hexapole aberration correctors