Needs["hirota`"]

n=3;

B[f_, g_]:=Dx[4,x,f,g]+Dxt[1,1,x,t,f,g]+3*Dx[2,y,f,g];

name="Kadomtsev-Petviashvili";

myfile="h_kp.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,2,True,False];

CloseLog[];