function(x, y, alfa)
{
	z1 <- (x - trim.mean(x, alfa))/
		trim.std(x, alfa)
	z1 <- z1 + (y - trim.mean(y, alfa))/
		trim.std(y, alfa)
	z2 <- (x - trim.mean(x, alfa))/
		trim.std(x, alfa)
	z2 <- z2 - (y - trim.mean(y, alfa))/
		trim.std(y, alfa)
	v1 <- (trim.std(z1, alfa))^2
	v2 <- (trim.std(z2, alfa))^2
	rs <- (v1 - v2)/(v1 + v2)
	rs
}
