typo?

Axel Vogt — Sat, 21 Jan 2012 20:57:28 Z

Your last line means, you want to simplify subject to 1/m^2 = 0,
which probably is not what you want.

Impossible equality

Kitonum — Sat, 21 Jan 2012 21:20:49 Z

The equality (1/m)^(N+1) = 0 is impossible for any m!

typo

samantasubho — Sun, 22 Jan 2012 00:30:32 Z

Impossible equality

samantasubho — Sun, 22 Jan 2012 00:29:00 Z

I understand the error. I have corrected it by defining 1/m=epsilon. Now everything works fine, but instead of showing all the three equations after series substitution MAPLE is showing only Eq [3]. But it should show all equations together. Can you please check if I am logically incorrect somewhere?

restart; N := 2

alias(eta = e, theta = t)

Eq[1] := 5*(diff(F(e), `$`(e, 3)))+(1+3/m)*F(e)*(diff(F(e), `$`(e, 2)))-(2+1/m)*(diff(F(e), e))^2-(4+2/m)*H(e)-(1-2/m)*e*(diff(H(e), e)) = 0

Eq[2] := diff(H(e), e) = t(e)

Eq[3] := 5*(diff(t(e), `$`(e, 2)))/Pr+(1+3/m)*F(e)*(diff(t(e), e))-5*(diff(F(e), e))*t(e)

pars := 1/m = epsilon

for i to 3 do eq[i] := subs(pars, Eq[i]) end do

eqs := eq[1], eq[2], eq[3]

F1 := unapply(add(epsilon^i*F[i](e), i = 0 .. N), e)

F2 := unapply(add(epsilon^i*H[i](e), i = 0 .. N), e)

F3 := unapply(add(epsilon^i*t[i](e), i = 0 .. N), e)

deqn := simplify(collect(subs({F(e) = F1(e), H(e) = F2(e), t(e) = F3(e)}, eqs), epsilon), {epsilon^(N+1) = 0})

Thanks in advance.Perturbation.mw

still

Axel Vogt — Sun, 22 Jan 2012 01:33:49 Z

1/m = epsilon does not matter, You still say 1/m = epsilon = 0.

MaplePrimes - answers and comments on Question, Error using simplify

typo?

Impossible equality

typo

Impossible equality

still