Based on the macroscopic description of a type II anisotropic superconducting material, this paper presents numerical methods to model the magnetic induction of these materials in a critical state. In this case, an elliptical model of the critical state is used to describe a type II anisotropic superconductor in terms of parallel geometry. It describes how the mathematical problem is solved in order to find the distribution of magnetic induction in the material. This consists of solving a system of ordinary differential equations. In order to do this, the standard libraries of MATLAB software were used. Two routines were especially used, the first of which gives an initial rough result which is then refined with the second routine. The authors have ample experience in the macroscopic study of the magnetic properties of type II superconductors as well as international publications in this specialized study area.