Temas

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;
}