x = rand (100, 100); bm1 = bm2 = 0; for i = 1:30 tic; y1 = poissinv (x, i); bm1 += toc; tic; y2 = poissinv3 (x, i); bm2 += toc; if (! isequal (y1, y2)) printf ("Failed at index %d\n", i); keyboard; endif endfor