الگوريتم روش نيتون در محاسبه ريشه (درس محاسبات عددي)

[كیانی]

کد:

#include<iostream.h>
#include<math.h>
#include<conio.h>
double f( double x ) { return 100 * sin(x) - x;}
double fprime( double x ) { return 100 * cos(x) - 1;}
void main(){
clrscr();
double eps,x0,x;
int n=1;
cout<< "Enter x0,eps:";
cin>> x0 >> eps;
x = x0 - f(x0) / fprime(x0);
while ( fabs(f(x) ) >= eps){
cout<< x << "\n";
x0 = x;
x = x0 - f(x0) / fprime(x0);
n++;
}
cout<< "ROOT = "<< x;
cout<< "\nITERATION = "<< n;
getch();
}