

function y = f(theta)

     y = sin(4*theta).^2;
