function z=hx(x) %usage: z=hx(x); z=sum(x.^2);