mathematicianS

5 Reputation

One Badge

9 years, 40 days

MaplePrimes Activity


These are replies submitted by mathematicianS

@Preben Alsholm 

This is the coding for Maple... How to write the IC for u(0)? That boundary condition contains of f'''(0)... Thanks for your help! =)


restart

Shootlib := "C:\Users\Asus\\maple/"

"C:UsersAsus\maple/"

(1)

libname := Shootlib, libname

with(Shoot)

with(plots)

 

`ς` := 10

10

(2)

a := 1

1

(3)

b := -1

-1

(4)

ss := 2

2

(5)

Pr := 1

1

(6)

lambda := 0

0

(7)

sigma := -2

-2

(8)

FNS := {f(eta), s(eta), t(eta), u(eta), v(eta)}

ODE := {(diff(t(eta), eta))/Pr+f(eta)*t(eta)-u(eta)*s(eta) = 0, diff(v(eta), eta)+f(eta)*v(eta)-u(eta)^2+lambda*s(eta) = 0, diff(f(eta), eta) = u(eta), diff(s(eta), eta) = t(eta), diff(u(eta), eta) = v(eta)}

{diff(t(eta), eta)+f(eta)*t(eta)-u(eta)*s(eta) = 0, diff(v(eta), eta)+f(eta)*v(eta)-u(eta)^2 = 0, diff(f(eta), eta) = u(eta), diff(s(eta), eta) = t(eta), diff(u(eta), eta) = v(eta)}

(9)

IC := {f(0) = ss, s(0) = 1, t(0) = beta, u(0) = sigma+a*alpha+b*(diff(v(0), 0)), v(0) = alpha}

Error, invalid input: diff received 0, which is not valid for its 2nd argument

 

NULL

BC := {s(`ς`) = 0, u(`ς`) = 0}

{s(10) = 0, u(10) = 0}

(10)

infolevel[shoot] := 1

S := shoot(ODE, IC, BC, FNS, [alpha = .35, beta = 20])

Error, (in Shoot:-shoot) invalid input: process_required_args expects its 2nd argument, IC, to be of type {list, set, equation}, but received IC

 

``


Download 2013.mw2013.mw

Page 1 of 1