#include
#include
double bisek(double a,double b,double(*f)(double x),double eps,int *it);
double fx(double x);
double fm(double m);
double pp,gam;
void main()

{ double a,b,eps,x,m; int k;
//printf("podaj pp i gam:\n");
//scanf("%lf%lf",&pp,&gam);
gam=1.4;
pp=0.2724;
printf("podaj a i b:\n");
scanf("%lf%lf",&a,&b);
printf("podaj eps:\n");
scanf("%lf",&eps);
//x=bisek(a,b,fx,eps);
m=bisek(a,b,fm,eps,&k);
printf("m=%lf l=%d\n",m,k);
}
double bisek(double a,double b,double(*f)(double x),double eps,int *it)
// PROCEDURA BISEKCJI DO ROZWIAZYWANIA ROWNAN NIELINIOWYCH f(x)=0
// f- ZADANA FUNKCJA a,b-PRZEDZIAL PIERW. eps-DOKLADNOSC
{
double x;int l;
l=0;
do
{
x=0.5*(a+b);
l=l+1;
if(f(a)*f(x)<0.0) b=x;
else a=x;
*it=l;
}
while(fabs(b-a)>eps);

return x;
}



double fx(double x)
{double f;
f=cos(x)-x;
return f;
}
double fm(double m)
{ double r,q;
r=1.+(gam-1.)/2.*m*m;
q=-gam/(gam-1.);
return pow(r,q)-pp;
}
double sieczne(double a,double b,double(*f)(double x),double eps)
// PROCEDURA BISEKCJI DO ROZWIAZYWANIA ROWNAN NIELINIOWYCH f(x)=0
// f- ZADANA FUNKCJA a,b-PRZEDZIAL PIERW. eps-DOKLADNOSC
{
double x;
do
{
x=b-f(b)8(b-a)/(f(b)-f(a));
if(f(a)*f(x)<0.0) b=x;
else a=x;
}
while(fabs(b-a)>eps);
return x;
}