f(1) = ∫[0,∞] Y^(3/2)・e^(-Y) dY
= ∫[0,∞] Z^3・e^(-Z^2) 2Z dZ ; Z = √Y で置換
= 2 ∫[0,∞] Z^4・e^(-Z^2) dZ
部分積分で、
∫[0,∞] Z^4・e^(-Z^2) dZ
= ∫[0,∞] (1/2)Z^3・2Ze^(-Z^2) dZ
= [ (1/2)Z^3・{-e^(-Z^2)} ]_(0,∞) - ∫[0,∞] (3/2)Z^2・{-e^(-Z^2)} dZ
= 0 + (3/2) ∫[0,∞] Z^2・e^(-Z^2) dZ
∫[0,∞] Z^2・e^(-Z^2) dZ
= ∫[0,∞] (1/2)Z・2Ze^(-Z^2) dZ
= [ (1/2)Z・{-e^(-Z^2)} ]_(0,∞) - ∫[0,∞] (1/2)・{-e^(-Z^2)} dZ
= 0 + (1/2) ∫[0,∞] e^(-Z^2) dZ
最後の
∫[0,∞] e^(-Z^2) dZ
= (1/2)√π
は、面倒くさいから Gauβ に頼ろう。(Wolfram じゃなく)
参考↓
http://ja.wikipedia.org/wiki/%E3%82%AC%E3%82%A6%E3%82%B9%E7%A9%8D%E5%88%86
よって、f(1) = (3/4)√π。
補足
すみません。0~∞を抜かしていました。申し訳ありません。