function h = parl(a,b)
	ab = a.*b;
	h = ab./(a + b);