A regra da quadratura gaussiana é uma aproximação da integral de uma função, geralmente estabelecido como um somatório com pesos dos valores assumidos pela função em pontos específicos dentro do domínio de integração. (Descrição adaptada de Wikipédia)

Você pode encontrar facilmente uma tabela de pontos na internet para usar com este programa.

• Download quad_gauss.c
• Enviar nova versão

Esconder código-fonte

#include <stdio.h>
#include <stdlib.h>

float eval_function(float x){
   return (exp(-(x^2)/2))
}

int main(void){

float a, b, sum=0, point, weight, h, jc;
int n i;
FILE *p;

printf("Extremo esquerdo do intervalo: ");
scanf("%f", &a);
printf("Extremo direito do intervalo: ");
scanf("%f", &b);
printf("Numero de pontos (1<=n<=8): "); /* com apenas 8 pontos temos uma aproximação excelente da integral da função acima */
scanf("%d", &n);

p=fopen("tabela.txt", "r");
if(!p){
  printf("Erro!");
  exit(1);
}

jc=(b-a)/2;
h=(b+a)/2;

for(i=0; i<n; i++){
   fscanf(p, "%f\t%f", &point, &weight);
   sum=sum+ (weight*jc*eval_function(jc*point+h));

}

system("pause");
return (0);
}

Manipulando árvores.

Agenda eletrônica feita em C

Crud em C++ orientado a objetos com banco de dados MySQL

Calculadora

[C] Parâmetros para printf