Raiz cúbica pelo método de bissecção
Publicado por Gabriel (última atualização em 15/05/2010)
[ Hits: 13.234 ]
Resolve numericamente a equação x³ - n = 0 utilizando o método de bissecção. O objetivo é calcular a raiz cúbica de n.
Obs.: Foi utilizada precisão dupla e como critério de parada foi utilizado o DBL_EPSILON (2.220446e-16).
Problema: para -1 <= n <= 1 não é possível calcular a raiz cúbica (com esse script)... Alguém sabe porquê?
#include <stdio.h>
#include <float.h>
#include <math.h>
double cubica(double n);
int main(){
double n = 0;
printf("Entre com um número: ");
scanf("%lf", &n);
printf("%lf\n", cubica(n));
return 0;
}
// calcular raiz cúbica de n e para isso
// encontrar x tal que x³ - n = 0
double cubica(double n) {
double x1 = 0, x2 = 0, xm = 0;
if (n > 0){
x1 = -n;
x2 = n;
} else {
x1 = n;
x2 = -n;
}
do {
xm = (double) (x2 + x1) / 2;
if (xm*xm*xm - n > 0) {
if (x1*x1*x1 - n > 0) x1 = xm;
else x2 = xm;
} else {
if (x1*x1*x1 - n < 0) x1 = xm;
else x2 = xm;
}
} while (fabs(x2 - x1) > DBL_EPSILON);
return xm;
}
Passando uma matriz para funcao
Desenhando uma curva Dragão ou o Fractal Jurassic Park
Papagaiando o XFCE com temas e recursos
WhatsApp com Chamadas no Linux via Waydroid
XFCE - quase um Gnome ou Plasma mas muito mais leve
LXQT - funcional para máquinas pererecas e usuários menos exigentes
Usando o --load-average no Gentoo para evitar OOM Kill
Instalando o COSMIC DE no Gentoo (systemd)
Centralizar Logo com Transparência via ImageMagick
Problemas com o "startx&... no Slackware (11)
Ajudante para criar dicas e artigos no VOL (3)









