
(* d_flie2.m *)

debug = False;

(* Example 2 from preprint by David Fairlie *)

eq[1][x,t] = D[u[1][x,t],t]-
 (u[4][x,t]*u[2][x,t]+u[2][x,t]*u[3][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[1][x,t]*u[3][x,t]+u[3][x,t]*u[4][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[2][x,t]*u[4][x,t]+u[4][x,t]*u[1][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[3][x,t]*u[1][x,t]+u[1][x,t]*u[2][x,t]+u[2][x,t]*u[3][x,t])*D[u[4][x,t],x];

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

formrho[x,t] = {};

(* d_flie2.m *)
