(* d_flie4.m *) debug = False; (* Example 4 from preprint by David Fairlie, system (4.1) *) eq[1][x,t] = D[u[1][x,t],t]- ((n1-1)*u[1][x,t]+(n2)*u[2][x,t]+(n3)*u[3][x,t])*D[u[1][x,t],x]; eq[2][x,t] = D[u[2][x,t],t]- ((n1-1)*u[1][x,t]+(n2-1)*u[2][x,t]+(n3)*u[3][x,t])*D[u[2][x,t],x]; eq[3][x,t] = D[u[3][x,t],t]- ((n1)*u[1][x,t]+(n2)*u[2][x,t]+(n3-1)*u[3][x,t])*D[u[3][x,t],x]; noeqs = 3; name = "Fairlie4 System"; parameters = {n1,n2,n3}; weightpars = {}; formrho[x,t] = {}; (* d_flie4.m *)