Come suggerisce il titolo, il modo migliore per affrontare questo problema è il metodo delle cariche immagine.
Fisso un sistema di assi cartesiani

tale che il piano conduttore abbia equazione

e che la carica

si trovi in
)
. Grazie ai teoremi di unicità dell'elettrostatica, se riesco a trovare una distribuzione di cariche (fittizie) tale che il potenziale
)
dovuto a esse e alla carica

soddisfi la condizione

sull'intera superficie del conduttore (in quanto è collegato a terra), allora il potenziale elettrico nella parte di spazio vuota sarà proprio
)
.
Introduco pertanto tre cariche immagine:

in
)
,

in
)
e

in
)
. Scrivo l'espressione del potenziale:
=\frac{1}{4 \pi \epsilon_0}\bigg[\frac{q}{\sqrt{x^2+y^2+(z-d)^2}}-\frac{q}{\sqrt{x^2+y^2+(z+d)^2}}+\frac{qa/d}{\sqrt{x^2+y^2+(z+a^2/d)^2}}-\frac{qa/d}{\sqrt{x^2+y^2+(z-a^2/d)^2}}\bigg])
Per

, la superfice del conduttore è descritta da

,

, mentre per

l'equazione della superficie è

. Inserendo la prima condizione in
)
, il primo termine si cancella con il quarto e il secondo col terzo, inserendo la seconda invece il primo termine si elide col secondo e il terzo col quarto, dunque

,

e

costituiscono la distribuzione cercata.
L'espressione
)
ovviamente vale solo per il potenziale elettrico sulla superficie del conduttore e nello spazio ad esso esterno, in quanto all'interno è noto che

.
La densità di carica indotta avrà modulo massimo nel punto del conduttore più vicino a

, cioè
)
. Per trovarla, ricorro al Teorema di Coulomb:

Dove

è la densità superficiale di carica,

è la differenza fra il campo elettrico appena fuori e appena dentro il conduttore, e

è il versore normale alla superficie e uscente.
In

è

. Il campo elettrico appena dentro il conduttore è identicamente nullo, mentre quello appena fuori si trova facilmente conoscendo la posizione delle quattro cariche oppure applicando l'operatore gradiente a
)
:
=-\frac{q(3d^2+a^2) \hat z}{2\pi \epsilon_0 (d^2-a^2)^2})
Quindi:
 \cdot \hat z =-\frac{q(3d^2+a^2)}{2 \pi (d^2-a^2)^2})
Riporto anche, per chi fosse interessato, il mio risultato per

su tutta la superficie del conduttore:
 =
\Bigg \{
\begin{array}{rl}
-\frac{q(d^2-a^2)}{4 \pi a} \big[(a^2+d^2-2dz)^{-3/2}-(a^2+d^2+2dz)^{-3/2}\big] & (x^2+y^2 \leq a^2) \\
-\frac{qd}{2 \pi} \big [(x^2+y^2+d^2)^{-3/2}-a^3 (d^2x^2 +d^2 y^2 +a^4)^{-3/2} \big] & (x^2+y^2 \geq a^2) \\
\end{array})