#include <stdio.h>
#include <math.h>
#include <conio.h>

//numero maximo de iteraciones
#define N      50
//tolerancia maxima
#define TOL 0.001
//funcion a aproximar
//f(v)=((1+(3.592/x^2))*(x - 0.04267) -(0.082054*300))
double f(double v)
 {
 //VALOR DE LA PRESION Y LA TEMPERATURA
 int P=100,T=300;
return ((P+(3.592/(pow(v,2))))*(v-0.04267))-(0.082054*T);
  }
 //derivada de la funcion f(x)
//f'(v)=1-(3.592/x^2)+(2*(3.592*0.04267)/x^3)
double df(double v){
//VALOR DE LA PRESION
  int P=100;
 return ((P-(3.592/(pow(v,2))))+((2*(3.592*0.04267))/(pow(v,3))));

}
int main(void)
{
double v, v0;
int i;
clrscr();
i=1;

//aproximacion inicial de la raiz f(v)
v0=24.6162;
printf("\n %d %15.12lf", i, v0);
while (i<N)
	{
	v=v0-f(v0)/df(v0);
	if (fabs(v-v0)<TOL)
		{
		//se alcanzo la tol pedida
		//y no se superaron las
		//de iteraciones, esta es la raiz.
		printf("\n Esta es la  raiz V=%15.5f", v);
		getch();
		return 0;
		};
	i+=1;
	v0=v;
	//muestra las sucesivas aproximaciones de la raiz.
	printf("\n %d %5.5f", i, v);
	}
// el metodo fallo.
printf("Se alcanzo el numero maximo de iteraciones");
getch();
return 0;
}