>>244
ちょっと反則臭いやり方だけど、
exp(ix) = cos(x) + isin(x)を使って、
c[n] = a[n] + ib[n]
= ∫[0,π/2]exp(inx)(sin(x))^n dx
= (1/(2i)^n)∫[0,π/2](exp(2ix)-1)^n dx と置くと、

c[n+1] = (1/(2i)^(n+1))∫[0,π/2](exp(2ix)-1)^(n+1) dx
= (1/(2i)^(n+1))∫[0,π/2]exp(2ix)*(exp(2ix)-1)^n dx - c[n]/(2i)

∫[0,π/2]exp(2ix)*(exp(2ix)-1)^n dx = (-2)^(n+1) /{(n+1)*2i} から、
c[n+1] = (i^n)/{2(n+1)} - c[n]/(2i)
a[n+1] = Re(c[n+1])
b[n+1] = Im(c[n+1])

テストでコレを書くと減点されるだろうから
とりあえず答えを知りたいときにだけ使って。