06/02/07 16:48:19
>>123
まずその前に・・・・・
C=3.0E+8
WL=1.55E-6
これだと単制度で定数が入るので、ここで数値が腐ってしまう。E→Dにせよ。
吸収係数かなんかの計算か？
発想が逆さまだ。もともとの公式を素直にプログラムするほうがよっぽど楽だぞｗ
ほいで
(1 + i) * EXP(iθ) = (1 + i) * (cosθ+ isinθ)
= (cosθ+ isinθ) + (icosθ - sinθ)
= (cosθ - sinθ) + i(cosθ +sinθ)
ゆえに
ZE = EXP(-geff * h) * (1.0d0, 1.0d0) * EXP(CMPLX(0.0d0, th))
こうじゃねーか?
間違っていてもシラネｗ