reset;
model zermelo.mod;
data  zermelo.dat;
option solver ipopt;
#option solver snopt;
solve;
print: {s in 1..5} m > zermelo.out; 
print {s in 0..m}: t[s],u[s,1],x[s,1],x[s,2],x[s,3] > zermelo.out;

