Hi. I found two cases where the Mathematica code differs from the papers "code". See the attached, the changes I made are highlighted in yellow.
With the changes the T and TP values becomes close to the expected.
Attachments:
