Расчет крутящего момента цилиндрической магнитной муфты

[saberhagen]

Здравствуйте товарищи!
Пытаюсь выполнить расчет крутящего момента цилиндрической магнитной муфты по методике А.Я. Красильникова,
опубликованной в журнале "Вестник машиностроения 2009 №6", не совпадают результаты.
Прикладываю исх. код MATLAB
https://yadi.sk/i/lxVSj4tntq3eh
и статью
https://yadi.sk/i/7kptEL9Mtq3em
Разумеется, есть вероятность, что не вижу опечаток (хотя несколько раз перепроверил).