X1 = X2;
X2 = PHI(X1);
}
return X2;
}
#endif
#ifdef __NEWTON
double NEWTON(double X, double Eps, int &N) { extern double F1(double);
double Y, Y1, DX, Eps0; N = 0;
double m1 = 0.0, // наименьшее значение модуля 1-ой производной M2 = 0.0; // наибольшее значение модуля 2-ой производной
Eps0 = sqrt(2 * m1 * Eps / M2);
do {
Y = F(X);
if (Y == 0.0) { return X;
}
Y1 = F1(X);
if (Y1 == 0.0) {
puts("Производная обратилась в ноль\n"); exit(1);
}
DX = Y / Y1; X -= DX; N++;
} while (fabs(DX) > Eps0);
return X;
}
#endif
double HORDA(double Left, double Right, double Eps, int &N) { double FLeft = F(Left);
double FRight = F(Right);
11
double X, Y;
if (FLeft * FRight > 0.0) { puts("Неверное задание интервала\n"); exit(1);
}
if (Eps <= 0.0) {
puts("Неверное задание точности\n"); exit(1);
}
N = 0;
if (FLeft == 0.0) { return Left;
}
if (FRight == 0.0) { return Right;
}
do {
X = Left - (Right - Left) * FLeft / (FRight - FLeft); Y = F(X);
if (Y == 0.0) { return X;
}
if (Y * FLeft < 0.0) { Right = X;
FRight = Y; } else {
Left = X;
FLeft = Y;
}
N++;
} while (fabs(Y) >= Eps);
return X;
}
12
#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; long int s;
float a1, b1, eps1, delta1; double a, b, eps, x_B,x_H;
a = 0.5; b = 1;
printf("eps\t\tdelta\t\ta\t\tb\t\tx_B\t\tx_H\t\tn_m\t\tn\t\t\tk_B\t k_H\tk_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);
printf("%lf\t%lf\t%lf\t%lf\t%lf\t%lf\t%lf\t%lf\t%d\t%d\t%d\t%d\n", eps, delta, a, b, x_B, x_H, eps / delta, 1 / (double)FF(x_H), k_B, k_H, (int)round(log2((b - a) / eps)), eps / delta> 1 / FF(x_H));
}
}
return 0;
}
13