/* --------------------------------- */
/* save as c16a.c                    */
/* --------------------------------- */
#include "x_hfile.h"
#include      "fa.h"
/* --------------------------------- */
int main(void)
{

double  a = .5;
double  b = .5;
pt2d    p = {a,b};
int     n = 4;
double  h = .1;

 clrscrn();
 p = i_pt2d(a,b); /* For info */
 printf(" Use Newton's method to approximate,        \n");
 printf(" the solutions of the following system :\n\n\n");

 printf("  | %s = 0    \n", feq);
 printf("  | %s = 0\n\n\n", geq);
 
 printf(" As a first approximation x = %.1f y = %.1f \n\n", a, b);
 
 stop();
 
  clrscrn();       
  p_newton_fxy(  n, f, g, h, p);   
  stop();        
   

  clrscrn();
  p = newton_fxy(  n, f, g, h, p); 
  printf(" the solutions of the following system is :\n\n\n");
  printf("         x = %f  y = %f  \n\n\n",p.x,p.y);

  printf(" f(%f,%f) = %f  \n",p.x,p.y, f(p.x, p.y));
  printf(" g(%f,%f) = %f\n\n",p.x,p.y, g(p.x,p.y) );
  stop();

 return 0;
}