We determine the range function by determining the x coordinate for the time t (Solve[ y[t,theta]==0,t]) at which the projectile strikes the ground.

Input := 


range[theta_] =  x[t,theta]/.

	Solve[ y[t,theta]==0,t] [[2]]

Output =


5625 Cos[theta] Sin[theta]