-
-
Save parzibyte/1c291f9669d4e72b30c50865282d8aec to your computer and use it in GitHub Desktop.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
/* | |
____ _____ _ _ _ | |
| _ \ | __ \ (_) | | | | |
| |_) |_ _ | |__) |_ _ _ __ _____| |__ _ _| |_ ___ | |
| _ <| | | | | ___/ _` | '__|_ / | '_ \| | | | __/ _ \ | |
| |_) | |_| | | | | (_| | | / /| | |_) | |_| | || __/ | |
|____/ \__, | |_| \__,_|_| /___|_|_.__/ \__, |\__\___| | |
__/ | __/ | | |
|___/ |___/ | |
Blog: https://parzibyte.me/blog | |
Ayuda: https://parzibyte.me/blog/contrataciones-ayuda/ | |
Contacto: https://parzibyte.me/blog/contacto/ | |
Copyright (c) 2020 Luis Cabrera Benito | |
Licenciado bajo la licencia MIT | |
El texto de arriba debe ser incluido en cualquier redistribución | |
=========================================== | |
Método de Newton en C++ | |
=========================================== | |
*/ | |
#include <cmath> | |
#include <iostream> | |
using namespace std; | |
float resolverEcuacion(float valor) { | |
// Con la ecuación: x^3 + 4x^2 -10 | |
return pow(valor, 3) + 4 * pow(valor, 2) - 10; | |
} | |
float resolverDerivada(float valor) { | |
// Con la ecuación: 2x^2 + 8x | |
return pow((2 * valor), 2) + 8 * valor; | |
} | |
int main() { | |
int iterador = 0; | |
float xi = 0.75; | |
float xi_xi = 0; | |
float ultimo_xi = 0; | |
printf("Método de Newton\n"); | |
printf("Problema: x^3 + 4x^2 -10\n"); | |
printf("Derivado: 2x^2 + 8x\n\n\n"); | |
printf("+----+-------------+-------------+-------------+-------------+-------" | |
"-----+\n"); | |
printf("+ i | xi | f(xi) | f'(xi) | xi + 1 | " | |
"xi_xi |\n"); | |
printf("+----+-------------+-------------+-------------+-------------+-------" | |
"-----+\n"); | |
while (1) { | |
float fxi = resolverEcuacion(xi); | |
float _fxi = resolverDerivada(xi); | |
float xi_1 = xi - (fxi / _fxi); | |
printf("|%3d |%12.8f |%12.8f |%12.8f |%12.8f |%12.8f|\n", iterador, xi, fxi, | |
_fxi, xi_1, xi_xi); | |
iterador++; | |
ultimo_xi = xi; | |
xi = xi_1; | |
xi_xi = abs(xi - ultimo_xi); | |
if (xi_xi == 0) { | |
printf("+----+-------------+-------------+-------------+-------------+---" | |
"---------+\n"); | |
break; | |
} | |
} | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment