Método de Newton-Raphson
Publicado por Renan Birck Pinheiro (última atualização em 27/05/2010)
[ Hits: 35.488 ]
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); }
Exceções em C através de try-throw-catch
Resolver problemas de Internet
Como compartilhar a tela do Ubuntu com uma Smart TV (LG, Samsung, etc.)
Descritores de Arquivos e Swappiness
Solução rápida para o problema do Network Manager conectar mas não navegar
Como instalar no Linux Jogos da Steam só para Windows
Instalando o Team Viewer no Debian Trixie - problema no Policykit
Como rodo essa suinaria? [RESOLVIDO] (7)
Montando e usando iso de um sistema dentro do outro (2)
Criar atalho para uma pasta na area de trabalho no Linux Mint? (0)
Precisa de ajuda com monitoramento ou problemas com Zabbix? (0)