(* d_stdnls.m *) u[1][n_]'[t]:= aa*(u[1][n+1][t]-2*u[1][n][t]+u[1][n-1][t]) + 2*u[1][n][t]^2*u[2][n][t]; u[2][n_]'[t]:= -aa*(u[2][n+1][t]-2*u[2][n][t]+u[2][n-1][t]) - 2*u[1][n][t]*u[2][n][t]^2; noeqs = 2; name = "NLS Equation (Standard Discretization)"; parameters = {}; weightpars = {aa}; weightu[1] = weightu[2]; formrho = 0; (* d_stdnls.m *)