Hi there, I am ultimately trying to solve a system of differential equation numerically that reads

diff(y(x),x,x) + (2/x)diff(y(x),x) + y(x)(1-(y(x))^2) = (F(x) - G(x)/x^2)
(1/2 - Xi(x))F(x) = -diff(G(x),x) - (1/x)G(x)
(1/2 + Xi(x))G(x) = diff(F(x),x) - (1/x)F(x)