/* Copyright (C) 1994 by Thomas Glen Smith.	All Rights Reserved. */
/* binomx APL2 V1.0.0 **************************************************
* Binomial for complex numbers = (!n)%(!m)X!n_m.                       *
***********************************************************************/
#define INCLUDES STDIO
#include "includes.h"
void binomx(m,n,ret)
double *m,*n,*ret;
{
	Dividex; Gammax; Minusx; Timesx;
	double wrka[2],wrkb[2],wrkc[2];

	minusx(n,m,wrka);		/* wrka == n_m */
	gammax(wrka,wrkb,1);	/* wrkb == !n_m */
	gammax(m,wrka,1);		/* wrka == (!m) */
	timesx(wrka,wrkb,wrkc);	/* wrkc == (!m)X!n_m */
	gammax(n,wrka,1);		/* wrka == (!n) */
	dividex(wrka,wrkc,ret);	/* ret  == (!n)%(!m)X!n_m */
}