Question: Error in numeric dsolve

Hello,

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

But I received errors.

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.

Please Wait...