// Source: https://www.stata.com/manuals/m-1first.pdf page 12 numeric matrix tanh(numeric matrix u) { numeric matrix eu, emu eu = exp(u) emu = exp(-u) return( (eu-emu):/(eu+emu) ) }