#define MAXPREC 50

double rad(int rd,double num,double in,double fi,int prec)
{
 double m=(fi+in)/2;
 double test;
 if(num<0)return 0.0;
 if(prec>=MAXPREC)return m;
 if((test=pow(m,rd))==num)return m;
 if(test<num)return rad(rd,num,m,fi,prec+1);
 return rad(rd,num,in,m,prec+1);
 

}

double pow(double base,int exp)
{
	if((base==0)||(base==1)||(exp==1))return base;
	if(exp==0)return 1;
	return base * pow(base,exp-1);
}

/*Esempio di chiamata per la radice di due...
  
  rad(2,2.0,0,2,0);

rd   = La base della radice (es: 2 per la radice quadrata 3 per la cuba ecc.)
num  = Il radicando
in   = Inizio del intervallo di ricerca della radice
fi   = Fine del intervallo...
prec = Variabile sentinella per le ricorsioni
MAXPREC = Valore massimo di ricorsioni (nel caso di radici irrazzionali interrompe
la ricerca ad un valore approssimato)
*/