(* d_flie3.m *) debug = False; (* Example 3 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 = "Fairlie3 System"; parameters = {}; weightpars = {}; (* Possible weight that could be given in this data file *) (* weight[t] = 4; *) formrho[x,t] = {}; (* d_flie3.m *)