*Ts,double *Tinitial,double *counter1);
void meanvar(double *pmatrix,int *g,int *nr,int *nc,double *storage);
void start(double *Imatrix,int *g,int *nr,int *nc,double *testing,
*Ts,double *Tinitial,double *counter1);
void meanvar(double *pmatrix,int *g,int *nr,int *nc,double *storage);
void start(double *Imatrix,int *g,int *nr,int *nc,double *testing,