痔になりやすい生活習慣とは?

Fortranとgnuplotについて困っています

こんにちは。

現在、Fortranプログラミングでgnuplotを使用した座標出力プログラムを考えているのですが、いまいちわかりません。

中心(0,0,0),半径1の級の球面上の座標を出力するプログラムを考えています。その際
・z>=0のときだけでよい
・gnuplotで確認できる
・Fortranのプログラムファイルでまとめる。

以上を満たすものを考えています。内容自体は基本的なことだとは思いますが、いまいち勝手がわからず困っています。

拡張子.f90でどのようなファイルを作ればいいのか教えていただければ幸いです。

A 回答 (3件)

#2です。



途中で送信しちゃいました。
最後のサンプルプログラムに、x, y, zを計算する処理が抜けてます。
printする前に、

x=sin(th)*cos(phi)
y=sin(th)*sin(phi)
z=cos(th)

って計算が必要かも。
    • good
    • 0

> ・gnuplotで確認できる



のためには、x, y, zが空白区切りで出力されていればOKです。


座標の表現を、

同志社大学 工学部電気系 応用数学研究室 - 2005年度講義 - 解析学II - 2.23 3 次元空間の極座標
http://gandalf.doshisha.ac.jp/~kon/lectures/2005 …

に倣い、球面上の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

とか?
    • good
    • 0

「中心(0,0,0),半径1の級の球面上の座標」といわれてもそれがどういうものか全くわからないし「gnuplotで確認できる」が何

を意味するのかさっぱり分からない.
    • good
    • 0

お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!

このQ&Aを見た人が検索しているワード


人気Q&Aランキング