where, J - Angular Momentum Number g - Lande g-factor, Value 1 or 2 B - Bohr Magneton (9.274X10^-24 J/T) N - Number of spins K - Boltzmann Constant (1.381X10^-23 J/K)
I'm not sure if I understood correctly, but I assume that you want to create a function of the Brillouin Function for curve fitting, you may refer to the following tutorial showing how to create a user defined fitting function: