三点高斯积分程序

最新推荐文章于 2024-01-03 02:59:19 发布

Y___Y

最新推荐文章于 2024-01-03 02:59:19 发布

阅读量7.9k

点赞数

分类专栏：数值计算程序专栏文章标签： fun math.h ini

本文链接：https://blog.csdn.net/Y___Y/article/details/1788593

版权

这是一个C语言实现的高斯积分程序，适用于一维和二维的积分计算。程序利用三点高斯公式，通过调用用户定义的函数进行数值积分，具有高效、精确的特点。示例中展示了如何对自定义函数进行积分验证。

摘要由CSDN通过智能技术生成

//高斯积分程序，高斯积分具有计算速度快，精度高，能计算反常积分等优点
#include <stdio.h>
#include <math.h>

#define PI 3.14159265
static const double gp3[3]={-0.77459666924148338, 0.0                , 0.77459666924148338};//高斯点
static const double gc3[3]={ 0.55555555555555556, 0.88888888888888889, 0.55555555555555556};//系数

//三点Guass积分
double Guass_Integral(double x1,double x2,double (*fun)(double),int n)
{
    int i;
    double dx,k,l,r=0.0;
    dx=(x2-x1)/n;
    k=dx/2;
    for(i=0;i<n;i++)
    {
        l=x1+(i+0.5)*dx;
        r+=gc3[0]*fun(k*gp3[0]+l)+gc3[1]*fun(k*gp3[1]+l)+gc3[2]*fun(k*gp3[2]+l);
    }
    return r*k;
}

double Guass_Integral(double u1,double u2,double v1,double v2,void (*fun)(double,double,double *,double *,double *),int m,int n)
{