Código fonte

Código Fonte Métodos 2° Regra de Simpson e Trapézio para Integral Dupla

Código Runge – Kutta
Para Sistemas de Equação Diferencial

Pontifícia Universidade Católica de Minas Gerais
Unidade São Gabriel
Engenharia de Computação

Aluno: Thiago Corrêa da Silveira
Professor: Paulo Fernando

Disciplina: Cálculo Numérico

Código Fonte Métodos 2° Regra de Simpson e Trapézio para Integral Dupla

/*PONTIFÍCIA UNIVERSIDADE CATÓLICA DE MINAS GERAIS

Trabalho de Métodos Numéricos para a computação Integração Dupla

ALUNO: Thiago Corrêa da Silveira
MATRÌCULA: 383737
*/

//----------------------------------- BIBLIOTECAS--------------------------------------
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<ctype.h>
#include<conio.h>

//-------------------PROTÓTIPOS DOS MÉTODOS E FUNÇÕES UTILIZADOS-----------------------
//double f(double x,double y);
//double Trapezio(double subint,int coef[50],double x[50],double matriz[50][50],double h,int posi, bool boInt);
//double Simpson(int regra,double subint,int coef[50],double x[50],double matriz[50][50],double h,int posi, bool boInt);
//-------------------------------------------------------------------------------------

//--------------------------------VARIÁVEIS GLOBAIS------------------------------------- double X[50],Y[50],resp[50][50],ptsX[50],ptsY[50];//vetores para guardar resultados das operações e resultados parciais que serão exibidos na tabela

int coeficientes[50],coeficientesX[50]; //vetores para armazenar os coeficientes que multiplicarão os termos do somatório de um método (Trapézio ou Simpson
//------------------------------------------------------------------------------------------------

//*************************** FUNÇÃO A SER INTEGRADA ****************************

double f(double x, double y)
{
return( sin(x+y) );//Informe a função a ser integrada dentro do "return"
}

double Trapezio(double subint,int coef[50],double x[50],double matriz[50][50],double