Numerical simulation of axisymmetric spheromak merging