"Not sure what I'm doing wrong here."
You need to use element-wise division, not matrix division:
phi = 0:0.1:6*pi;
I2 = sin(phi).^2 ./ sin(phi/2).^2;
% ^^
plot(I2)
Tip: if you are not doing linear algebra you probably should be using array operations:
