(* d_ablnls.m *)

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

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

noeqs = 2;
name = "NLS Equation (Ablowitz-Ladik Discretization)";
parameters = {};
weightpars = {aa};

weightu[1] = weightu[2];

formrho = 0;

(* d_ablnls.m *)