In this paper, analysis and modeling of bovine babesiosis disease are designed with fractional calculus. The solution for a bovine babesiosis disease and tick populations fractional order system is determined using the Caputo and Atangana-Baleanu-Caputo (ABC) fractional derivatives. Applying the homotopy analysis method and the Laplace transform with polynomial homotopy, the analytical solution of the bovine babesiosis disease has obtained. Furthermore, using an iterative scheme by the Laplace transform, and the Atangana-Baleanu fractional integral, special solutions of the model are obtained. Uniqueness and existence of the solutions by the fixed-point theorem and Picard-Lindel of approach are studied. Numerical simulation has been established for both Caputo and ABC fractional derivative of the proposed system is carried out. The numeric replications for diverse consequences are carried out, and data attained are in good agreement with theoretical outcomes, displaying a vital perception about the use of the set of fractional coupled differential equations to model babesiosis disease and tick populations.