#include <stdio.h>
#include <math.h>
#include "protos.h"

#define epsilon 1.e-6 /* precisao */
#define PI 3.14159

void funcd(float x, float *f, float *df){

  *f = cos(x)-x;
  *df = -sin(x)-1.;

}

int main() {

  float x1, x2; /* extremos do intervalo de busca */
  float raiz;   /* armazena raiz procurada */

  /* inicializa intervalor de busca */
  x1 = 0.;
  x2 = 0.5*PI;

  /* Chamada da rotina rtnewt */
  raiz = rtnewt(x1, x2, epsilon);

  /* Imprime resultado na tela se a raiz foi encontrada */
  if (raiz != 0.0)
    printf("Raiz de f(x): %10.8f\n", raiz);

  return 0;

}