Calculate coshm(x) = cosh(x) - 1
elemental real function coshm_fun(x) real, intent(in) :: x coshm_fun = cosh(x) - 1 end function coshm_fun