Fortranで直交座標から極座標変換のプログラム

Fortranで直交座標から極座標変換のプログラム

FDTD法を用いて、散乱電場を求める際、最初Ex(i,j,k), Ey(i,j,k), Ez(i,j,k)を求めましたが、
それから座標をr方向に座標変換したく、プログラムを作ろうと思っているのですが、どのように書いてよいのか悩んでいます。
単位ベクトル　r = (x,y,z)=(sinθcosφ,sinθsinφ,cosθ)と定義できるのですが、これを
どのように極座標のプログラムとして書いてよいのかわかりません。
どなたかわかる方がいらっしゃたら教えて下さい。よろしくお願いします。

No.1 ベストアンサー 回答者： spring135 回答日時： 2010/07/27 12:29

座標変換（デカルト座標から極座標）に伴う単位ベクトルの変換またはベクトル成分の変換を行おうということなら下記URL参照。

参考URL：http://www-d.ige.solan.chubu.ac.jp/goto/docs/mat …