Integração numérica - Método da Quadratura Gaussiana
Publicado por Daniel Moreira dos Santos (última atualização em 02/09/2009)
[ Hits: 14.515 ]
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); }
Enviar mensagem ao usuário trabalhando com as opções do php.ini
Meu Fork do Plugin de Integração do CVS para o KDevelop
Compartilhando a tela do Computador no Celular via Deskreen
Como Configurar um Túnel SSH Reverso para Acessar Sua Máquina Local a Partir de uma Máquina Remota
Configuração para desligamento automatizado de Computadores em um Ambiente Comercial
Compartilhamento de Rede com samba em modo Público/Anônimo de forma simples, rápido e fácil
Cups: Mapear/listar todas as impressoras de outro Servidor CUPS de forma rápida e fácil
Criando uma VPC na AWS via CLI
Erro no upgrade: Sub-process /usr/bin/dpkg returned an error code (1) (2)
Falta pacotes de suporte ao sistema de arquivos (Gerenciador de discos... (6)
A area de trabalho ficou preta (16)
Agora temos uma assistente virtual no fórum!!! (244)
Servidor said: 530 5.7.0 Must issue a STARTTLS command first (in r... (0)