> 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.