#include<stdio.h>
#include<conio.h>
#include<math.h>
#include<stdlib.h>
#define a -1.70555
#define b 13181.2
#define c -2.58768
#define d 2.4475E-03
#define e -3.1144E-07

double min (double aa, double bb)
{
double cc;
if (aa<bb)
cc=aa;
else cc=bb;
return (cc);
}

double bisection (double x[3], double K, double p)
{
int i;
static double m, n, N, F, Fm, y;
m=-x[3]/2;
n=min((x[1]/2), x[2]);
i=0;
N=fabs(log(fabs((m-n)/pow10(18))))/log(2);
printf("%.18lf",n);
do
  {y=m+(n-m)/2;
   F=exp(y)-2;
   Fm=exp(m)-2;
   if ((F*Fm)<0)
   n=y;
   else m=y;
   i++;
  }
while (i<N);
printf("\n%d",i);
printf("\n%.17lf",y);
return (y);
}

void main (void)
{
int i, j, xx, yy;
static double x[3], y, T, p, K;
printf("Enter the temperature: ");
scanf("%lf",&T);
printf("Enter the pressure: ");
scanf("%lf",&p);
printf("Enter the quantities of reagents: \n\nNO\tO2\tNO2\n\n");
for (i=1; i<=3; i++)
{xx=wherex();
 yy=wherey();
 gotoxy(xx, yy-1);
 scanf("%lf",&x[i]);
 for (j=1; j<=i; j++)
 printf("\t");
}
K=a+b/T+c*log(T)+d*T+e*T*T;
K=exp(K);
y=bisection(x,K,p);
}