@Carl Love 

Thanks for the reply and noticing. It seems I have been to mix up between bvp and rkf. 

when I remove Eq3, there are other errors come up

Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging

Meanwhile, when I remove Eq2, the error come up

Error, (in dsolve/numeric/bvp/convertsys) too many boundary conditions: expected 7, got 8

Do I have issue on ode itself?