Find Hirota Bilinear form using Maple

Ask a Question

Create a Post

Question:Find Hirota Bilinear form using Maple

Posted: alaloush 35 Product: Maple 2018

February 25 2020

I have read the article attached in the following

https://drive.google.com/open?id=1BLzklSFx3nPovLC_2L1muAWY1kFnbwLV

I wrote the code mentioned in the article in order to find the Hirota formula This is the code that I wrote For Hirota Method But when running the program it does not work well

Could you help me to fix the mistake and run the program,plz ?

(1)

BD := proc (FF, DD) local f, g, x, m, opt; if nargs = 1 then return `*`(FF[]) end if; f, g := FF[]; x, m := DD[]; opt := args[3 .. -1]; if m = 0 then return procname(FF, opt) end if; procname([diff(f, x), g], [x, m-1], opt)-procname([f, diff(g, x)], [x, m-1], opt) end proc

(2)

(3)

(4)

print_BD := proc (FF, DD) local f, g, x, m, i; f, g := FF[]; f := cat(f, '', g); g := product(D[args[i][1]]^args[i][2], i = 2 .. nargs); if g <> 1 then f := g*(`*`(f)) end if; f end proc

printHBT := proc (uf, u, f, i, j, PL, alpha := 1) local DD, g, C, tmp, pl; pl := printHBF(PL); if 0 < j then print(u = 2*alpha*('diff')(ln(f), `$`(x, j))) else print(u = 2*alpha*ln(f)) end if; if 0 < i then print(('diff')(pl/f^2, `$`(x, i)) = 0) else print(pl/f^2 = 0) end if; NULL end proc

$guessalpha := proc (Res, uf, u, f, i, j, PL, alpha) local tmp, res, pl, flag, k; flag := 1; tmp := [op(Res)]; tmp := map(numer, tmp); tmp := gcd(tmp[1], tmp[-1]); if type(tmp, `*`) then tmp := remove(has, tmp, f) end if; if tmp <> 0 and has(tmp, {alpha}) then tmp := solve(tmp/alpha^difforder(uf), {alpha}); if tmp <> NULL and has(tmp, {alpha}) then lprint(tmp); for k to nops([tmp]) while flag = 1 do res := collect(expand(subs(tmp[k], Res)), f, factor); if res = 0 then pl := subs(tmp[k], PL); printHBT(uf, u, f, i, j, pl, rhs(tmp[k])); flag := 0 end if end do end if end if; PL end proc$

eq := diff(u(x, t), t, t)-(diff(u(x, t), x, x))-3*(diff(u(x, t), x, x))^2-(diff(u(x, t), x, x, x, x)) = 0

(5)

diff(diff(u(x, t), t), t)-(diff(diff(u(x, t), x), x))-3*(diff(diff(u(x, t), x), x))^2-(diff(diff(diff(diff(u(x, t), x), x), x), x)) = 0

Error, (in getFnumer) expected 1 st argument be type `+`.

Download Hirota_to_send.mw

227 views
More...

Contact Author

Subscribe

Generate PDF

Answer

Reply

Share

Branch

Create a new Post in reply

Ask a related Question

Flag

Answer this Question Add a Comment

Share via:

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

Ask a Question

Create a Post

Question: Find Hirota Bilinear form using Maple

Question:Find Hirota Bilinear form using Maple

227 views
More...

Contact Author

Subscribe

Generate PDF

Answer

Reply

Share

Branch

Create a new Post in reply

Ask a related Question

Flag

Preview

Preview

Save this setting as your default sorting preference?

Ask a Question

Create a Post

Generating PDF…

Save this setting as your default sorting preference?
Note: You can change your preference any time in your account settings
Don't show this again

Ask a Question

Create a Post

Question: Find Hirota Bilinear form using Maple

Question:Find Hirota Bilinear form using Maple

227 views More... Contact Author Subscribe Generate PDF Answer Reply Share Branch Create a new Post in reply Ask a related Question Flag

Preview

Preview

Save this setting as your default sorting preference?

Ask a Question

Create a Post

Generating PDF…

227 views
More...

Contact Author

Subscribe

Generate PDF

Answer

Reply

Share

Branch

Create a new Post in reply

Ask a related Question

Flag