Foro de elhacker.net

Hola, hace poco coloqué una duda sobre el método de Bisección, ahora tengo uno parecido, el de "Regla Falsa":
http://esimecu-anumerico.blogspot.mx/2011/06/metodo-de-la-regla-falsa.html (http://esimecu-anumerico.blogspot.mx/2011/06/metodo-de-la-regla-falsa.html)

El problema es que al momento de dar un valor mayor a 1 en el extremo inferior (xi) me marca error para cualquier número de iteraciones, este es el código:

#include <iostream>
#include <stdlib.h>
#include <math.h>
 
using namespace std;
 
double eval_funcion(double);
 
int main(){
 float TOL, xi, xs, xr, xra;
 float fxi , fxs, fxr; 
 float validador, exr, er;
 float A, B, C, D;
 int N, i;
 
 cout<< "Ingrese el extremo inferior (Xi) : " << endl;
 cin >> xi; 
 cout << "Ingrese el extremo superior (Xs) : " << endl;
 cin >> xs; 
 cout << "Ingrese la tolerancia : " << endl;
 cin >> TOL; 
 cout << "Ingrese el n\xA3mero de iteraciones : " << endl;
 cin >> N; 
 
 fxi=eval_funcion(xi); 
 fxs=eval_funcion(xs);
 A = xi - xs;
 B = fxs * A;
 C = fxi - fxs;
 D = B/C;   
 
 xr = xs - D;   
 cout<<"\nAproximaci\xA2n de la ra\xA1z buscada: " << xr;    fxr = eval_funcion(xr);
 
 cout<<"\n"<<"n"<<"\t"<<"Ra\xA1z"<<"\t\t"<<"f(ra\xA1z)"<<"\t\t\t"<<"Error";
 
 while(i <= N)
 {
  xra = xr; 
  validador = fxi * fxr;
 
  if(validador < 0) {   
  	 xs=xr;    
  }
  else{
    xi=xr;   
  } 
 
  fxi=eval_funcion(xi); 
  fxs=eval_funcion(xs);
  A = xi - xs;
  B = fxs * A;
  C = fxi - fxs;
  D = B/C;
  xr = xs - D;
  fxr=eval_funcion(xr);    
 
  //error relativo
  exr = xr - xra;  
  er = exr/xr; 
  if(er<0){
     er=er*(-1);
  }
 
  //resultados iteracion a iteracion
  cout<<"\n"<<i<<"\t"<<xr<<"\t\t"<<fxr<<"\t\t"<<er;
 
  if((fxr==0)||(er<TOL)){
  	cout<<"\nEl proceso se realiz\xA2 de manera correta."; 
	  system("pause"); exit(1);
  }
 
  i++;
 } 
 
 if((i>N)&&(er>TOL)) {
   cout<<"\nEl m\x82todo fracas\xA2 despu\x82s de "<< N <<" iteraciones\n";
 }
   system("pause");     
}
/******************************************************************************/
 
double eval_funcion(double v)
{
 double r , ee, emx, logx;
 ee = pow(2.718281828,v); emx = 1/ee;
 logx = log(v);
 
 r = emx - logx; 
 return r;    
}
 

Por qué sólo funciona dando el valor de 1 en xi?, gracias de antemano.

intentaste hacer un debug del programa para ver linea a linea que pasa y si se comporta como esperas?

Sí, lo hice pero no logro captar ningún error, no entiendo por qué no funciona el programa, debería funcionar perfecto así como está :s

cuando hiciste el debug, calculaste en que linea el calculo matemático deja de ser correcto?

Pues, tal parece que el error se encontraba en esta línea:


y ya la modifiqué por esta otra:


Pero sigo sin poder solucionarlo, creo que así lo dejaré, igual y debo comenzar todo desde 0.