>> theta=-pi:0.01:pi;
>> y=f(theta);
>> plot(theta,y);
>> grid

>> format long
>> fzero(@f,1.5)

ans =

  1.547945496012254

>> fzero(@f,2.3)

ans =

  2.298884445898739