p((Subscript[r, i])\[VerticalSeparator]θ)=(1/(Sqrt[2 π (Subscript[σ, i]^2)]))(e^((-(((d((Subscript[a, i]),θ)-(Subscript[r, i])))^2))/(2(Subscript[σ, i]^2))))