function [arith, geom] = twomeans(a,b) 
arith = (a + b)./2; %! Arithmetic mean is $\frac{a+b}{2}$
geom  = sqrt( a.*b ); %! Geometric mean is $\sqrt{a b}$
