Numerical study of M=1 resistive internal kink mode in a cylindrical Tokamak