#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;
double f(double x) {
return pow(x, 3) - x - 3;
}
int main() {
double x1, x2, x3, f1, f2;
const double tolerance = 1e-6;
const int max_iterations = 100;
int iteration = 1;
cout << "x1 এর মান প্রবেশ করুন: ";
cin >> x1;
cout << "x2 এর মান প্রবেশ করুন: ";
cin >> x2;
cout << fixed << setprecision(6);
cout << "\nIteration\tx1\t\tx2\t\tx3\t\tf(x1)\t\tf(x2)\n";
while (iteration <= max_iterations) {
f1 = f(x1);
f2 = f(x2);
if (fabs(f2 - f1) < 1e-12) {
cout << "বিভাজকের মান খুব ছোট। শূন্য দ্বারা বিভাজনের ঝুঁকি।\n";
break;
}
x3 = (f2 * x1 - f1 * x2) / (f2 - f1);
cout << iteration << "\t\t" << x1 << "\t" << x2 << "\t" << x3 << "\t" << f1 << "\t" << f2 << "\n";
if (fabs((x3 - x2) / x3) < tolerance) {
cout << "\nআনুমানিক মূল = " << x3 << endl;
break;
}
x1 = x2;
x2 = x3;
iteration++;
}
if (iteration > max_iterations) {
cout << "সর্বাধিক পুনরাবৃত্তির মধ্যে পদ্ধতিটি সন্নিবেশিত হয়নি।\n";
}
return 0;
}
I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y21hdGg+CiNpbmNsdWRlIDxpb21hbmlwPgoKdXNpbmcgbmFtZXNwYWNlIHN0ZDsKCmRvdWJsZSBmKGRvdWJsZSB4KSB7CiAgICByZXR1cm4gcG93KHgsIDMpIC0geCAtIDM7Cn0KCmludCBtYWluKCkgewogICAgZG91YmxlIHgxLCB4MiwgeDMsIGYxLCBmMjsKICAgIGNvbnN0IGRvdWJsZSB0b2xlcmFuY2UgPSAxZS02OwogICAgY29uc3QgaW50IG1heF9pdGVyYXRpb25zID0gMTAwOwogICAgaW50IGl0ZXJhdGlvbiA9IDE7CgogICAgY291dCA8PCAieDEg4KaP4KawIOCmruCmvuCmqCDgpqrgp43gprDgpqzgp4fgprYg4KaV4Kaw4KeB4KaoOiAiOwogICAgY2luID4+IHgxOwogICAgY291dCA8PCAieDIg4KaP4KawIOCmruCmvuCmqCDgpqrgp43gprDgpqzgp4fgprYg4KaV4Kaw4KeB4KaoOiAiOwogICAgY2luID4+IHgyOwoKICAgIGNvdXQgPDwgZml4ZWQgPDwgc2V0cHJlY2lzaW9uKDYpOwogICAgY291dCA8PCAiXG5JdGVyYXRpb25cdHgxXHRcdHgyXHRcdHgzXHRcdGYoeDEpXHRcdGYoeDIpXG4iOwoKICAgIHdoaWxlIChpdGVyYXRpb24gPD0gbWF4X2l0ZXJhdGlvbnMpIHsKICAgICAgICBmMSA9IGYoeDEpOwogICAgICAgIGYyID0gZih4Mik7CgogICAgICAgIGlmIChmYWJzKGYyIC0gZjEpIDwgMWUtMTIpIHsKICAgICAgICAgICAgY291dCA8PCAi4Kas4Ka/4Kat4Ka+4Kac4KaV4KeH4KawIOCmruCmvuCmqCDgppbgp4Hgpqwg4Kab4KeL4Kaf4KWkIOCmtuCnguCmqOCnjeCmryDgpqbgp43gpqzgpr7gprDgpr4g4Kas4Ka/4Kat4Ka+4Kac4Kao4KeH4KawIOCmneCngeCmgeCmleCmv+ClpFxuIjsKICAgICAgICAgICAgYnJlYWs7CiAgICAgICAgfQoKICAgICAgICB4MyA9IChmMiAqIHgxIC0gZjEgKiB4MikgLyAoZjIgLSBmMSk7CgogICAgICAgIGNvdXQgPDwgaXRlcmF0aW9uIDw8ICJcdFx0IiA8PCB4MSA8PCAiXHQiIDw8IHgyIDw8ICJcdCIgPDwgeDMgPDwgIlx0IiA8PCBmMSA8PCAiXHQiIDw8IGYyIDw8ICJcbiI7CgogICAgICAgIGlmIChmYWJzKCh4MyAtIHgyKSAvIHgzKSA8IHRvbGVyYW5jZSkgewogICAgICAgICAgICBjb3V0IDw8ICJcbuCmhuCmqOCngeCmruCmvuCmqOCmv+CmlSDgpq7gp4LgprIgPSAiIDw8IHgzIDw8IGVuZGw7CiAgICAgICAgICAgIGJyZWFrOwogICAgICAgIH0KCiAgICAgICAgeDEgPSB4MjsKICAgICAgICB4MiA9IHgzOwogICAgICAgIGl0ZXJhdGlvbisrOwogICAgfQoKICAgIGlmIChpdGVyYXRpb24gPiBtYXhfaXRlcmF0aW9ucykgewogICAgICAgIGNvdXQgPDwgIuCmuOCmsOCnjeCmrOCmvuCmp+Cmv+CmlSDgpqrgp4HgpqjgprDgpr7gpqzgp4PgpqTgp43gpqTgpr/gprAg4Kau4Kan4KeN4Kav4KeHIOCmquCmpuCnjeCmp+CmpOCmv+Cmn+CmvyDgprjgpqjgp43gpqjgpr/gpqzgp4fgprbgpr/gpqQg4Ka54Kav4Ka84Kao4Ka/4KWkXG4iOwogICAgfQoKICAgIHJldHVybiAwOwp9Cg==