Question:Error in numeric dsolve

Question:Error in numeric dsolve

Maple

Hello,

I want to solve three coupled differential equations with initial and boundary conditions numerically and get the plots of solutions.

Could you please help me to solve the error and get all three plots that I need?

My codes:

restart;

sys := {diff(phi(eta), eta\$2) + 5.261282735*f(eta)*diff(phi(eta), eta) - 2.630641368*phi(eta) = 0, 1.059704409*diff(theta(eta), eta\$2) + 6.176017503*f(eta)*diff(theta(eta), eta) + 21.03607964*diff(f(eta), eta\$2) + 0.5*phi(eta) = 0, diff(f(eta), eta\$4) - 1.052256547*diff(f(eta), eta)*diff(f(eta), eta\$2) + 1.052256547*f(eta)*diff(f(eta), eta\$3) + 5.165076420*diff(theta(eta), eta) + 5.261282735*diff(phi(eta), eta) = 0, eval(diff(phi(eta), eta), {eta = 0}) = 1 + 0.5*eval(diff(f(eta), eta\$2), {eta = 0}), eval(diff(phi(eta), eta), {eta = 1}) = 0.5*eval(diff(f(eta), eta\$2), {eta = 1}), f(0) = -0.5, f(1) = 0.5, phi(0) = 1, phi(1) = 0, theta(0) = 1, theta(1) = 0};

dsol:=dsolve(sys,numeric);
Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system
plots[odeplot]((dsol),eta=0..1);
Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

Thank you.

﻿