Please, i need help USING ODE1 and ODE2 with given BCS and Pr=0.714

it is needed to to generates

[-0.2], [0.51553], [0.4000]

[-0.1], [0.57000], [0.4371]

[0.], [0.62756], [0.4764]

[0.1], [0.68811], [0.5176]

[0.2], [0.75153], [0.5609]

but it is generarting

[-0.2], [0.51553], [0.42342]

[-0.1], [0.57000], [0.46114]

[0.], [0.62756], [0.50088]

[0.1], [0.68811], [0.54261]

[0.2], [0.75153], [0.58628]

the values of D(theata)(0) is wrong. Please i need HELP. this the code below that i use:

>restart;

>with (plots):ode1:=diff(f(eta),eta,eta,eta)+f(eta)*diff(f(eta),eta,eta)-M*diff(f(eta),eta)=0:

>ode2:=diff(theta(eta),eta,eta)+Pr*f(eta)*diff(theta(eta),eta)=0:

>bcs1:= f(0)=w,D(f)(0)=1,D(f)(10)=0:

>bcs2theta(10)=0,theta(0)=1:

>fixedparameter1:=[M=0.0]:

>ode3:=eval(ode1,fixedparameter1):

>fixedparameter2:=[Pr=0.714]:

>ode4:=eval(ode2,fixedparameter2):

>G:=[-0.2,-0.1,0.0,0.1,0.2]:

>for ode3 and ode4:

for k from 1 to 5 do

sol_All:=dsolve(eval({ode3,ode4,bcs1},w=G[k]), [f(eta),theta(eta)],numeric,output=listprocedure);

Y_sol||k:= -rhs(sol_All[4]);

YP_sol||k:=-rhs(sol_All[6]);

end do:

>Digits:=5:

>for k from 1 to 5 do

evalf([G[k]]),evalf([(Y_sol||k(0))]),evalf([YP_sol||k(0)]);

od;

[-0.2], [0.51553], [0.42342]

[-0.1], [0.57000], [0.46114]

[0.], [0.62756], [0.50088]

[0.1], [0.68811], [0.54261]

[0.2], [0.75153], [0.58628]