質問: fortranとgnuplotについて困っています fortranとgnuplotについて困っています こんにちは。 現在、fortranプログラミングでgnuplotを使用した座標出力プログラムを考えているのですが、いまいちわかりません。 中心(0,0,0),半径1の級の球面上の座標を出力するプログラムを考えています。その際 ・z>=0のときだけでよい ・gnuplotで確認できる ・fortranのプログラムファイルでまとめる。 以上を満たすものを考えています。内容自体は基本的なことだとは思いますが、いまいち勝手がわからず困っています。 拡張子.f90でどのようなファイルを作ればいいのか教えていただければ幸いです。
回答: > ・gnuplotで確認できる のためには、x, y, zが空白区切りで出力されていればokです。 座標の表現を、 同志社大学 工学部電気系 応用数学研究室 - 2005年度講義 - 解析学ii - 2.23 3 次元空間の極座標 http://gandalf.doshisha.ac.jp/~kon/lectures/2005.calculus-ii/html.dir/node31.html に倣い、球面上の1点って事なら、 x = sinθcosφ y = sinθsinφ z = cosθ で計算され、 > ・z>=0のときだけでよい の条件ならば、 θは0~π/2までループ φは0~2πまでループ で、それぞれの状況でのx, y, zの値を出力します。 fortranはしばらく見ていないのでうろ覚えですが、ざっくり書くと、 double pi, th, phi, x, y, z pi=3.14159265358979 do th=0.0, pi/4, pi/60 do phi=0.0, pi*2, pi/30 print *, x, y, z end do end do とか?