function [y,acc]=hz_average1(u,acc)% returns the scalar value for % avaergae of matrix uT=numerictype(1,31,15);%N=size(u,1)*size(u,2);for i=1:Nacc(1) = acc(1) + u(i);end;;y=divide(T,acc,N);