Hola, me podrían ayudar con el algoritmo de Pegasus, llevo varios días tratando y mi cabeza ya no da para más. Esto es.lo que llevo, espero puedan ayudarme

#include <stdio.h>
#include <math.h>
 
int main(void) {
 
 function{}=pegasus(f,a,b,delta,epsilon,maxi)
 fx=inline(f);
 ya=feval(fx,a);
 yb=feval(fx,b);
 if ya*tb>0{
 printf("f(a)+f(b)>0");
 break;
}
for (K=1,max1)
dx=yb+(b-a)/yb-ya);
c=b-dx;
alpha=b/(b+c);
printf("alpha=%015f\n",alpha);
 
ac=c-a;
yc=feval(fx,c);
if(yc==0){
break;
else
if yb+yc>0
b=c;
ya=(alpha)+ya;
yb=yc;
else
a=c;
ya=yc;
yb=(alpha)+yb;
}
dx=min(abs(dx),ac);
if (abs(dx)<delta){break}
if(yc)<epsilon
break;
err=abs(b-a)/2;
yc(feval)(fx,c);
printf("La solucion es %.15f, y se ha alcanzado en %d de interaciones\n", c, k);
printf("Con un error de %e\n", err );
printf("Con f(raiz)=%e\n", yc);
 
 
 return 0;
}