Watch this Author's files
Worked great as an example of using nlinfit. Made it work for a truncated Poisson. Many thanks.
Good stuff. As usual, be prepared to edit C-code to reflect locations of standard headers in your operating system, and to comment out data types not supported for your operating system.
Reliable and robust, but it takes 6.6 seconds for two 10,000 point-functions with ~20 intersections, with the robust flag off. May be I can figure out how to make it less robust, but a bit faster.