楕円とシンプソン公式から円周率を求める

Question

#include<math.h>

double f(double,double,double);

int main(void){

int i;
double a,b,s,x0,dx,y0,y1,y2,xi,s1;
x0=0;
dx=0.1;
s=0.0;
a=2.0;
b=3.0;
xi=x0;

for(i=1;i<=10;i++){
y0=f(xi,a,b);
xi=x0+dx*i;
y1=f(xi,a,b);
xi=x0+2*dx*i;
y2=f(xi,a,b);
s=s+(y0+4*y1+y2)*dx/3;
}

s1=4*s;
printf("楕円の面積→%f \n",s1);

return 0;
}

double f(double x,double a,double b){
double y;
a=2.0;
b=3.0;
y=b/a*sqrt(a*a-x*x);

return y;
}

```楕円ｘ＾２/9+ｙ＾２/4＝１を用いてシンプソンの公式を用いてπの値を求める問題なんですが刻み０.1として
どうしたらいいか全然わからないので教えて欲しいです
シンプソンから楕円をどうやったら使えるのかを知りたいです

Comments

[tanakacge]

書き直しました

Answer 1

何人かから指摘されているように、わかっていないことを丸投げした質問に思われます。

あなたがわかっていないのは「シンプソンの公式」ですか？　「コードの書き方」ですか？
コンパイルや実行時にはどのようなエラーが出ているのでしょうか？

質問する際には、最低限「何がわかっていないか」と「どのような状況か」を明示すると回答者にも親切な質問になります。
（あとコードはバッククオートで囲んだ方が見やすいですよ）

おせっかいはこのぐらいにして、プログラムが動かない原因はfor文の後のカッコ({)が抜けているからです。
ただそれを差し置いてもπを求めるプログラムにはなっていないと思うので、まずは公式の理解から初めてみてはいかがでしょうか？

個人的には以下のサイトなどわかりやすいように感じました
http://www.aoni.waseda.jp/ichiji/2009/first-term/ruby-12-1.html

Comments

[tanakacge]

サイト先はc言語で表すとどうなっているのですか？