计算方法实验7--牛顿下山法.cpp

#include “stdio.h“
#include “conio.h“
#include “math.h“
#include “stdlib.h“

#define Et 1e-3  /* 下山因子下界 */
#define E1 1e-3  /* 根的误差限 */
#define E2 1e-3   /* 残量精度 */

double f（double x） { return x*x*x - x - 1;}   /* f（x）的表达式 */
double f1（double x） { return 3*x*x - 1;}  /* f（x）的导数 */

void main（）
{
    int b;
    
    double f（double）;
    double f1（double）;
    
    void errormess（int）;   /* 根据计算结果，输出判断情况 */
    
    double x0;
    
    int Newton（double（*）（double）double（*）（double）double &x0）;
    
    x0 = 0.6;
    
    b = Newton（ff1x0）;
    
    if（ b = 1 ）   printf（“\n\nThe root of is x=%f\n\n“x0）;
      
    else errormess（b）;
}

int Newton（ double（*f）（double） double（*f1）（double） double &x0 ）
{
	int k = 0;
    double t = 1.0;
    double xk;
	printf（“k\t\t   t\t\t   xk\t\t  f（xk）\n“）;
	printf（“--------------------------------------------------------------------------------“）;
    printf（“%d\t\t\t\t%f\t%f\t\t\t“kx0f（x0））;
    
    xk = x0 - t*（  f（x0）/f1（x0） ）