@Preben Alsholm thank you so much.. actually i dont really understand the step you show me above..

did i need to find initial for f'',g' and theta' first then we can solved it normal way??

@Preben Alsholm thank you Preben :)

@Carl Love thank you carl.. its working :)

@Preben Alsholm im agree with u.. after figure out.. the equations of F is double prime so it need initial condition.. huhu..

@Preben Alsholm actually the F doesnt have initial condition.. before this i have tried the same problem and its ok..  but for this problem its doesnt run well.

@Carl Love here the code.. i cant upload the file.. dont know why.. if u want the file.. could u mind give ur email please.. :)

restart; with(plots); Pr := 1.0; B := .5; Ec := 2; l := 1; y := 1; Bt := 1; N := 10;
Eq1 := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta, eta))-(diff(f(eta), eta))^2-M*(diff(f(eta), eta))+l*B*(diff(F(eta), eta)-(diff(f(eta), eta))) = 0;

Eq2 := (diff(F(eta), eta))^2-F(eta)*(diff(F(eta), eta, eta))+B*(diff(F(eta), eta)-(diff(f(eta), eta))) = 0;
print(`output redirected...`); # input placeholder
                 
Eq3 := diff(theta(eta), eta, eta)+Pr*((diff(theta(eta), eta))*f(eta)-2*(diff(f(eta), eta))*theta(eta))+l*y*Bt*Pr*(theta1(eta)-theta(eta))+l*B*Pr*Ec*(diff(F(eta), eta)-(diff(f(eta), eta)))^2 = 0;
 
Eq4 := 2*(diff(F(eta), eta))*theta1(eta)-F(eta)*(diff(theta1(eta), eta))+Bt*(theta1(eta)-theta(eta)) = 0;

bcs1 := f(0) = 0, (D(f))(0) = 1, (D(f))(N) = 0, F(N) = f(N), theta(0) = 1, theta(N) = 0, theta1(N) = 0;
 
L := [0, 5, 10];

for k to 3 do R := dsolve(eval({Eq1, Eq2, Eq3, Eq4, bcs1}, M = L[k]), [f(eta), F(eta), theta(eta), theta1(eta)], numeric, output = listprocedure) end do;
%;
Error, (in dsolve/numeric/bvp/convertsys) too few boundary conditions: expected 8, got 7

@Carl Love thank you carl.. :) its working..

@Carl Love its work but it does not printout the graph but this

do you have any suggestion?

@Preben Alsholm Its ok .. Thank you Preben for your help .. :)

@Preben Alsholm its working thank you very much.. but when i change again the parameter.. like B=2.. it give the same error again..Error, (in dsolve/numeric/bvp) Newton iteration is not converging i try changing the value of N but its still not working.. may i know how you choosing the N ?

@J4James if you have any idea may i know why it not converge for certain valu such as if i change B=2 or 3.. or theta2=-5 or -1.5 .. it cannot work.. hydronew2.mw

@J4James thank yo so much :)

@mehdi jafari i solve the Equations using Keller Box method.. thats its second order finite difference(Newton mwthod).. since i still new in Maple i dont know much how much MAPLE can help me.. :)

@J4James may i know how do can i print the value of the graph.. actually i just i want to know the initial value of the graph.. how can i get the value? i write a comand (odeprint(eta,F(0)) but it doesnt come out the value.. i just learned to used MAPLE two days ago.. so many things i lack about it.. could you help me how can i get the value from the graph.. your coding before it really helpfull.. thank you :)

@mehdi jafari yup.. the boundary conditions are actually uncomplete.. so i do this computation to find its initial condition.. after i get the initial i will use matlab to compute the results.. I will try as your suggestion.. Thank you so much.. :)