Question: time dependent variable

> restart;
> libname = [shootlib, libname];
> with(shoot);
Error, invalid input: with expects its 1st argument, pname, to be of type {`module`, package}, but received shoot
> with(plots);
Pr := 10; s = -.1; lambda := 0; Gr := 1.0; Gm := 1.0; beta := -1.20;
10
s = -0.1
0
1.0
1.0
-1.20
> M := 0.; z := .1; Xi := .5; Nt := .5; Nb := .2; l := 5; Nr := .5; epsilon1 := .2; epsilon2 := .2;
0.
0.1
0.5
0.5
0.2
5
0.5
0.2
0.2
> Prff := Pr/(1+4.*N*(1/3));
10
-----------------
1 + 1.333333333 N
> FNS := {f(eta), h(eta), r(eta), u(eta), v(eta), theta(eta), `ϕ`(eta)};
{f(eta), h(eta), r(eta), u(eta), v(eta), theta(eta), ϕ(eta)}
> ODE := {diff(h(eta), eta)+.75*l*f(eta)*h(eta)-(1/4)*l*u(eta)*epsilon2-Nt*(.75*f(eta)*r(eta)-(1/4)*u(eta)*epsilon1+Nb*r(eta)*h(eta)+Nt*r(eta)*r(eta))/Nb = 0, .75*f(eta)*r(eta)+diff(r(eta), eta)-(1/4)*u(eta)*epsilon1+Nb*r(eta)*h(eta)+Nt*r(eta)*r(eta) = 0, diff(v(eta), eta)+3*(f(eta)*v(eta)-u(eta)*u(eta))/(4*Pr)-(M+lambda)*u(eta)+theta(eta)-Nr*`ϕ`(eta) = 0, diff(f(eta), eta) = u(eta), diff(u(eta), eta) = v(eta), diff(theta(eta), eta) = r(eta), diff(`ϕ`(eta), eta) = h(eta)};
/ / d \
{ 0.75 f(eta) r(eta) + |----- r(eta)| - 0.05000000000 u(eta)
\ \ deta /

2 / d \ 3
+ 0.2 r(eta) h(eta) + 0.5 r(eta) = 0, |----- v(eta)| + -- f(eta) v(eta)
\ deta / 40

3 2 / d \
- -- u(eta) + theta(eta) - 0.5 ϕ(eta) = 0, |----- h(eta)|
40 \ deta /

+ 3.75 f(eta) h(eta) - 0.1250000000 u(eta) - 1.875000000 f(eta) r(eta)

2
- 0.5000000000 r(eta) h(eta) - 1.250000000 r(eta) = 0,

d d d
----- f(eta) = u(eta), ----- u(eta) = v(eta), ----- theta(eta) = r(eta),
deta deta deta

d \
----- ϕ(eta) = h(eta) }
deta /
> IC := {f(0) = s, h(0) = xi, r(0) = tau, u(0) = 0, v(0) = alpha(0), theta(0) = 1-(1/4)*epsilon1, `ϕ`(0) = (1/4)*epsilon2};
{f(0) = s, h(0) = xi, r(0) = tau, u(0) = 0, v(0) = alpha(0),

theta(0) = 0.9500000000, ϕ(0) = 0.05000000000}
> L := 2;
2
> BC = {u(L) = 0, theta(L) = 0, `ϕ`(L) = 0};
BC = {u(2) = 0, theta(2) = 0, ϕ(2) = 0}
> S := Shoot(ODE, IC, BC, FNS, [alpha = .42453091564332, tau = -.21166705749821127, xi = -.4944583739651814]);
/ / / d \
Shoot|{ 0.75 f(eta) r(eta) + |----- r(eta)| - 0.05000000000 u(eta)
\ \ \ deta /

2 / d \ 3
+ 0.2 r(eta) h(eta) + 0.5 r(eta) = 0, |----- v(eta)| + -- f(eta) v(eta)
\ deta / 40

3 2 / d \
- -- u(eta) + theta(eta) - 0.5 ϕ(eta) = 0, |----- h(eta)|
40 \ deta /

+ 3.75 f(eta) h(eta) - 0.1250000000 u(eta) - 1.875000000 f(eta) r(eta)

2
- 0.5000000000 r(eta) h(eta) - 1.250000000 r(eta) = 0,

d d d
----- f(eta) = u(eta), ----- u(eta) = v(eta), ----- theta(eta) = r(eta),
deta deta deta

d \
----- ϕ(eta) = h(eta) }, {f(0) = s, h(0) = xi, r(0) = tau, u(0) = 0,
deta /

v(0) = alpha(0), theta(0) = 0.9500000000, ϕ(0) = 0.05000000000}, BC,

{f(eta), h(eta), r(eta), u(eta), v(eta), theta(eta), ϕ(eta)}, [

alpha = 0.42453091564332, tau = -0.21166705749821127,

\
xi = -0.4944583739651814]|
/
RungeKutta(ODE, BC, alpha = .42453091564332, tau = -.21166705749821127, xi = -.4944583739651814, output=plot);
/ / / d \
RungeKutta|{ 0.75 f(eta) r(eta) + |----- r(eta)| - 0.05000000000 u(eta)
\ \ \ deta /

2 / d \ 3
+ 0.2 r(eta) h(eta) + 0.5 r(eta) = 0, |----- v(eta)| + -- f(eta) v(eta)
\ deta / 40

3 2 / d \
- -- u(eta) + theta(eta) - 0.5 ϕ(eta) = 0, |----- h(eta)|
40 \ deta /

+ 3.75 f(eta) h(eta) - 0.1250000000 u(eta) - 1.875000000 f(eta) r(eta)

2
- 0.5000000000 r(eta) h(eta) - 1.250000000 r(eta) = 0,

d d d
----- f(eta) = u(eta), ----- u(eta) = v(eta), ----- theta(eta) = r(eta),
deta deta deta

d \
----- ϕ(eta) = h(eta) }, BC, alpha = 0.42453091564332,
deta /

\
tau = -0.21166705749821127, xi = -0.4944583739651814, output = plot|
/
>

Dear sir 

in the above problem im geiitng the problem with , with(shoot) command and even it is not executing at

S := Shoot(ODE, IC, BC, FNS, [alpha = .42453091564332, tau = -.21166705749821127, xi = -.4944583739651814]) this command, here alpha,tau and zi variable should change.