function y = f(theta)

% y = 1 + sin(theta);
% y = 2 + 0.5*sin(theta) - 0.1*cos(theta) + (sin(theta)).^2 + cos(4*theta);
% y = cos(theta) + 4 * sin(theta).^3;
y = exp(-(theta - pi).^2) - 3*exp(-2*(theta - 3*pi/2).^2);

