式をC言語で立てれません

式をC言語で立てれません

こんにちは。私は今実験で円板振動子の中心軸上の音圧分布を求めています。
そこで写真のような図をプロットするために与えられた式をC言語で書いています。
式の詳細は
P/2ρcV ＝　|sinπ(√(r/λ)^2 + (a/λ)^2)(√はここまで)-(r/λ)^2)|
です。
しかしうまくプロットできません。
プログラムは
-------------------------------------
#include <stdio.h>
#include <math.h>

#define F 120 //周波数(Hz)
#define C 1500 //音速(m/s)
#define A 0.2 // 半径(m)
int main()
{
double x, m, s, p, o, w, l;
int i;

for (i=0; i<100; i++)
{ x= (double)i;　//観測点までの距離を表しました。距離が変化することで音圧が変わるためです。
l = C/F; //λをlとして音速÷周波数で定義できます。
m = x/l;　//x軸です。
o = m*m + (A/l)*(A/l);　//式のsinの中身です。
s = sin(M_PI*(sqrt (o) - m)); //sin全体を定義します。
w = fabs(s);　//sin全体を絶対値に置き換えます。
p = w;　//新しい変数に入れます。
printf("%9.9f %9.9f＼n", m, p);
}
return 0;
}
------------------------------
といった感じです。先生からはx軸のプログラムだけいじれば勝手に
y軸も出てくるとヒントを頂いているのでx軸についてのプログラムを
書いています。コンパイルはできても写真のようにいきません。
半径やλや周波数はa/λ=2.5と書いてあったので値は推測です。
C言語はあまり得意ではないので困っています。
お分かりになる方、ご教授お願い致します。

No.2 ベストアンサー 回答者： nitho_t 回答日時： 2010/05/21 19:21

C/Fはint

なので精度が出ない。
科学技術計算は（メモリが許す限り）精度をとっておくか、（きちんと理解していれば）不要な精度を落とすかどちらか

l=(double)(C)/F;
で問題ないと思います。

この回答へのお礼

返事が遅れてしまい大変申し訳ございません。
double型も関係ありました！ありがとうございました。

お礼日時：2010/05/22 23:07

No.1 回答者： magicalpass 回答日時： 2010/05/21 17:06

あまり良く分かりませんが、とりあえず

> s = sin(M_PI*(sqrt (o) - m)); //sin全体を定義します。
は
　s = sin(M_PI*(sqrt (o) - m*m));
じゃないかと……

この回答への補足

すみません！
式はP/2ρcV ＝　|sinπ(√(r/λ)^2 + (a/λ)^2)(√はここまで)-(r/λ))|です！
最後の二乗いりませんでした。
これでよろしくお願いします。

補足日時：2010/05/21 17:39

この回答へのお礼

プロットできました！回答して頂きありがとうございました。

お礼日時：2010/05/22 23:08