あ、もっと簡単になる
e^((1+t^2)/2) = 1/t
より
t = √W(1/e) (WはランベルトのW関数)
S = 8π(1+2W(1/e))/(1+W(1/e))^2