//---------------------------------------------------------------------------
#include <vcl.h>
#pragma hdrstop

#include "mathengine.h"
#include <math.h>

//---------------------------------------------------------------------------
#pragma package(smart_init)
//---------------------------------------------------------------------------
double __stdcall Distance(POINT &p1, POINT &p2){
  int REPEAT=8;
  int x = p1.x-p2.x;
  int y = p1.y-p2.y;
  //return sqrt(x*x + y*y);

  double a = x*x + y*y;      //    from "A Book On C", chapter 4, p198, excercise 36
  double r=0.5*(1.0+a);        //       maybe speed-up........
  for(int i=0;i<REPEAT;i++){     //
    r = 0.5*(r+a/r);               //
  }                                  //
  return r;                            //
}
//---------------------------------------------------------------------------