a=load('test4b.in');
[sum,av]=mysum(a)

