Como suponía y no desmientes, la función
q(x^2+y^2) es una función de dos variables, llemémosla f(x, y)
f(x,y) = q(x^2+y^2)
Calculamos las derivadas parciales segundas que se necesitan para estudiar si es una función armónica.
fx(x,y) = q'(x^2+y^2)·2x
fxx(x,y) = q''(x^2+y^2)·2x·2x+q'(x^2+y^2)·2 = 4x^2·q''(x^2+y^2) + 2q'(x^2+y^2)
Análogamente
fyy(x,y) = 4y^2·q''(x^2+y^2)+2q'(x^2+y^2)
La suma de ambas debe ser cero
4(x^2+y^2)q''(x^2+y^2)+4q'(x^2+y^2) = 0
(x^2+y^2)q''(x^2+y^2)+q'(x^2+y^2) = 0
Que haciendo el cambio de variable z=x^2+y^2 nos da esta ecuación diferencial con q la función y z la variable independiente
zq'' + q' = 0
hacemos el cambio de variable p = q' con lo cual quedará
zp' + p = 0
zp' = -p
p' / p = -1/z
Integrando en ambos lados
ln p = -ln z + ln C
ln p = ln(C/z)
p = C/z
q = $(C/z) dz = Cln(z) + C2
q(z) = C·ln(z) + C2
Asi debe ser la función q, con dos constantes C y C2 cualesquiera de R.
Como no da gusto usar subíndices con este editor llamare A y B a las constantes
Luego la función armónica será
f(x,y) = A·ln(x^2+y^2) + B
Lo dejo aquí de momento, esto no lo he dado y me cuesta. Haz el favor de puntuar el trabajo y manda el otro apartado en una pregunta nueva.