- 124 名前：デフォルトの名無しさん mailto:sage [2006/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))
こうじゃねーか?
間違っていてもシラネｗ
