> Hi David, No, you didn't make an error - the error is in the equation, for the case when epsilon_surface isn't 1. A factor 3 in the denominator should be replaced with (2*epsilon+1). Then the correction should be ok for arbitrary epsilon_surface. Do you agree? Cheers, Erik