Here's a more efficient version: gcdn=: _ q:inv i.@>:&.(p:inv) x:@<[EMAIL PROTECTED] ] G=:gcdn%((]+:@x:@[EMAIL PROTECTED]>:*.-:@[<])(-.]-.]&.(p:inv))@(1 _1+2^2<[EMAIL PROTECTED]))
ts=:6!:2,7!:2 ts 'F 1e4' 0.35108 6.83871e7 ts 'G 1e4' 0.136732 1.17188e7 (F-:G)1e4 1 #":G 1e4 4345 -- Raul ---------------------------------------------------------------------- For information about J forums see http://www.jsoftware.com/forums.htm
