MATLAB: How to solve nonlinear coupled ode by Shooting method .

Respected sir ,

I am getting problem to solve non linear coupled BVP by shooting method . Can you help me to solve that problem?

My problem is

F'''=(1/(1+epsilon1*(1-G)))*(A*F'+.5*t*A*F''+(F')^2-F*F''+epsilon1*G'*F''-lambda*G-delta*H+M*F'-(-1+epsilon1*G-epsilon1)*bita*F');

G''=(1/(1+epsilon2*G+Nr))*(Pr*(2*A*G+.5*A*t*G'+F'*G-F*G'-M*Ec*(F')^2-Ec*(1-epsilon1*G'+epsilon1)*(F'')^2)-epsilon2*(G')^2);

H''= Sc*(2*A*H+.5*A*t*H'+F'*H-F*H'+Rex*Zai*H);

where A=0 ; epsilon1= 0 ; epsilon2=1; lambda= 1; delta=1; bita=0; Nr=.1; Pr=5; M=.5;

Ec=.1 ; Sc=1; Rex= .3 ; Zai= .1 ;

and F(0)=0 , F'(0)=1 F'(infity)=0 G(0)=1 G(infity)=0 H(0)= 1 H(infity)=0

and F''(0) , G'(0) & H'(0) we have to guess

so tell me how to solve by shooting method with using rk -4 method .