Integração numérica - Método da Quadratura Gaussiana
Publicado por Daniel Moreira dos Santos (última atualização em 02/09/2009)
[ Hits: 14.768 ]
Homepage: http://www.danielmoreira.wordpress.com
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.
#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 argumentos com getopt_long
Teste de desempenho com números primos em C (corrigido)
Modo Simples de Baixar e Usar o bash-completion
Monitorando o Preço do Bitcoin ou sua Cripto Favorita em Tempo Real com um Widget Flutuante
Instalando partes faltantes do Plasma 6
Adicionar botão "mostrar área de trabalho" no Zorin OS
Como montar um servidor de backup no linux
Estou tentando ser legalista, mas tá complicado! (9)
espelhar monitores nao funciona (2)
SQLITE não quer funcionar no LINUX LMDE6 64 com Lazaruz 4.2 64bit (n... (1)









