functions in fermii.i -

             fdi12(x, b)  
 
    return incomplete Fermi-Dirac integral of order 1/2,  
       fdi12(x, b) = integral[b to inf]{ dt * t^0.5 / (exp(t-x)+1) }  
    default accuracy to about 1e-10  

             fdi32(x, b)  
 
    return incomplete Fermi-Dirac integral of order 1/2,  
       fdi32(x, b) = integral[b to inf]{ dt * t^1.5 / (exp(t-x)+1) }  
    default accuracy to about 1e-10  

             fdi52(x, b)  
 
    return incomplete Fermi-Dirac integral of order 1/2,  
       fdi32(x, b) = integral[b to inf]{ dt * t^2.5 / (exp(t-x)+1) }  
    default accuracy to about 1e-10  

             fdim12(x, b)  
 
    return incomplete Fermi-Dirac integral of order -1/2,  
       fdim12(x, b) = integral[b to inf]{ dt * t^-0.5 / (exp(t-x)+1) }  
    default accuracy to about 1e-10  

             #include "fermii.i"  
 
  Incomplete Fermi-Dirac integrals of orders -1/2, 1/2, 3/2, 5/2  
    
  after M. Goano, Algorithm 745, ACM TOMS 21, #3, pp 221-232 (Sept. 95)