Integrable full discretization of the multi-component short pulse equation

We propose a new formulation of the multi-component short pulse (MCSP) equation that includes the coupled complex short pulse (CCSP) equation as a reduction.

Using Hirota's bilinear method, we construct its $N$-soliton solutions in Pfaffian form.

We then derive integrable semi-discrete and fully discrete analogues of the MCSP equation admitting Pfaffian $N$-soliton solutions.

The resulting fully discrete system provides a practical self-adaptive moving mesh scheme for numerical simulations.

For the parameter sets considered, numerical simulations demonstrate excellent agreement between the numerical and exact solutions, confirming the robustness and high accuracy of the proposed scheme.