Mediante simulaciones de Monte Carlo, se estudia el comportamiento magnético de un sistema de Ising ferromagnético, sobre una red cuadrada compuesta por dos subredes A y B, con espines SiA = 3/2, y s jB = 5/2, y campos cristalinos D A y D B , respectivamente. La interacción de intercambio J1 , entre espines SiA y s jB , es ferromagnética. Se calcula el diagrama de estados base, con diversos puntos donde coexisten tres fases. Se analiza el efecto de los campos sobre las magnetizaciones de las subredes, el calor específico y la temperatura crítica. Se hallan posibles transiciones de primer orden para D A = D B , D A< 0 y D B< 0. Cuando alguno de los campos es nulo y el otro se incrementa, la temperatura crítica se incrementa hasta un valor límite. El efecto de D B> 0 es más determinante en el ordenamiento del sistema, que el efecto de D A> 0.