(* d_flie1.m *)

debug = False;

(* Example 1 from preprint by David Fairlie *)

eq[1][x,t] = D[u[1][x,t],t]-(u[2][x,t]+u[3][x,t]+u[4][x,t])*D[u[1][x,t],x];
eq[2][x,t] = D[u[2][x,t],t]-(u[3][x,t]+u[4][x,t]+u[1][x,t])*D[u[2][x,t],x];
eq[3][x,t] = D[u[3][x,t],t]-(u[4][x,t]+u[1][x,t]+u[2][x,t])*D[u[3][x,t],x];
eq[4][x,t] = D[u[4][x,t],t]-(u[1][x,t]+u[2][x,t]+u[3][x,t])*D[u[4][x,t],x];

noeqs = 4;
name = "Fairlie1 System";
parameters = {};
weightpars = {};

formrho[x,t] = {};

(* d_flie1.m *)