/* Various equations in bilinear form */
/* This compilation of command files produces all the output files */

writefile("h_kdv.out");
batchload("hir_sing.max")$  
n:3$
B(f,g):= Dxt[1,1](f,g) + Dx[4](f,g)$ 
name:KdV$
hirota(b,name,n,2,3,true,true)$
kill(all)$
closefile();

writefile("h_sk.out");
batchload("hir_sing.max")$  
n:3$
B(f,g):= Dxt[1,1](f,g) + Dx[6](f,g)$ 
name:Sawada_Kotera$
hirota(b,name,n,3,2,true,true)$
kill(all)$
closefile();

writefile("h_7kdv.out");
batchload("hir_sing.max")$  
n:3$
B(f,g):= Dxt[1,1](f,g) + Dx[8](f,g)$ 
name:Seventh_KdV$
hirota(b,name,n,2,2,true,true)$
kill(all)$
closefile();

writefile("h_m5kdv.out");
batchload("hir_sing.max")$  
n:3$
B(f,g):= Dxt[3,1](f,g) + Dx[6](f,g)$ 
name:modified_5th_KdV$
hirota(b,name,n,1,3,true,true)$
kill(all)$
closefile();

writefile("h_itob4.out");
batchload("hir_sing.max")$  
n:3$
B(f,g):= Dxt[1,1](f,g) + Dxt[2,2](f,g)$ 
name:Ito_b4$
hirota(b,name,n,3,1,true,true)$
kill(all)$
closefile();

writefile("h_itob6.out");
batchload("hir_sing.max")$  
n:3$
B(f,g):= Dxt[0,2](f,g) + Dxt[3,1](f,g) +Dy[2](f,g)$ 
name:Ito_b6$
hirota(b,name,n,2,2,true,true)$
kill(all)$
closefile();

writefile("h_itob8.out");
batchload("hir_sing.max")$  
n:3$
B(f,g):= - Dxt[2,2](f,g) - Dxt[0,2](f,g) + Dx[2](f,g)$ 
name:Ito_b8$
hirota(b,name,n,1,1,true,true)$
kill(all)$
closefile();

writefile("h_itob9.out");
batchload("hir_sing.max")$  
n:3$
B(f,g):= Dxt[0,2](f,g) + Dx[6](f,g)$ 
name:Ito_b9$
hirota(b,name,n,1,1,true,true)$
kill(all)$
closefile();

writefile("h_kp.out");
batchload("hir_sing.max")$  
n:3$
B(f,g):= Dxt[1,1](f,g) + Dx[4](f,g) + 3*Dy[2](f,g)$ 
name:Kadomtsev_Petviashvili$
hirota(b,name,n,1,1,true,false)$
kill(all)$
closefile();

writefile("h_bous.out");
batchload("hir_sing.max")$  
n:3$
B(f,g):= Dxt[0,2](f,g) - Dx[2](f,g) - Dx[4](f,g)$ 
name:Boussinesq$
hirota(b,name,n,1,1,true,false)$
kill(all)$
closefile();

writefile("h_shwat.out");
batchload("hir_sing.max")$  
n:3$
B(f,g):= Dxt[3,1](f,g) - Dx[2](f,g) - Dxt[1,1](f,g)$ 
name:Shallow_Water$
hirota(b,name,n,1,1,true,false)$
kill(all)$
closefile();

writefile("h_param.out")$
batchload("hir_sing.max")$ 
N:3$
/* kpar:-1/5; */
B(f,g):=kpar*Dt[2](f,g)+Dxt[3,1](f,g)+Dx[6](f,g)$ 
name:Parameter_Dependent$
hirota(B,name,N,2,2,true,true)$
kill(all)$
closefile();

quit();

/* end of command file h_vareqs.com for various data */