/*trial.h*/
#ifndef _EXT_VARS_H_
	#define _EXT_VARS_H_
	#include "trial.h"
#endif

void start_upcrmv(DVECTOR wdose1,DVECTOR wdose2,DVECTOR parm, DVECTOR priori_alpha,int *level1,int *level2,int n_r,int n_c);
double hv(double dose1, double dose2, DVECTOR parm); /*start_upcrm*/
void prep_crmv(DVECTOR crm_parm,DMATRIX toxicp,IVECTOR y,IVECTOR y_lev1,IVECTOR y_lev2,int tot_pid,int n_r,int n_c,double *r);
void get_crmparm(IVECTOR y, IVECTOR level1,IVECTOR level2,DVECTOR wdose1, DVECTOR wdose2,int pi,DVECTOR priori_alpha,DMATRIX toxicp, int n_r, int n_c);
void compute_muparm(long double f (IVECTOR,IVECTOR,IVECTOR,DVECTOR,DVECTOR,int,DVECTOR),IVECTOR y, \
IVECTOR level1,IVECTOR level2,DVECTOR wdose1,DVECTOR wdose2,int pi,DVECTOR priori_alpha,DMATRIX toxicp, int n_r, int n_c);
long double prod_phi(IVECTOR y,IVECTOR level1,IVECTOR level2,DVECTOR wdose1,DVECTOR wdose2,int pi, DVECTOR parm);
double phi(int res, double dose1,double dose2, DVECTOR parm);
void recommand_dosecrmv(IVECTOR y,IVECTOR y_lev1,IVECTOR y_lev2,int tot_pid,DVECTOR crm_parm,DVECTOR mtd,DVECTOR bias,DVECTOR mse,int *count,int n_r,int n_c);
double atanh(double x);