前提・実現したいこと
台形公式を用いて[1,exp(4)]の範囲で7xlog(x)の数値積分を行いたいのですが、logが読み込まれません(?)
お手数ですが助言をいただけると幸いです。よろしくお願いします。
発生している問題・エラーメッセージ
exit status 1 /tmp/cc83RQzh.o: In function `f': main.c:(.text+0x1f5): undefined reference to `log' collect2: error: ld returned 1 exit status
該当のソースコード
C
1#include <stdio.h> 2#include <math.h> 3 4#define nmax 100 5#define pi atan(1.0)*4.0 6#define xa 1.0 // 積分範囲_下限 7#define xb exp(4) // 積分範囲_上弦 8#define h 0.5 // 刻み幅 9 10double f( double ); 11double Sdai( int ni,double y[] ); 12 13int main(void) 14{ 15 double S,Sd,Err_d; 16 double xt[nmax],yt[nmax]; 17 int i,j,l,ni; 18 19 /** きざみ幅 **/ 20 ni=(xb-xa)/h; 21 22 /** 初期値の代入 **/ 23 for(i=0; i<=ni; i++){ 24 xt[i] = xa + i*h; 25 yt[i] = f(xt[i]); 26 printf("%12.9f, %12.9f \n", xt[i], yt[i]); 27 } 28 29 /** 面積の計算 **/ 30 S=15859.7769; // 真値 31 Sd=Sdai(ni,yt); 32 Err_d=fabs(Sd-S); 33 34 /** 計算結果の出力 **/ 35 printf(" きざみ幅 真値 台形 \n"); 36 printf("%12.9f %12.9f %12.9f %12.9f \n",h,S,Sd,Err_d); 37} 38 39/** 積分関数 **/ 40double f(double xx) 41{ 42 return 7*xx*log(xx); 43} 44 45/** 台形公式による積分 **/ 46double Sdai(int ni,double y[]) 47{ 48 double S; 49 int l; 50 51 S=0.0; 52 for(l=1; l<ni; l++) S+=2.0*y[l]; 53 S=(S+y[0]+y[ni])*h/2.0; 54 return S; 55}
試したこと
7*xを省いてlogのみでの計算も試したのですが同様のエラーメッセージが出ました。
補足情報(FW/ツールのバージョンなど)
ウェブコンパイラのRepl.itを使用しています。
回答2件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
退会済みユーザー
2019/01/24 17:32
2019/01/24 20:34
2019/01/25 01:21 編集
退会済みユーザー
2019/01/25 05:10