Método de Newton-Raphson
Publicado por Renan Birck Pinheiro (última atualização em 27/05/2010)
[ Hits: 35.531 ]
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); }
Sintetizador sonoro que toca redez vouz 4
Aprenda a Gerenciar Permissões de Arquivos no Linux
Como transformar um áudio em vídeo com efeito de forma de onda (wave form)
Como aprovar Pull Requests em seu repositório Github via linha de comando
Aplicativo simples para gravar tela
Quebra de linha na data e hora no Linux Mint
Firefox não abre em usuário não administradores (0)
Sempre que vou baixar algum pacote acontece o erro dpkg (8)
tentando instalar em um notebook antigo o Linux LegacyOS_2023... [RESO... (8)
Problema com Conexão Outlook via Firewall (OpenSUSE) com Internet Fibr... (5)