Question: Problem with shooting method

Dear all,

I would like to ask you to help me with the following problem that I got error.

Thanks in advance

 

restart; Shootlib := "C:/Shoot9"; libname := Shootlib, libname; with(Shoot);
with(plots);
N1 := 1; N3 := .5; N2 := 5; Bt := 6; AA := N1*Bt; gamma1 := .2; blt1 := 1;
FNS := {f(eta), fp(eta), fpp(eta), g(eta), gp(eta), h(eta), hp(eta), i(eta), ip(eta), fppp(eta)};
ODE := {diff(f(eta), eta) = fp(eta), diff(fp(eta), eta) = fpp(eta), diff(fpp(eta), eta) = fppp(eta), diff(g(eta), eta) = gp(eta), diff(gp(eta), eta) = N1*(2*g(eta)+eta*gp(eta)+2*g(eta)*fp(eta)-2*f(eta)*gp(eta)+2*N2*N3*(h(eta)*ip(eta)-i(eta)*hp(eta))), diff(h(eta), eta) = hp(eta), diff(hp(eta), eta) = AA*(h(eta)+eta*hp(eta)-2*f(eta)*hp(eta)+2*h(eta)*fp(eta)), diff(i(eta), eta) = ip(eta), diff(ip(eta), eta) = AA*(2*i(eta)+eta*ip(eta)-2*f(eta)*ip(eta)+2*N2*h(eta)*gp(eta)/N3), diff(fppp(eta), eta) = N1*(3*fpp(eta)+(eta-2*f(eta))*fppp(eta)-(2*N2*N2)*(diff(hp(eta), eta)))};
IC := {f(0) = 0, fp(0) = gamma1*fpp(0), g(0) = 1+gamma1*gp(0), gp(0) = beta, h(0) = 0, hp(0) = beta1, i(0) = 0, ip(0) = beta2, fppp(0) = alpha};
BC1 := {f(blt1) = .5, fp(blt1) = gamma1*fpp(blt1), g(blt1) = gamma1*gp(blt1), h(blt1) = 1, i(blt1) = 1};
infolevel[shoot] := 1;
S := shoot(ODE, IC, BC1, FNS, [alpha = .1, beta = .2, beta1 = .3, beta2 = .4], maxfun = 50000);
 
Please Wait...