series expand[(sqrt(4*k*T+1)-1)/(2*k)]