Boundary value problem with bvp4c

Ive attached my bvp4c code. I want to genare the figure im showing here, but im am not getting the desired graph .

I am attaching both desired graph and coupled equations with boundary condition.

I would be very helpfull for your help.

Answer 1

[Ameer Hamza]

There were issues in your ODE function and the BC function. Replace them with these defintions.

function res = bc(ya,yb)
res = [ ya(1); ya(2)-1; ya(5)-1; yb(2); yb(3); yb(5)];
end

function dydx =jeffrey_fluid(x,y)
M=1;ga=5;la=1;beta=2;Pr=0.7;

dydx = [y(2); 
        y(3); 
        y(4);
        -1/y(1)*(-1/beta*((1+la)*(y(1)*y(3)-y(2)^2-y(2)*(ga+M)+la*y(5))-y(4))-y(3).^2); 
        y(6); 
        -Pr*(y(1)*y(6)-y(2)*y(5))];
end

[Didarul Ahasan Redwan]

thanks a lot ...<3

[Ameer Hamza]

Glad to be of help.