Hi Labog,
The definition itself seems correct. I don't know the range of H (which is, I guess, an independent variable, not dependent, and M is), but for example, when H=1e-50, your function returns 8.7431693989071E42 .
Sinh may easily overflow such that:sinh(700)=5.071160273675E303
sinh(750)=--
using tanh might help to avoid overflowing.
Hope this helps.
--Hideo Fujii
OriginLab