FRight = Y; } else {
Left = X;
FLeft = Y;
}
N++;
} while (fabs(Y) >= Eps);
return X;
}
#define FF(x) ( (M_PI * pow(x, M_PI + 7) + 2 * M_PI*pow(x, M_PI + 3) + M_PI * pow(x, M_PI - 1) + 4 * pow(x, 3)) / (pow(x, 8) + 2 * pow(x, 4) + 1) )
int main()
{
int k_B,k_H,k_N; long int s;
float a1, b1, eps1, delta1; double a, b, eps, x_B, x_H, x_N;
a = 0.5; b = 1.25;
double x = 0.867086; printf("eps\t\tdelta\t\ta\t\tb\t\tx_B\t\tx_H\t\tx_N\t\t\tk_B\tk_H\t
k_N\tDx\t\tc\n");
for (delta = 0.1; delta >= 0.000001; delta /= 10)
{
for (eps = 0.1; eps >= 0.000001; eps /= 10)
{
x_B = BISECT(a, b, eps, k_B); x_H = HORDA(a, b, eps, k_H); x_N = NEWTON(b, eps, k_N);
printf("%lf\t%lf\t%lf\t%lf\t%lf\t%lf\t%lf\t%d\t%d\t%d\t%lf\t%d\n", eps, delta, a, b, x_B, x_H, x_N, k_B, k_H, k_N, abs(x-x_N), eps>=abs(x - x_N));
}
}
return 0;
}
16