Hi,
I found an article by Mark C Williams on Simulation and Modelling of Non Linear Magnetics.
There is a code written in Verilog on the Jiles Atherton Model, but I want a simple code to be able to produce a hysteresis curve through MATLAB.  

Is there a way in which the code be redesigned? I tried but there are certain symbols or syntax that I am not familiar with.
I have attached the paper for the same.

Thanks