大学の課題です。積分を区分求積法やシンプソン則などの方法を用いて計算させろというもので、分割数をコマンドライン引数を用いて3つ読み込ませるというものです。
コンパイルが通り、実行したところで動作を停止しましたというウィンドウが出てきてしまいました。
多分ですがコマンドライン引数の使い方が間違っているのではないかと思っています。
よろしくお願いします。
C言語
1#include<stdio.h> 2#include<math.h> 3#include<stdlib.h> 4 5const double Ans = pow(log(10),2)/2; 6 7double function(double x) 8{ 9 return log(x)/x; 10} 11 12void rect(double a, double b, int n) //区分求積 13{ 14 double h = (b-a)/n; 15 double s = 0.0; 16 17 for(int i=0;i<n;i++){ 18 s += function(a + h*i); 19 } 20 s *= h; 21 22 printf("Rectangle: %g(%g)\n", s,s/Ans); //()内に誤差を表示 23} 24 25void daikei(double a, double b, int n) 26{ 27 double h = (b-a)/n; 28 double s = 0.0; 29 30 for(int i=0;i<n;i++){ 31 s += function(a+h*i)+function(a+h*(i+1)); 32 } 33 s *= 0.5*h; 34 35 printf("Trapezoid: %g(%g)\n", s,s/Ans); 36 37} 38 39void simpson(double a, double b, int n) //シンプソン則 40{ 41 double h = (b-a)/n; 42 double s = 0.0; 43 44 for(int i=0;i<n;i++){ 45 s += function(a+h*i)+4*function(a+h*(i+0.5))+function(a+h*(i+1)); 46 } 47 s *= h/6; 48 49 printf("Simpson: %g(%g)\n", s,s/Ans); 50} 51 52int main(int argc, char* argv[]) 53{ 54 const double a = 1; 55 const double b = 10; 56 57 int n1 = atoi(argv[1]); 58 int n2 = atoi(argv[2]); 59 int n3 = atoi(argv[3]); 60 61 printf("Analytical Solution = %g\n",Ans); 62 63 printf("Division Number = %d\n",n1); 64 rect(a,b,n1); 65 daikei(a,b,n1); 66 simpson(a,b,n1); 67 68 printf("\nDivision Number = %d\n",n2); 69 rect(a,b,n2); 70 daikei(a,b,n2); 71 simpson(a,b,n2); 72 73 printf("\nDivision Number = %d\n",n3); 74 rect(a,b,n3); 75 daikei(a,b,n3); 76 simpson(a,b,n3); 77 78 return 0; 79} 80 81 82 83 84
回答2件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
2017/04/26 03:29