(* d_kdv.m *) u[1][n_]'[t]:= u[1][n][t]*(u[1][n+1][t]-u[1][n-1][t]); noeqs = 1; name = "Kac-van Moerbeke Equation"; parameters = {}; weightpars = {}; formrho = 0; (* d_kdv.m *)