牛顿迭代法 POJ 2868

数值计算刚学了牛顿迭代法求方程零点，今天正好看到POJ有关于这方面的题

这题水过·，关键要知道公式

g(x) = x – (x – A/x) / 2 根据这条公式迭代就可以

终止条件是前后两次迭代出来的值少于1e-6

#include <iostream>
#include <iomanip>
using namespace std;

int main()
{
int n;
double x1, x2;
int count;
double err;
//freopen("C:\\Users\\Haojian\\Desktop\\test.txt", "r", stdin);
while (cin >> n)
{
x1 = 1;
count = 0;

do{
x2 = x1 - (x1 - n/x1)/2;
count++;
err = abs(x2 - x1);
x1 = x2;
}while (err > 1e-6);

cout << count << ' '
<< setprecision(2)
<< setiosflags(ios::fixed)
<< x2 << endl;
}
return 0;
}