(1) Hu.m function a=hu(x,y,m) S = zeros(1,2*m+1); T = zeros(m+1,1); for k = 1:2*m+1 S(k) = sum(x.^(k-