Sorry

 endv=100 not 5 (as shown below)

I have used dsolve to obtain the solution of the inserted two systems below. The obtained results for the first system is accurate, where the obtaind solution of the second one is not accurate. I think the problem is in approximating the first system solution to be used in the second one. I would appreciate any advice on how to  obtain accurate results for the second system

I have used dsolve to obtain the solution of the inserted two systems below. The obtained results for the first system is accurate, where the obtaind solution of the second one is not accurate. I think the problem is in approximating the first system solution to be used in the second one. I would appreciate any advice on how to  obtain accurate results for the second system

 restart; endv := 5; sys0 := diff(f(x), `$`(x, 2...

I used dsolve to obtain the solution of the attached two systems of ODEs, where the second system depends on the solution of the first one. The obtained results for the first system is accurate, where the obtaind solution of the second one is not accurate. I think the problem is using cubic spline approximation. I would appreciate any advice on how to  obtain accurate results for the second system

I have solved  the following BVP system in Maple

restart; e := 12; Sc := .1; Pr := .7; Gr := 1; M := .1; sys := diff(f(t), `$`(t, 3))+2*f(t)*(diff(f(t), `$`(t, 2)))-(diff(f(t), t))^2-M^2*(diff(f(t), t))+Gr*theta(t) = 0, diff(theta(t), `$`(t, 2))+2*Pr*f(t)*(diff(theta(t), t)) = 0, diff(phi(t), `$`(t, 2))+Sc*(2*f(t)*(diff(phi(t), t))-phi(t)) = 0; bcs := f(0) = .1, (D(f))(0) = 0, theta(0) = .5, phi(0) = 1, (D(f))(e) = 0, theta(e) = 0, phi(e) = 0; sol := dsolve({bcs, sys}, numeric); plots[odeplot...