- 締切済み
関数の積分を求めるプログラムで質問です。
シンプソンの公式を用いて積分を求めるプログラムで、 「 y=1/(1+x*x) のように±∞で0に収束するような関数は以下のような無限積分を求めることができる。 ∫(-∞→∞){1/(1+x*x)}dx …(a) ただし、無限区間を分割することはできないので、コンピュータを用いた計算では ∫(-a→a){1/(1+x*x)}dx のaに大きな値を入れることで代用する。この積分を求めるプログラムをシンプソンの式を用いて作成し、以下の3ケースについて値を求めて真値と比較せよ。真値は式(a)を解析的に積分することで求めよ。 ・a=100、N=1000 ・a=100、N=500 ・a=200、N=1000 (Nは以下のプログラムと対応しているものです) 」 というもので参考としてS=∫(1→2){1/x}dxをシンプソンの公式を使い求めるものは以下のものなのですが #include <stdio.h> double f(double x) { return 1/x; } int main() { int i,N; double a,b,h,S; double x[1000],y[1000]; a=1.0; b=2.0; N=10; h=(b-a)/N; for(i=0;i<=N;i++) x[i]=a+i*h; y[0]=f(x[0]); y[N]=f(x[N]); S=h*(y[0]+y[N])/3.0; for(i=1;i<=N-1;i=i+2) { y[i]=f(x[i]); S=S+h*y[i]*4.0/3.0; } for(i=2;i<=N-2;i=i+2) { y[i]=f(x[i]); S=S+h*y[i]*2.0/3.0; } printf("N=%d,S=%15.10lf\n",N,S); } これを生かしてプリグラムをつくりたいのですが、分からなくて困っています。助けてください。
- みんなの回答 （1）
- 専門家の回答
みんなの回答
- επιστημη（@episteme）
- ベストアンサー率46% (546/1184)
関連するQ&A
- 数値積分について
私が作成した下記のJavaの数値解析の矩形法のプログラムで、区間分割数nを10から100まで10ずつ増やして計算値を求め、面積の計算誤差を(計算値－真値(0.68269)）/真値(0.68269)＊100 (%)として計算するプログラムがまったくわかりません。Javaで数値解析をするのは初めてなのでどこが足りないのか、どこを直したら動くのか教えてください public class Kukei { static double f(double x) { // ここに任意の被積分関数を記述 double y = Math.exp(- x * x / 2) / Math.sqrt(2.0 * Math.PI); return y; } public static void main(String[] args) { double a = - 1.0, b = 1.0; // 積分範囲 int n = 10; // 区間分割数 double suti= (n-0.68269)/(0.68269*100); for(int j=0; j<n; j++){ double h = (b - a) / (double)n; // 分割幅 double s = 0.0; n=n; for (int i=0; i < n; i++) { s += f(a + i * h); } s *= h; System.out.println("区間分割数 ＝" + n); System.out.println("矩形法による計算値 ＝" + s); System.out.println("矩形法による計算誤差 ＝" +suti+"\n"); } } }
- 締切済み
- Java
- javaのプログラムについて
質問させていただきます // 数値積分 class Integral1 { public static void main(String[] args) { final int n = 100; // 区間数 final double a = 0; // 始点 final double b = 1; // 終点 final double h = (b-a) / n; // 区間の幅 double x, y; double S = 0; // Sを0で初期化 int i; // 区間 for(i = 1; i <= n ; i++) { // 区間1からnまで x = a + h * i; // 区分積分 y = Math.sqrt(1-x*x); // yを計算 S += y * h; } System.out.println("S = " + S + " 4S = " + 4*S); } } これは区分積分法で計算する数値積分のプログラムなのですが、 これをシンプソン法に改造したプログラムに直していただきたいのですが・・・ どなたかお願いします
- 締切済み
- Java
- シンプソン法の出力結果について
シンプソン法の区間分割数nを10～100まで10ずつ増やして計算値と計算誤差を求めるプログラムを書きに作成したのですが、出力結果に-7.341865999405491E-5などとあります。この「E-数字」とはJavaではどういうこと示しているのですか？ また、計算誤差の求め方は下記のプログラムでいいのですか？ public class Simpson { static double f(double x) { // ここに任意の被積分関数を記述 double y = Math.exp(- x * x / 2) / Math.sqrt(2.0 * Math.PI); return y; } public static void main(String[] args) { double a = - 1.0, b = 1.0; // 積分範囲 int n = 100; // 区間分割数 for(n=10; n <=100; n+=10){ double h = (b - a) / (double)n; // 分割幅 double s, s1 = 0.0, s2 = 0.0; for (int i = 1; i <= n / 2; i++) { s1 += f(a + (2 * i - 1) * h); } for (int i = 1; i <= n / 2 - 1; i++) { s2 += f(a + 2 * i * h); } s = h / 3.0 * (f(a) + 4.0 * s1 + 2.0 * s2 + f(b)); double suti = (s-0.68269)/0.68269*100; //計算誤差＝（計算値ー真値）/真 値×100 System.out.println("区間分割数 ＝" + n); System.out.println("シンプソン法による計算値 ＝" + s); System.out.println("シンプソン法による計算誤差 ＝"+suti+"\n"); } } }
- 締切済み
- Java
- 矩形法プログラムの計算値と計算誤差の求め方
下記の矩形法の計算値を求めるプログラムで、区間分割数nを10から100まで10ずつ増やして計算値を求めるプログラムを作成したいのですが、どのようにしたらいいのか分かりません。Java初心者なのでJavaプログラムが詳しい人は至急ご連絡をお願いします。 *ちなみにこの矩形法プログラムのの計算誤差も分かる人はお願いしたいのですが・・・・・・ 計算誤差の求め方は 計算誤差＝(計算値ー真値)/真値×100（％） (真値=0.68269とする） public class Kukei { static double f(double x) { // ここに任意の被積分関数を記述 double y = Math.exp(- x * x / 2) / Math.sqrt(2.0 * Math.PI); return y; } public static void main(String[] args) { double a = - 1.0, b = 1.0; // 積分範囲 int n = 100; // 区間分割数 double h = (b - a) / (double)n; // 分割幅 double s = 0.0; for (int i=0; i < n; i++) { s += f(a + i * h); } s *= h; System.out.println("区間分割数 ＝" + n); System.out.println("矩形法による計算値 ＝" + s); } }
- ベストアンサー
- Java
- γ関数のプログラム（初心者です）
以下のようにγ関数のプログラムを組みました。 とりあえず整数値を入力すれば、正しい値は返しているということがprintfの４で確認できました。 もとはfortranで組んだプログラムをCに置き換えました。 ですが、実際走らせてみると、４で値は確認できますがsegmentation faultが出てしまいます。 ですからサブルーチンファイル（ユーザー関数？）として利用できません。 何がいけないのでしょうか？ 正しくｙが帰ってくるようにどうなおしたらよいのか教えてください。 #include <stdio.h> #include <math.h> double gamma(double x) { double c[8],y,a,r,b,s; int i; a=1.; r=1.; c[1]=5.771916e-01; c[2]=9.882058e-01; c[3]=8.970569e-01; c[4]=9.182068e-01; c[5]=7.567040e-01; c[6]=4.821993e-01; c[7]=1.935278e-01; c[8]=3.586834e-02; printf("0 %f\n",x); while(1){ if(x>2.){ x=x-1.; a=a*x; printf("1 %f %f\n",x,a); } else if(x<1.){ a=a/x; x=x+1.; printf("2 %f %f\n",x,a); } else{ break; } } x=x-1.; for(i=1;i<8;i++){ b=(double)(i); s=(c[i]*((double)(pow(-1,b))) *((double)(pow(x,b)))); printf("3 %d %f\n",i,c[i]); r=r+s; } y=a*(r+(0.03586834*((double)(pow(-1,8)))*((double)(pow(x,8))))); printf("4 %f\n",y); return y; } main() { double x,y; printf("数字を入力してください。"); scanf("%lf",&x); printf("メインプログラム %lf\n",x); y=gamma(x); printf("%f\n",y); }
- ベストアンサー
- C・C++・C#
- 数値解析の矩形法について
区間分割数nを10から100まで10ずつ増やして計算値を求め、面積の計算誤差が区間分割数によってどのように変化するかを求める矩形法プログラムを下記に作成したのですが、 計算誤差を求める式で ・計算誤差＝（計算値ー真値）/真値×100 （％） （真値は0.68269） と計算するのですが、プログラムでは分母の真値×100 （％）を（）を付ける場合の計算の答えと（）を付けない計算の答えとが全然違います。どうしてこのようなことが起こるのですか？ また、この計算誤差の求め方は（）を付ける場合とつけない場合のどちらが正しいのですか？ public class Kukei { static double f(double x) { // ここに任意の被積分関数を記述 double y = Math.exp(- x * x / 2) / Math.sqrt(2.0 * Math.PI); return y; } public static void main(String[] args) { double a = - 1.0, b = 1.0; // 積分範囲 int n = 10; // 区間分割数 double suti= ((n+10)-0.68269)/0.68269*100; //真値 for(int k=1; k <= 10; k++){ double h = (b - a) / (double)n; // 分割幅 double s = 0.0; n=k*10; for (int i=0; i < n; i++) { s += f(a + i * h); } s *= h; suti= (s-0.68269)/0.68269*100; System.out.println("区間分割数 ＝" + n); System.out.println("矩形法による計算値 ＝" + s); System.out.println("矩形法による計算誤差 ＝" +suti+"\n"); } } }
- 締切済み
- Java
- 定積分のプログラミングについて
刻み幅ｄｘを入力し、ｆ（ｘ）＝４/（１+ｘ＾２）を0～1まで定積分するというプログラムを作ったんですが、なぜか、刻み幅を入力してＥＮＴＥＲを押すと、 ずら～～～～～っと無限に結果の文字列が出てきてしまいます。 書いたプログラムを載せるので、間違っている箇所を教えてもらえませんか。お願いします。 #include <stdio.h> #include <math.h> double f(double m); main(){ double x,F,dx; printf("Numerical integration of y=4/(1+x^2).\n"); /*刻み幅の入力*/ scanf("%lf,&dx:"); /*aとbと刻み幅*/ printf("The lowwer bound:0.0\n"); printf("The upper bound:1.0\n"); printf("The meth width:dx\n"); /*関数*/ for(x=0.0;x<1.0;x+=dx) { F=f(x); /*計算結果を出力*/ printf(" f(%lf)= 4/(1+x^2) =%lf\n",x,F); } } double f(double m) { double y; y=4/(1+m*m); return y; }
- ベストアンサー
- その他（インターネット・Webサービス）
- ６重積分の発散のオーダーの問題です。
S_n =∫_{-∞→∞}dx1∫_{-∞→∞}dx2∫_{-∞→∞}dx3∫_{-∞→∞}dy1∫_{-∞→∞}dy2∫_{-∞→∞}dy3 H(x,y) ただし x とyは３次元ベクトルで x=(x1, x2, x3), y=(y1, y2, y3), H(x,y)=Σ_{i=1, 2, 3}H_i(x,y). H_i(x, y)={I_[a, n](|x|) I_[a, n](|y|)}/{f(x)f(y)}(xi/E(x)^2+yi/E(y)^2)(1/E(x)+1/E(y)){1/G(x,y)}, ただし|x|=√(x1^2+x2^2+x3^2), I_[a, n](|x|)は定義関数で、 a≦|x|≦nのとき I_[a, n](|x|)=1, それ以外のとき I_[a, n](|x|)=0となる関数である。(a>0, n>0) また、f(x)=√(|x|^2+c^2), c>0, E(x)=|x|^2/(2m)+f(x), m>0, G(x,y)={|x|^2+|y|^2+2(x1y1+x2y2+x3y3)}/(2m)+f(x)+f(y)であるとする。 n→∞のときのS_nの発散のオーダーを求めてください。よろしくお願いします。
- ベストアンサー
- 数学・算数
- シンプソン則について
大学の課題で ∫0～5 e^(-3t) * cos^2 *(2t) をシンプソン則を用いて求めるプログラムを作成せよ。積分の刻みは0.01毎 という課題がでました。 私は授業で作ったプログラムを元に、以下のプログラムを作ったのですが、答えが大きすぎて（5より大きいと間違いと言っていました）受け取ってもらえません。どこがいけないのでしょうか？ あと、ｎには どのような数字をいれればよいのでしょうか？ｎは使わないのでしょうか。 お願いします。 #include "stdio.h" #include "math.h" void main() { int n; double i, a, b, h, S, S1, S2, S3; double f (double); a = 0.; b = 5.; printf ("N="); scanf ("%d", &n); h = (b - a) / (double)n; S1 = 0.; S2 = 0.; S3 = f (a) + f (b); /* 両端の値*/ /* 奇数番目の値*/ for (i = 0.01; i < 5; i = i + 0.02) { S1 = S1 + 4. * (f (a + (double)i * h)); } /* 偶数番目の値*/ for (i = 0.0; i < 5; i = i + 0.02) { S2 = S2 + 2. * (f (a + (double)i * h)); } S = (h / 3.) * (S1 + S2 + S3); printf ("%d %lf\n", n, S); } double f (double t) { double y; y = (exp(-3. * t) / 2.) + ((exp(-3.* t) * cos(4. * t)) / 2.); return y; }
- 締切済み
- C・C++・C#
- ３重積分の発散のオーダーを求めてください。
S_n=∫_{-1→1}dz∫_{a→n}dx∫_{a→n}h(x,y,z)dy, ただしa>0, nは実数、 h(x,y,z)=[x^2y^2/{f(x)f(y)}](x^2/{E(x)}^4+y^2/{E(y)}^4)[1/L(x,y,z)], f(x)=√(x^2+c^2), c>0, E(x)=x^2/(2m)+f(x), m>0, L(x,y,z)=(x^2+y^2+2xyz)/(2m)+f(x)+f(y)とします。 この積分のn→∞としたときの発散のオーダーを求めてください。 よろしくお願いします。
- ベストアンサー
- 数学・算数
- ポケトーク同時通訳を契約したが、文字入力による翻訳結果が出ない問題について質問。
- ソースネクストのマニュアルを熟読しても翻訳結果が出ないため、設定に問題があるのか疑問。
- Wi-Fi通信を使用しているが、翻訳ができない原因と解決法を教えて欲しい。
