((Subscript[I, r])((Subscript[r, i]),(Subscript[θ, i]))=(Subscript[I, i])(A/(r^2))cos((Subscript[θ, i]))BRDF((Subscript[θ, i])).(3))