
Needs["hirota`"]

debugwilly=False;

n=3;

B[f_, g_]:=Dxt[0,2,x,t,f,g]-Dx[2,x,f,g]-Dx[4,x,f,g];

name="Boussinesq";

myfile="h_bous.o";

OpenLog[myfile];

(* Do not run this case with the symbolic test for 4 solitons set to true *)
(* unless you have a lot of patience or a large computer *)

hirota[B,name,n,2,1,True,False];

CloseLog[];

