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

