Solo quería completar la respuesta dada por Siddhant Das usando las condiciones de contorno dadas.
Encontró una solución particular de la EDO en [math] f (\ theta) [/ math] obtenida de ansatz [math] \ psi (r, \ theta, \ phi) = f (\ theta) / r ^ 2 [ /matemáticas],
[matemáticas] f ” + f ‘\ cot \ theta + 2f = \ sin ^ 2 \ theta [/ matemáticas]
La solución particular que encontró fue
[matemática] f_p = \ frac {\ cos \ theta} {4} (\ cos \ theta – \ log \ cot \ theta / 2) [/ math]
Necesitamos obtener una solución que sea finita para todos [math] \ theta [/ math] para que coincida con la condición de límite que [math] \ psi [/ math] va a cero en el infinito. Hay dos constantes arbitrarias de la solución homogénea,
[matemáticas] f_h = c_1 \ cos \ theta + c_2 (\ cos \ theta \ log \ cot \ theta / 2 – 1) [/ matemáticas]
Evidentemente, podemos cancelar la parte infinita de [math] f_p [/ math] estableciendo [math] c_2 = – \ frac {1} {4} [/ math], que da
[matemáticas] f = \ frac {1} {4} (1 + \ cos ^ 2 \ theta) + c_1 \ cos \ theta [/ matemáticas]
La constante restante [matemática] c_1 [/ matemática] no se puede determinar a partir de las condiciones de contorno ya que [matemática] \ frac {\ cos \ theta} {r ^ 2} [/ matemática] es armónica y se define en todas partes, excepto el origen. Sin embargo, en una situación física ( es decir, donde [math] \ frac {\ sin ^ 2 \ theta} {r ^ 4} [/ math] representa una distribución de carga) debemos elegir el valor de [math] c_1 [/ math ] para que el potencial tenga la misma simetría que la fuente. La fuente es simétrica con respecto a una reflexión en el plano xy; la solución será si y solo si elegimos [math] c_1 = 0 [/ math] ya que el término [math] \ cos \ theta [/ math] es impar con respecto a esta reflexión y los otros términos son pares. Esto da la respuesta final,
[matemáticas] \ psi (r, \ theta, \ phi) = \ frac {1 + \ cos ^ 2 \ theta} {4r ^ 2} [/ matemáticas]