NAME:
  TXY

 PURPOSE:
  Compute the synchronziation measure T(x|y).
  For details of the measure, please see
  M. A. Kramer, E. Edwards, M. Soltani, M. S. Berger, R. T. Knight and A. J. Szeri
  "Synchronization measures of bursting data: Application to the electrocorticogram of an auditory event-related experiment."
  Physical Review E , 70: 011914 (2004).

 CALLING SEQUENCE:
  res = txy(ensemble_1, ensemble_2, embedding_dimension, delay_time, eps_radius_factor)

 INPUTS:
    ensemble_1:  An array constituting the first ensemble of time series, e.g. s[# of ensemble members, length of each ensemble member].
    ensemble_2:  An array constituting the second ensemble of time series, e.g., r[# of ensemble member, length of each ensemble member].
    embedding_dimension:  An integer specifying the embedding dimension for the data.
    delay_time:  An integer specifying the delay time [in indices!!] for the data.
    eps_radius_factor:  The amount by which to scale the epsilon neighborhood of a point.

 DEPENDENT ROUTINES:
    embed_data.pro --> embed an ensemble of data given the embedding dimension and delay time.

 EXAMPLE:
    Please see the file txy_example.pro.

 MODIFICATION HISTORY:
    MAK (BU).  November 29, 2006.