function [r1,e]=t1(r0)
 
r1 = r0+0.1;
r1 = t2(r1);
n=length(r1);
e=sum(r1)./n;
