高斯-勒让德公式求积分

 1 #include <iostream>
 2 #include <cmath>
 3 using namespace std;
 4 //计算f(x)=1/x 在[1,3]上的积分
 5
 6 double f(double x){
 7     return 1.0/(x+2);
 8 }
 9
10 double f1(double x){
11     return 1.0/(x+5);
12 }
13 double f2(double x){
14     return 1.0/(x+7);
15 }
16
17 double f3(double x){
18     return 1.0/(x+9);
19 }
20 double f4(double x){
21     return 1.0/(x+11);
22 }
23
24 int main(){
25     //运用三点高斯公式积分值
26     printf("%2.8f\n",    (0.5555556*f( -0.7745967)+0.8888889*f( 0 )+0.5555556*f( 0.7745967) ));
27
28     //运用五点高斯公式求积分值
29     printf("%2.8f\n",  0.2369269*f(-0.9061799)+ 0.2369269*f(0.9061799)+0.4786287*f(-0.5384693)+0.4786287*f(0.5384693) +
30     0.5688889*f(0) );
31     printf("\n");
32
33     //将区间4等分，分别利用二点高斯公式，然后累加得积分值
34         printf("%2.8f\n",f1( -1.0/sqrt(3.0))+f1(1.0/sqrt(3.0))  );
35         printf("%2.8f\n",f2( -1.0/sqrt(3.0))+f2(1.0/sqrt(3.0))  );
36         printf("%2.8f\n",f3( -1.0/sqrt(3.0))+f3(1.0/sqrt(3.0))  );
37         printf("%2.8f\n",f4( -1.0/sqrt(3.0))+f4(1.0/sqrt(3.0))  );
38
39         printf("sum=%2.8f\n",f1( -1.0/sqrt(3.0))+f1(1.0/sqrt(3.0)) +f2( -1.0/sqrt(3.0))+f2(1.0/sqrt(3.0)) +f3( -1.0/sqrt(3.0))+f3(1.0/sqrt(3.0)) +f4( -1.0/sqrt(3.0))+f4(1.0/sqrt(3.0))  );
40
41 }

转载于:https://www.cnblogs.com/liugl7/p/4963500.html

高斯-勒让德公式求积分相关推荐

复合梯形公式与复合辛普森公式求积分
一实验目的 1. 掌握复合梯形公式与复合辛普森公式的基本思想. 2. 编程实现用复合梯形公式与复合辛普森公式求积分. 3. 熟悉matlab软件的使用. 二实验内容 1.用复合梯形公式计算积分 I ...
重复高斯勒让德法则(gauss-legendre)求积分(python,数值积分)
第四十一篇高斯勒让德法则牛顿-科特斯法则很方便使用,因为样本点在积分范围内均匀分布,而且权重系数很容易记住(至少到辛普森法则).详情可见重复牛顿-科特斯法则高斯法则采取样本点在积分范围内的最佳间 ...
高斯-勒让德求积分的Python程序
理论知识: 勒让德多项式及性质第四节高斯(Gauss)求积公式高斯-勒让德求积公式及Matlab实现目标函数求积区间为[-1,1]时代码: #高斯-勒让德求积公式 from sympy im ...
高斯勒让德求积公式matlab通用程序,数值分析——Gauss-Legendre 求积分(Matlab实现)...
2020 4.4 题目:编写Gauss求积法计算积分的程序(Gauss点数取1,2,3,4,5即可)并用于计算积分 ∫10sinxxdx∫01sin⁡xxdx\int_{0}^{1} \frac{\s ...
数据积分-牛顿科茨法与高斯勒让德法对比及示例
在结点xi上插值f的次数最多是n次的多项式 p(x)=∑f(xi)li(x) 拉格朗日插值多项式: li(x)= ∏(x-xj)/(xi-xj) 基于插值的数值积分: ∫f(x)dx≈∫p(x)dx= ...
高斯-勒让德积分——MATLABPython混合实现
高斯-勒让德积分前言一种很NB的数值积分.原理我们就不多说了,这里只介绍怎么直接拿来用. 核心问题是如何获取系数表. 参数介绍对于函数 f(x),积分区间为 a 到 b,求在这一区间上的定积分. ...
GSL科学计算库——计算高斯-勒让德积分
相关文章: Windows系统Qt5配置GSL科学计算库 Windows系统下Cygwin+Dev C ++ 配置GSL科学计算库假设计算下列积分: ∫0πexcos(x)dx\int_0^\pi ...
MATLAB高斯-勒让德积分求解
MATLAB高斯-勒让德积分求解积分函数 function [ y ] = CalcuFunctionValue( x ) y=exp(-x^2);%积分函数 end 求积分算法函数 functio ...
高斯-勒让德积分中不同阶数下最大高斯节点间距的关系
高斯-勒让德计算方法的高斯节点是勒让德多项式的根值,很多教材都把常用的根值和其对应的加权系数表示在如下的表格(P403 SEC. 7.5 GAUSS-LEGENDRE INTEGRATION (OPT ...

高斯-勒让德公式求积分

高斯-勒让德公式求积分相关推荐

最新文章

热门文章

高斯-勒让德公式 求积分

高斯-勒让德公式 求积分相关推荐

最新文章

热门文章

高斯-勒让德公式求积分

高斯-勒让德公式求积分相关推荐