restart;

Pr:=0.71: n:=-1:

eta0:=0.0699;

EQ1:=diff(H(x), x ) - x*diff(F(x), x ) ;

EQ2:=(1+x^2)*diff(F(x), x$2) + (3*x + x*F(x)-H(x))*diff(F(x), x) + F(x)^2 + G(x)^2 +2*P(x) + x*diff(P(x), x) ;

EQ3:=(1+x^2)*diff(G(x), x$2) + (3*x + x*F(x)-H(x))*diff(G(x), x) ;

EQ4:=(1+x^2)*diff(H(x), x$2) + (3*x + x*F(x)-H(x))*diff(H(x), x) + (1+F(x))*H(x)- diff(P(x), x);

EQ5:=(1+x^2)*diff(theta(x), x$2) + x*(1-2*n)*diff(theta(x), x) + n^2*theta(x) - Pr*( n*F(x)*theta(x) + ( H(x)-x*F(x) )*diff(theta(x), x) ) ;

EQ:={EQ1=0, EQ2=0,EQ3=0,EQ4=0 ,EQ5=0}:

IC:={ F(0)=0, G(0)=12, H(0)=0, theta(0)= 1, F(eta0)=0, G(eta0)=12, H(eta0)=0, theta(eta0)= 0, P(0)=0};

sol:= dsolve(EQ union IC,numeric,output=Array([0,0.0699]));