function y = H(x,i)
y = x(i).^2.*sin(3*x(i));