(* d_toda.m *)

(* Toda Lattice *)

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

noeqs = 2;
name = "Toda Lattice";
parameters = {};
weightpars = {};

formrho = 0;

(* d_toda.m *)