Método de Newton-Raphson
Publicado por Renan Birck Pinheiro (última atualização em 27/05/2010)
[ Hits: 35.222 ]
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); }
Ajuda para calcular 5 Notas! Não Consigo colocar os gets na formula
Exceções em C através de try-throw-catch
Raizes reais e complexas de uma equação de 2º grau
Imprime a soma dos numeros positivos e negativos
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
Efeito "livro" em arquivos PDF
Como resolver o erro no CUPS: Unable to get list of printer drivers
Flatpak: remover runtimes não usados e pacotes
Mudar o gerenciador de login (GDM para SDDM e vice-versa) - parte 2
Como atualizar o Debian 8 para o 10 (10)
Dica sobre iptables ACCEPT e DROP (6)
NGNIX - Aplicar SNAT para evitar roteamento assimetrico (29)
[Python] Automação de scan de vulnerabilidades
[Python] Script para analise de superficie de ataque
[Shell Script] Novo script para redimensionar, rotacionar, converter e espelhar arquivos de imagem
[Shell Script] Iniciador de DOOM (DSDA-DOOM, Doom Retro ou Woof!)
[Shell Script] Script para adicionar bordas às imagens de uma pasta