（再質問）デバイ温度の式変形
以前回答していただいたのですが、後日回答を見直してみるとよくわからない点が出てきてしまい、回答していただいた方には申し訳ないのですが再質問させていただきます。また、その時の回答を再掲させていただきます。
この中で（１）,(2)より、dk=(Kb*T/h*vo)dx, k=0～xmax=0～Td/Tを(6.4.29)に代入すると、
積分の前の係数として
(3Vo/(2π)^3)*4π*(KbT/h*vo)^3*Kb を得る。
という部分の (KbT/h*vo)^3*Kb となる部分が分かりませんでした。
（再掲）
添付したファイルで、式6.4.29を間の関係式を使って最も下の式を導けるようなのですが、よく分かりません。特になぜ９Rという変数があらわれてくるのかよく分かりません。
---------------------------------
次の様に導出できます。
文字の制約上条件式を書き換えます。
x=h*vo*kmax/Kb (1)
xmax=Td/T (2)
Td=h*vo*kmax/Kb=(h*vo*(6π)^1/3)/Kb*a (3) とすると、
（１）,(2)より、dk=(Kb*T/h*vo)dx, k=0～xmax=0～Td/Tを(6.4.29)に代入すると、
積分の前の係数として
(3Vo/(2π)^3)*4π*(KbT/h*vo)^3*Kb を得る。
(3)より、Kb/(h*vo)=((6π^2)^1/3)/(a*Td) として、係数を整理すると。
3Vo(T/a*Td)^3*Kb となる。aを()外に出して整理すると
9(Vo/a^3)*Kb*(T/Td)^3となり、Vo/a^3=N N:アボガドロ数であるから
9N*Kb(T/Td)^3=9R(T/Td)^3
aは調和振動子の一方向への振幅、a^3は１振動子の占める体積。
Voを１モルの調和振動子系の占める体積と考える。
-------------------------
