(35S) - Quadratic root finder -- high performance

[Douglas De Boer]

Quote:Case 1:

a = 1E-13 b = -2 c = 1

R1 = 2E13 R2 = 0.5

Jeff, Thanks for pointing out this "Cadillac" program. The first case is solved incorrectly by my program--an outright bug. My program gives R2 = 0. I'll look into it and fix that when I find time.

Most of the other cases show some discrepancies compared to my program. These look like matters of numerical accuracy. Since the "Cadillac" program has been around for a while, I think it is fair to trust that one rather than mine at this time. I'll keep looking at this, mainly since it is an interesting puzzle for me.

Thanks for your comments!

Regards,

--DDB