@janhardo you answer all question you are going to modify your Ai i like that , just last my question he can't answer thank you it work i have to check it then

@dharr let me give you a hint maybe already you know this, but lineaer part of pde is the same of bilinear equation, and 5 month before macdara write a program for me about how open this derivative and it is totally corect for (1+1) dimension when it is (x,t) but for (x,y,t) and (x,y,z,t) i don't know why not working , i want to check that if i open bilinear equation is the same as linear part or just the shape of bilinear form it is , i will explain more if you needed , 

i though equation 25 is linear part of my pde, the bilinear form it is i am sure %100 percent, but when i open bilinear equation it is equal to linear part or not i am not sure, i have to open by hirota derivstive in here i try it but if watch and make it for 3+1 dimension the prcecure we will be sure about that it is bilinear form or not 

also i check the result i don't know something must be mistake otherwise your code are working so good 

Download test-Blinear.mw

Download system-pdetest.mw

@janhardo i have a lot problem not just this, so i need a systematic progress

@janhardo replaced thus parameter in L2 then L2 in H then do pdetest if is zero then it is correct.

@ecterrab maybe  question is not related to dirac but i ask more than 10 time about this question and i didn't get any answer i want apply long wave limit which i am quit sure it is a technique, i asked a lot time to a lot people but i didn't get my answer i hope you have  idea about  this question,

https://www.mapleprimes.com/questions/240218-How-Apply-Substitution-And-Apply-Limit

i know every structure if have any question please let me know!.

@ecterrab  i have some question in a topic maybe you faced before can i ask the question ?

@janhardo  i know every structure of the topic in theory but when i come to apply by maple is making problem in 2 part finding long wave limit of series and finding paicewise function  of Bij, which have a trick i don't know how reach is easy  but how reach that i don't know yet

@janhardo  i will ask in here each week maybe someone can figure out

@janhardo  wrong -this two theta ^2 are extra

@janhardo  this is wrong the function don't have exponential even! 

@janhardo your ai define the function inside the proc, you have to use N-soliton series then have to idea first do change of variables then apply limit epsilon-->0 or k[1]-k[2]-->0

@janhardo  all wrong except F[3] which already define in procedure, thanks john take some rest

@janhardo i watched is totaly wrong becuase you write the function we already define and just add so it is wrong we have to apply something to N-soliton series to get F[1],F[2],F[3],F[4],F[5] in here we didn't apply long wave limit if you watch the procedure it use fucntion which already write down so don't have mean like that we don't reach out point i will search more and try to apply limit and tylor series

@janhardo  for N=4 we have this function 

@janhardo  this function is true i want other be like that can we do that?

LumpSolitonTauFunction(3);
 ((B23 + theta2) theta1 + B13 B23 + B13 theta2 + B12) exp(eta3)

    + theta1 theta2 + B12