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]]
