Programación en C++: método de iteración de Newton para resolver la ecuación xe^x-1=0

#include "iostream.h"

#include

#include

# incluir

#define N 100

#define PS 1e-5

#define TA 1e-5

float Newton(flotador (*f)(flotador),flotador(*f1)(flotador),flotador x0 )

{ float x1,d=0;

int k=0 ;

hacer

{ x1= x0-f(x0)/f1(x0);

if((k++>N)||(fabs( f1(x1))

{ cout<<"\n¡Error!";

getch();

salir(0);

}

d=(fabs(x1)<1?x1-x0:(x1-x0)/x1);

x0=x1;

cout<<"x("<

}

mientras((fabs(d) )>PS&&fabs(f(x1))>TA) ;

return x1;

}

flotador f(flotador x)

{ return x*exp(x)-1 }

float f1(float x)

{ return x*exp(x)+exp(x }

void main()

{ float f(flotador);

flotador f1(flotador);

flotador x0,y0; p>

cout<<"Entrada x0: ";

cin>>x0;

cout<<"x(0)="<

y0=Newton(f,f1,x0);

cout<<"\nLa raíz es x="<

getch();

}