(* d_dual.m *) u[1][n_]'[t]:= (aa+u[1][n][t]^2)*(u[2][n][t]-u[2][n-1][t]); u[2][n_]'[t]:= (aa+u[2][n][t]^2)*(u[1][n+1][t]-u[1][n][t]); noeqs = 2; name = "Self-Dual Network Equations"; parameters = {}; weightpars = {aa}; formrho = 0; (* d_dual.m *)