Dear MRMC developers,<br><br><div>I believe I&#39;ve found a bug in mrmc. When using transient analysis, MRMC outputs:<br><br><div>       &gt;&gt;&gt;&gt;ERROR: Fox-Glynn: lambda &gt;= 25, underflow. The results are UNRELIABLE.<br>
<br></div><div>which, I think, is the result of a wrong formula in foxglynn.c, line 199:<br><br>       k_prime = k_rtp * sqrt_2 + 3.0 / 2.0 * sqrt_lambda;<br></div><div><br></div><div>This line should be:</div><div><div><br>
       k_prime = k_rtp * sqrt_2 + 3.0 / (2.0 * sqrt_lambda);<br></div><div><br></div><div>as written in Corollary 3 of the &quot;Computing Poisson Probabilities&quot; paper. </div><div>Without this line, mrmc always thinks that the right Poisson tail will underflow, writing the warning message.</div>
<div><br></div><div><br></div><div>Amparore Elvio</div></div></div>