資源簡介
1.1 double gauss_ch1(double(*f)(double), int n);求積分∫_(-1)^1 f(x)dx/√(1-x^2 )
實現n點Gauss-Chebyeshev積分公式;返回積分的近似值。
在區間[-1,1]上關于權函數1/√(1-x^2 )的正交多項為T_n (x)=cos(narccos(x)),T_n (x)在[-1,1]上的n個根是x_k=cos?((2k-1)/2n π),k=1,…,n. n點Gauss-Chebyeshev積分公式為∫_(-1)^1 f(x)dx/√(1-x^2 )≈π/n ∑_(k=1)^n f(cos?((2k-1)/2n π))
1.2 double gauss_ch2(double(*f)(double), int n); 求積分∫_(-1)^1 √(1-x^2 ) f(x)dx
實現n點Gauss-Chebyeshev II型積分公式;返回積分的近似值。
在區間[-1,1]上關于權函數√(1-x^2 )的正交多項為U_n (x)=sin?((n+1)arccos?(x))/sin?(arccos?(x)) ,U_n (x)在[-1,1]上的n個根是x_k=cos?(kπ/(n+1)),k=1,…,n. n點Gauss-Chebyeshev II型積分公式為
∫_(-1)^1 √(1-x^2 ) f(x)dx≈π/(n+1) ∑_(k=1)^n sin^2 (kπ/(n+1))f(cos?(kπ/(n+1)))
1.3 double comp_trep(double (*f)(double), double a, double b);求積分∫_a^b f(x)dx
函數實現逐次減半法復化梯形公式;返回積分的近似值。
1.4 double romberg(double (*f)(double), double a, double b); 求積分∫_a^b f(x)dx
函數實現Romberg積分法;返回積分的近似值。
1.5 double gauss_leg_9(double (*f));求積分∫_(-1)^1 f(x)dx
實現9點Gauss-Legendre求積公式。
使用上面實現的各種求積方法求下面的積分:∫_(-1)^1 e^x √(1-x^2 ) dx (=∫_(-1)^1 (xe^x)/√(1-x^2 ) dx)
使用第3,4,5個函數求積分:∫_0^(π/2) sin?x dx (=1)
代碼片段和文件信息
- 上一篇:Esxi6.0 or 6.5 全套注冊機
- 下一篇:SmartCos工具
評論
共有 條評論