Algoritmo Para Aproximação do Valor de Derivadas

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

// Programa Para Cálculo da Derivada em Determinado Ponto
// Autor: Rafael Yuri Medeiros Barbosa - Graduando em Engenharia Química
//=========================================================================================
//=========================================================================================
// Inclusão Bibliotecas
//=========================================================================================
#include <iostream>
#include <fstream>
#include <math.h>
#include <stdlib.h>
//=========================================================================================
//=========================================================================================
// Definição de Constantes
//=========================================================================================
#define ep 1e-5
//=========================================================================================
//=========================================================================================
// Definição do Padrão C++
//=========================================================================================
using namespace std;
//=========================================================================================
//=========================================================================================
// Função Primitiva
//=========================================================================================
double fx ( double x )
{
 return ( pow ( x , 3 ) ) ;
}
//=========================================================================================
//=========================================================================================
// Função Derivada
//=========================================================================================
double derivada ( double x )
{
 double d ;
 d = ( fx ( x + ep ) - fx ( x ) ) / ep ;
 return d ;
}
//=========================================================================================
//=========================================================================================
// Função Derivada Segunda
//=========================================================================================
double derivadaseg ( double x )
{
 double d2 ;
 d2 = ( fx ( x + ep ) - 2 * fx ( x ) + fx ( x - ep ) ) / ( pow ( ep , 2 ) ) ;
 return d2 ;
}
//=========================================================================================
//=========================================================================================
// Função main ( )
//=========================================================================================
int main ()
{
 double x ;
 fstream arquivo ( "resultados.dat" , ios::in | ios::out | ios::trunc ) ;
 cout << "Digite o valor do ponto onde sera calculada a derivada: ";
 cin >> x ;
 arquivo << "Derivada primeira: " << derivada ( x ) << "\n" ;
 arquivo << "Derivada segunda: " << derivadaseg ( x ) << "\n" ;
 arquivo.close ( ) ;
 return 0 ;
}
//=========================================================================================
//=========================================================================================
// Fim do Programa
//=========================================================================================