Hi Yujia,

1st: you need to get the susceptibility data (X) from your EPR data (this can be simply achieved by integrating the spectra for each temperature);

2nd: you simulate the X vs T data by the Bleaney-Bowers equation using your favorite programing language (MatLab, python, etc)