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