We plot the range function to estimate the value of theta which gives a maximum range.

Input := 


Plot[range[theta],{theta,0, Pi/2}]

Output =


-Graphics-