poissinv.m: Overhaul function for accuracy and performance (bug #34363). * poissinv.m: Eliminate while(1) loop which may never terminate. Subsitute for loop with limit of 20. If solution has not converged, use Mike Gile's analytic_approx routine to finish calculating a solution. Introduce special case for scalar lambda which speeds up calculation 6X.