I've calculated as follows:

 

matusita<-function(f1,f2,lw,up){

 

 fmatu<-function(x) ( sqrt(f1(x)) - sqrt(f2(x)) )^2

 

 sqrt (integrate( fmatu, lower=lw, upper=up, subdivisions=1000000)$value)

 

}

 

 

 

Josué Almansa 
Unitat de Recerca en Serveis Sanitaris 
Institut Municipal d'Investigació Mèdica, IMIM 
c/ Doctor Aiguader, 80 
08003 Barcelona 
e-mail: jalmansa@imim.es 

 


	[[alternative HTML version deleted]]


