Accurate parameter estimation is essential for gravitational wave data analysis. In extreme mass-ratio inspiral binary systems, orbital eccentricity is a critical parameter for parameter estimation. However, the current software for the parameter estimation of the gravitational wave often neglects the direct estimation of orbital eccentricity. To fill this gap, we have developed the MatBYIB, a MATLAB-based software (Version 1.0) package for the parameter estimation of the gravitational wave with arbitrary eccentricity. The MatBYIB employs the Analytical Kludge waveform as a computationally efficient signal generator and computes parameter uncertainties via the Fisher Information Matrix and the Markov Chain Monte Carlo. For Bayesian inference, we implement the Metropolis–Hastings algorithm to derive posterior distributions. To guarantee convergence, the Gelman–Rubin convergence criterion (the Potential Scale Reduction Factor $\widehat{R}$) is used to determine sampling adequacy, with MatBYIB dynamically increasing the sample size until $\widehat{R}<1.05$ for all parameters. Our results demonstrate strong agreement between predictions based on the Fisher Information Matrix and full MCMC sampling. This program is user-friendly and allows for the estimation of the gravitational wave parameters with arbitrary eccentricity on standard personal computers.