METODO DE BISECCION

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

 void main(void)
 {
   clrscr();
   float  cont, iter, a,b, Pint,fb,fPi, Ea, error, ftal;
   float c,d, fPin, fd, ftal2;
   cont=0;
   cout<<"\n\n\t\t  INSTITUTO POLITECNICO NACIONAL";
   cout<<"\n\n\t\t         ESIME ZACATENCO";
   cout<<"\n\n\t\t       Torres Martinez Gustavo Alan";
   cout<<"\n\n\t\t        ANALISIS NUMERICO";
   cout<<"\n\n\t\t           GRUPO: 4CM1"<<endl<<endl;

   cout<<endl<<"\n\n\tPROGRAMA PARA CALCULAR LAs RAICES DE LA FUNCION f(x)=(Cosx*Cosx)-2(Senx)";
   cout<<"\n\t\t  POR MEDIO DEL METODO DE BISECCION";
   cout<<endl<<"\nINDIQUE EL INTERVALO (a,b)\n";
   cout<<endl<<"a= ";
   cin>>a;
   cout<<endl<<"b= ";
   cin>>b;
   cout<<endl<<"INDIQUE EL ERROR DE APROXIMACION: ";
   cin>>error;
   Ea=1000;
while(Ea>error)
      {
Pint=(a+b)/(2);
      iter=cos(Pint);
      fPi=(cos(Pint)*cos(Pint))-(2*(sin(Pint)));
      fb=(cos(b)*cos(b))-(2*(sin(b)));
      ftal=(fPi*fb);
      if(ftal>0)
      {
b=Pint;
Ea=(a-b);
Ea=(Ea*Ea);
Ea=pow(Ea,0.5);
      }
      else
{
a=Pint;
     Ea=(a-b);
   Ea=(Ea*Ea);
Ea=pow(Ea,0.5);
}
      }
      cout<<"\n\n\n\tLa raiz es: "<<Pint;
      cout<<"\n\tCon un error de aproximacion de: "<<Ea;
   cout<<endl<<"\n\nAHORA INDIQUE EL NUEVO INTERVALO (c,d)\n";
   cout<<endl<<"c= ";
   cin>>c;
   cout<<endl<<"d= ";
   cin>>d;
   Ea=1000;
while(Ea>error)
      {
Pint=(c+d)/(2);
      fPin=(cos(Pint)*cos(Pint))-(2*(sin(Pint)));
      fd=(cos(d)*cos(d))-(2*(sin(d)));
      ftal=(fPin*fd);
      if(ftal>0)
      {
d=Pint;
Ea=(c-d);
Ea=(Ea*Ea);
Ea=pow(Ea,0.5);
      }
      else
{
c=Pint;
     Ea=(c-d);
   Ea=(Ea*Ea);
Ea=pow(Ea,0.5);
}
}
      cout<<"\n\n\n\tLa raiz es: "<<Pint;
      cout<<"\n\tCon un error de aproximacion de: "<<Ea;
      getche();
}