function [r2,s]=t2(r1)
 
r2 = r1+0.01;
n=length(r2);
s=prod(r2).^(1./n);
