“…The spin-1 Blume-Capel model was studied by a variety of methods such as mean-field approximation (MFA) [1,2], effective-field approximation (EFA) [4][5][6][7][8][9], the Bethe lattice approximation [10], series expansion method (SE) [11,12], cluster variation method (CVM) [13], Monte Carlo (MC) simulations [14][15][16][17], renormalization-group (RG) method [18,19] and rigorous inequality correlation function [20][21][22]. Most of these studies were done on hexagonal and rectangular lattices.…”