Question: dsolve bvp systems

I solved the following system using fourth order finite difference and shooting methed and i obtained the results below

but when i tried to solve it using dsolve in Maple I did not obtain the accuract results, I hope to know how to obtain accurate results with this system(the system is uploaded)

 

 restart; Pr := .72; a := 1; c := 1; lambda := 1; M := .2; lambda1 := .2; beta := 0; Ec := 0; A := 0; B := 0; R := 0;

sys := diff(f(t), `$`(t, 3))-(diff(f(t), t))^2+f(t)*(diff(f(t), `$`(t, 2)))+(M^2+lambda1)*(a/c-(diff(f(t), t)))+a^2/c^2+lambda*theta(t)-beta*(2*(diff(f(t), t))*(diff(f(t), `$`(t, 3)))-f(t)*(diff(f(t), `$`(t, 4)))-(diff(f(t), `$`(t, 2)))^2) = 0, (1+(4/3)*R)*(diff(theta(t), `$`(t, 2)))+A*exp(-t)+B*theta(t)+Pr*(f(t)*(diff(theta(t), t))-(diff(f(t), t))*theta(t)+M^2*Ec*(diff(f(t), t))*(diff(f(t), t)-a/c)+Ec*(diff(f(t), `$`(t, 2)))) = 0;

bcs := f(0) = 0, (D(f))(0) = 0, theta(0) = 1, (D(f))(10) = a/c, theta(10) = 0; Digits := 10; Digits := 10; sol := dsolve({bcs, sys}, numeric, abserr = 1.*10^(-9)); sol(0);

shooting and finite difference results   f"(0)=0.3559

dsolve results  f"(0)=1.7765

 

Please Wait...