How can I solve simultaneously two nonlinear coupled BVP

F''' + F*F' - (F')*(F') + 1 - M*F' = 0

G'' + P*F*G'= 0

Boundary condition:

F(0) = 0, F'(0) = 0, F'(inf) = 1

G(0) = 1, G(inf) = 0 or

G'(0) = -Bi(1 - G(0)), G(inf) = 0

where N , Bi and P are constant

I was able to solve the first equation by using the following code, but I do not no how to modify it to solve the two equations simultaneously.

function F = eqn1(t,y)

M = 0.1;

F = zeros(3,1);

F(1) = y(2);

F(2) = y(3);

F(3) = y(2)*y(2) + M*y(2) - y(1)*y(3) - 1;

function res = bc(ya,yb)

res = [ya(1);ya(2);yb(2)-1];

solt = bvpinit([0,5],[0;0;0;]);

Ans = bvp4c(@eqn1,@bc,solt);

X = Ans.x;

Y = Ans.y;

plot(X,Y)

I awaits kind responses Thanks