Método de Newton-Raphson
Publicado por Renan Birck Pinheiro (última atualização em 27/05/2010)
[ Hits: 35.864 ]
Homepage: http://renanbirck.rocks
Implementação na linguagem C do método de Newton-Raphson, usado para obtenção numérica do zero de funções. Não pretendo explicar o funcionamento do método aqui; recomendo a leitura de um livro de Cálculo Numérico para tal fim.
float f(float x) é a função cujo zero queremos descobrir.
float der(float x) é a derivada da função f(x); uma derivada errada tornará a convergência lenta ou impossível.
#include <stdio.h>
#include <math.h>
float f(float x) {
return pow(x,2)-5*x+6;
}
float der(float x) {
return 2*x-5; // derivada de f(x).
}
int main() {
float eps, x0, *iter; // erro, ponto inicial, vetor iterações.
int i, numiter; // iteração atual, número de iterações.
printf("Método de Newton-Rhapson para o zero da função f(x).\n");
printf("Digite o erro: \n");
scanf("%f",&eps);
printf("Digite o número máximo de iterações?\n");
scanf("%d", &numiter);
printf("Digite o X0 inicial? \n");
scanf("%f", &x0);
// Alocar dinâmicamente memória para o vetor das iterações.
iter = malloc(sizeof(float) * numiter);
// Condições iniciais.
iter[0] = x0;
i = 0;
// Iterações.
while(f(iter[i]) > eps) {
// Excedeu o nosso limite de iterações.
if(i > numiter) {
printf("Não convergiu em %d iterações!!!\n", numiter);
printf("Provavelmente f'(x) está errada.\n");
}
iter[i+1] = iter[i] - f(iter[i])/der(iter[i]);
i++;
}
printf("X ~= %f ", iter[i]);
printf("\nForam feitas %d iterações.\n",i);
}
Calculo do número de Euler em C++
Crud em C++ orientado a objetos com banco de dados MySQL
Cirurgia para acelerar o openSUSE em HD externo via USB
Void Server como Domain Control
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
Como impedir exclusão de arquivos por outros usuários no (Linux)
Cirurgia no Linux Mint em HD Externo via USB
Anúncio do meu script de Pós-Instalação do Ubuntu
Alguém executou um rm e quase mata a Pixar! (5)
Formas seguras de instalar Debian Sid (4)
Duas Pasta Pessoal Aparecendo no Ubuntu 24.04.3 LTS (12)
Alguém pode me indicar um designer freelancer? [RESOLVIDO] (5)
Por que passar nas disciplinas da faculdade é ruim e ser reprovado é b... (6)









