#include "k.h"
/* 	Make data, a vector of root(chi-squares) length n             
*	Input n,R,S,gamma,F1,F2,oo, x (for results, length n) and xt (length n) workspace                          
*	Output x (length n) */
void datagen(int n,int R,int S,double gamma,double F1,double F2,int oo,double *x,double *xt,long *idum) {
int k,o[3],e[3],N;
N=R+S;
for (k=0;k<(n-oo);k++) {
	datnull(R,F1,idum,o);
	datnull(S,F2,idum,e);
	xt[k]=chi(o,e,R,N);
}
for (k=(n-oo);k<n;k++) {
	datdis(R,gamma,idum,o);
	datnull(S,F2,idum,e);
	xt[k]=chi(o,e,R,N);
}
k=n-oo;
sort(xt,&k,x);
sort(xt+k,&oo,x+k);
return;
}