Question: Graph of system of ODEs

Aslam o Alikum!

Hope everyone going fine with sound health I want to make the graph of following system of ODEs

Eq1 := diff(F(eta), eta, eta, eta, eta)-M*(eta*(diff(F(eta), eta, eta, eta))+3*(diff(F(eta), eta, eta))+(diff(F(eta), eta))*(diff(F(eta), eta, eta))-F(eta)*(diff(F(eta), eta, eta, eta)))-Ha^2*(diff(F(eta), eta, eta));
Eq2 := diff(G(eta), eta, eta)+Pr*M*(F(eta)*(diff(G(eta), eta))-eta*(diff(G(eta), eta)))+Pr*Ec*(diff(F(eta), eta, eta))^2+Nb*(diff(H(eta), eta))+diff(G(eta), eta)+Nt*(diff(G(eta), eta))^2;
Eq3 := diff(H(eta), eta, eta)+M*Sc*(F(eta)*(diff(H(eta), eta))-eta*(diff(H(eta), eta)))+Nt*(diff(F(eta), eta, eta))/Nb;
 

using the BC's

IC1 := F(0) = 0, ((D@@2)(F))(0) = 0, (D(G))(0) = 0, (D(H))(0) = 0;
     F(0) = 0, @@(D, 2)(F)(0) = 0, D(G)(0) = 0, D(H)(0) = 0
IC2 := F(1) = 1, (D(F))(1) = 0, G(1) = 1, H(1) = 1;
 

for Ec := .1; Nt := .1; Nb := .1; Sc := .5; Pr := 10; M := .5; and different values of Ha like Ha:=0, Ha:=2, Ha:=4, Ha:=6 and Ha:=8 in one coordinate. 

I am waiting your positive response.

Please Wait...