
(* h_bous.m *)

(* Last Updated: August 16, 2012 at 18:00 in Golden *)

debug = True;

(*** Boussinesq System ***)

eq[1][x,t] = D[u[1][x,t],t] + D[u[2][x,t],x];

eq[2][x,t] = D[u[2][x,t],t]+ bb*D[u[1][x,t],x]-3*u[1][x,t]*D[u[1][x,t],x]-
             aa*D[u[1][x,t],{x,3}];

(* 
eq[2][x,t] = D[u[2][x,t],t]+ D[bb*u[1][x,t]-3/2*u[1][x,t]^2-
             aa*D[u[1][x,t],{x,2}],x];
*)

noeqs = 2;
name = "Boussinesq System";
parameters = {aa};
weightpars = {bb};

formrho[x,t] = {};

(* h_bous.m *)

