Jiles–Atherton parameters

How can I calculate Jiles-Atherton parameters. What will be the equations and how to write them in MATLAB code. can anyone give me code or suggestions?

[Rik]

I don't think this is a Matlab problem yet. Would you be able to explain with pen and paper how to do it? You need to explain the actual Matlab problem, otherwise you reduce the number of people who can help you to only those who understand this specific problem as well.

Have a read here and here. It will greatly improve your chances of getting an answer.