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