Bjorn Thanks for the labs feedback.
I tried the below on both Linux64 and Win64 and it worked fine. Could you try again, and if there is still a problem, please examine the error line and try to determine what causes the domain error. Thanks. Chris ── ((4) 11 of 23) Bernoulli numbers ───────────────────────── In BON, the coefficients given in the equations defining the Bernoulli numbers are the alternating binomials, and are used as shown below. The function x: (used in coeffs) invokes the use of extended precision, presenting rationals in the form 1r2 and 1r6, etc.: ) bern=:(] mp numerator % denominator) , ] numerator=:2 with trim on coeffs denominator=:2 with from on - on coeffs coeffs=:abct on inc on x: on # |domain error: on | coeffs=: abct on inc on x:on# |[-3] abct=:(] BC ei) * (_1 with power on ei) on=:@: (this works ok in 32 bit win32 j801) ---------------------------------------------------------------------- For information about J forums see http://www.jsoftware.com/forums.htm