sharena2

50 Reputation

5 Badges

11 years, 86 days

MaplePrimes Activity


These are replies submitted by sharena2

@J4James ok, thank you.. im using maple12.. huhu..

@J4James may i know which version maple you used. becoz it not working with my version..

@J4James thank you..  I will try to do as suggested

@mehdi jafari thanks mehdi.. i try three times to upload the file but it failed.. but its ok.. actually i just need change dots (.) to multiplication signs (*) all over.

@Preben Alsholm thank you so much.. its working.. ;)

@Preben Alsholm OWh i will try as your suggestion.. thanks :)

@mehdi jafari 

restart

with*plots

with*plots

(1)

Eq1 := diff(f(eta), eta, eta, eta)+Typesetting:-delayDotProduct(f(eta), diff(f(eta), eta, eta))-(diff(f(eta), eta))^2-Typesetting:-delayDotProduct(M, diff(f(eta), eta))+Typesetting:-delayDotProduct(Typesetting:-delayDotProduct(B, H(eta)), F(eta)-(diff(f(eta), eta))) = 0

diff(diff(diff(f(eta), eta), eta), eta)+f(eta).(diff(diff(f(eta), eta), eta))-(diff(f(eta), eta))^2-M.(diff(f(eta), eta))+`.`(B, H(eta), F(eta)-(diff(f(eta), eta))) = 0

(2)

Eq2 := Typesetting:-delayDotProduct(G(eta), diff(F(eta), eta))+F(eta)^2+Typesetting:-delayDotProduct(B, F(eta)-(diff(f(eta), eta))) = 0

G(eta).(diff(F(eta), eta))+F(eta)^2+B.(F(eta)-(diff(f(eta), eta))) = 0

(3)

Eq3 := Typesetting:-delayDotProduct(G(eta), diff(G(eta), eta))+Typesetting:-delayDotProduct(B, f(eta)+G(eta)) = 0

G(eta).(diff(G(eta), eta))+B.(f(eta)+G(eta)) = 0

(4)

Eq4 := Typesetting:-delayDotProduct(G(eta), diff(H(eta), eta))+Typesetting:-delayDotProduct(H(eta), diff(G(eta), eta))+Typesetting:-delayDotProduct(F(eta), H(eta)) = 0

G(eta).(diff(H(eta), eta))+H(eta).(diff(G(eta), eta))+F(eta).H(eta) = 0

(5)

bcs1 := f(0) = R, (D(f))(0) = 1, (D(f))(10) = 0;

f(0) = R, (D(f))(0) = 1, (D(f))(10) = 0

(6)

fixedparameter := [M = 3, B = 0.2e-1, R = 0];

[M = 3, B = 0.2e-1, R = 0]

(7)

Eq5 := eval(Eq1, fixedparameter);

diff(diff(diff(f(eta), eta), eta), eta)+f(eta).(diff(diff(f(eta), eta), eta))-(diff(f(eta), eta))^2-3*(diff(f(eta), eta))+0.2e-1*(H(eta).(F(eta)-(diff(f(eta), eta)))) = 0

(8)

Eq6 := eval(Eq2, fixedparameter);

G(eta).(diff(F(eta), eta))+F(eta)^2+0.2e-1*F(eta)-0.2e-1*(diff(f(eta), eta)) = 0

(9)

Eq7 := eval(Eq3, fixedparameter);

G(eta).(diff(G(eta), eta))+0.2e-1*f(eta)+0.2e-1*G(eta) = 0

(10)

bcs2 := eval(bcs1, fixedparameter);

f(0) = 0, (D(f))(0) = 1, (D(f))(10) = 0

(11)

bcs3 := F(10) = 0;

F(10) = 0

(12)

bcs4 := G(10) = -f(10);

G(10) = -f(10)

(13)

bcs5 := H(10) = n;

H(10) = n

(14)

L := [.2];

[.2]

(15)

for k to 1 do R := dsolve(eval({Eq4, Eq5, Eq6, Eq7, bcs2, bcs3, bcs4, bcs5}, n = L[k]), [f(eta), F(eta), G(eta), H(eta)], numeric, output = listprocedure) end do

``


Download sharena.mw

1 2 3 Page 3 of 3