C言語　球の体積を求めたい

C言語初心者です。
スクラッチで球の体積を求めるプログラムをつくり、
それに沿ったプログラムを書きたいのですがうまくいきません。

コードは以下の通りです。（スクラッチでは添付画像のように作り、動きました）

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
       
int main(int argc,char *argv[])
{
      double r;//半径
      double V;//体積
      double X;//分割数
      double P;//床面の数
      double x;//sqの中心
      double y;//sqの中心
      double No_x;//番号x
      double No_y;//番号y
      double h;//高さ
      double c;//中心
      double D;//x方向の分割数
      r=1;
      V=0;
      X=1000000;
     
      P=r/X;
      c=P/2;
      D=X;
　　　No_x=No_y=0;
 　　 
　　　while(No_x<X){
           x=No_x*P+c;
           y=No_y*P+c;
           h=(sqrt(r*r)-sqrt((x*x)+(y*y)));
           V+=(V+((P*P)*h));
           No_x=No_x+1;
         }
         
      while(No_y<X){
          No_x=0;
          D=(sqrt(r*r)-((No_y*P)*(No_y*P)));
          No_y=No_y+1;
          }
      
      
      V=V*8;
      
      printf("球の体積は%0.10fです。\n",V);
      
      return 0;
      
      }

【表示される警告】
variable 'D'set but not usedと表示されDが機能していないようですが、
その原因も分かっておりません。
結果は”球の体積はinfです”と表示されてしまっています。

【開発環境】
Linuxのubuntu上で開発しています。

while文の使い方から間違えているのかと思い、
ループ位置の変更など試してみましたがうまくいかない状況です。

超初心者質問で申し訳ありませんが、教えていただけますと幸いです。

cateye

2020/05/09 06:41 編集

行番号を削除して<code>で出てくる ```と```の間に貼り付けて下さい。検証が出来ません。

cateye

2020/05/09 06:52

No_x = No_y = 0;とwhile(No_x < X)の前に全角空白が有ります。修正しましょう。

jimbe

2020/05/09 07:01 編集

> variable 'D'set but not used 変数 'D' が代入される(式の左辺に出てくる)だけで使われていない(右辺に出てくる式が無い)という意味ですが, Dは何のために求めているのでしょうか?

行動規範の内容に同意します

回答2件

ベストアンサー

それなりにテコ入れが必要でした。

元々の質問にあったループについては、xとyで二重ループにする必要があります。
yで求めたDは分割数ではなく、yの移動により短くなっていくxに対応する半径だと思います。
これに伴い、xはDを半径としてXで分割してPx、cxで計算するようにしました。
sqrt(rr)-sqrt(... ではなくsqrt(rr-...)と思います。
V+=V+と二重に足しこんでいました。
分割数は10回にしています。

C
1#include <stdio.h>
2#include <stdlib.h>
3#include <math.h>
4
5int main(int argc, char* argv[])
6{
7	double r;//半径
8	double V;//体積
9	double X;//分割数
10	double P;//床面の数
11	double Px;//床面の数
12	double x;//sqの中心
13	double y;//sqの中心
14	double No_x;//番号x
15	double No_y;//番号y
16	double h;//高さ
17	double c;//中心
18	double cx;//中心
19	double D;//x方向の分割数
20	r = 1;
21	V = 0;
22	X = 10;
23	P = r / X;
24	c = P / 2;
25	D = X;
26	No_x = No_y = 0;
27	while (No_y < X) {
28		No_x = 0;
29		y = No_y * P + c;
30		D = sqrt(r * r - y * y);
31		Px = D / X;
32		cx = Px / 2;
33		while (No_x < X) {
34			x = No_x * Px + cx;
35			h = sqrt(D * D - x * x);
36			V += Px * P * h;
37			No_x = No_x + 1;
38		}
39		No_y = No_y + 1;
40	}
41	V = V * 8;
42	printf("球の体積は%0.10fです。\n", V);
43	return 0;
44}