(* d_mkdv.m *)

u[1][n_]'[t]:= (aa+u[1][n][t]^2)*(u[1][n+1][t]-u[1][n-1][t]);

noeqs = 1;
name = "Modified KdV Equation";
parameters = {};
weightpars = {aa};

formrho = 0;

(* d_mkdv.m *)