#include <iostream>
#include <cmath>
using namespace std;
// definicja funkcji
double f(double x)
{
return x*x*x - 2;
}
int main()
{
double a, b, srodek;
double eps = 0.000001; // dokladnosc
// sprawdzenie warunku metody bisekcji
if (f(a) * f(b) >= 0)
{
cout << "Brak miejsca zerowego w tym przedziale." << endl;
return 0;
}
// petla bisekcji
while ((b - a) / 2 > eps)
{
srodek = (a + b) / 2;
if (f(srodek) == 0)
break;
else if (f(a) * f(srodek) < 0)
b = srodek;
else
a = srodek;
}
srodek = (a + b) / 2;
cout << "Przyblizone miejsce zerowe: " << srodek << endl;
return 0;
}
I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y21hdGg+Cgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoKLy8gZGVmaW5pY2phIGZ1bmtjamkKZG91YmxlIGYoZG91YmxlIHgpCnsKICAgIHJldHVybiB4KngqeCAtIDI7Cn0KCmludCBtYWluKCkKewogICAgZG91YmxlIGEsIGIsIHNyb2RlazsKICAgIGRvdWJsZSBlcHMgPSAwLjAwMDAwMTsgLy8gZG9rbGFkbm9zYwoKICAgIAogICAgLy8gc3ByYXdkemVuaWUgd2FydW5rdSBtZXRvZHkgYmlzZWtjamkKICAgIGlmIChmKGEpICogZihiKSA+PSAwKQogICAgewogICAgICAgIGNvdXQgPDwgIkJyYWsgbWllanNjYSB6ZXJvd2VnbyB3IHR5bSBwcnplZHppYWxlLiIgPDwgZW5kbDsKICAgICAgICByZXR1cm4gMDsKICAgIH0KCiAgICAvLyBwZXRsYSBiaXNla2NqaQogICAgd2hpbGUgKChiIC0gYSkgLyAyID4gZXBzKQogICAgewogICAgICAgIHNyb2RlayA9IChhICsgYikgLyAyOwoKICAgICAgICBpZiAoZihzcm9kZWspID09IDApCiAgICAgICAgICAgIGJyZWFrOwogICAgICAgIGVsc2UgaWYgKGYoYSkgKiBmKHNyb2RlaykgPCAwKQogICAgICAgICAgICBiID0gc3JvZGVrOwogICAgICAgIGVsc2UKICAgICAgICAgICAgYSA9IHNyb2RlazsKICAgIH0KCiAgICBzcm9kZWsgPSAoYSArIGIpIC8gMjsKCiAgICBjb3V0IDw8ICJQcnp5Ymxpem9uZSBtaWVqc2NlIHplcm93ZTogIiA8PCBzcm9kZWsgPDwgZW5kbDsKCiAgICByZXR1cm4gMDsKfQo=