At this line: https://github.com/ModelDBRepository/247310/blob/d010593af62583d22d06c4a6b23220ba01d5f0c7/mfm.py#L156 The code used `^` (bitwise xor), but I believe the intention was to use `**`. As a result, `gammasq = 125^2 = 127` instead of `gammasq = 125**2 = 15,625`.