
double k_noedge(               /* Calculates K(h) without edge correction */
		coord *points,       /* An array of positions */
		int n,               /* Number of points */
		double a,            /* Length of box in x-direction */
		double b,            /* Length of box in y-direction */
		double h);           /* Parameter of k */

double k_tor(              /* Calculates K(h) on a torus */
	     coord *points,       /* An array of positions */
	     int n,               /* Number of points */
	     double a,            /* Length of box in x-direction */
	     double b,            /* Length of box in y-direction */
	     double h);           /* Parameter of k */

double k_edge(               /* Calculates K(h) with edge correction */
	      coord *points,       /* An array of positions */
	      int n,               /* Number of points */
	      double a,            /* Length of box in x-direction */
	      double b,            /* Length of box in y-direction */
	      double h);           /* Parameter of k */

double k_rip(               /* Calculates Ripley's edge-correcting estimate 
			       for K(h) */
	     coord *points,       /* An array of positions */
	     int n,               /* Number of points */
	     double a,            /* Length of box in x-direction */
	     double b,            /* Length of box in y-direction */
	     double h);           /* Parameter of k */

double edge(coord *point,     /* Location of point */
	    double a,       /* size of rectangle in x-direction */
	    double b);      /* size of rectangle in y-direction */

double weight(coord *point,   /* Location of point */
	      double dist,    /* Radius of circle */
	      double a,       /* size of rectangle in x-direction */
	      double b);      /* size of rectangle in y-direction */
