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
