Three (3) constraint equations were added 
G0001.org  -----> GOO01.mod 
D0001.org  -----> DOO01.mod 
E0001.org  -----> EOO01.mod 

The constraints are at the bottom of the G and D matrices

Run delaytmR.m in Octave using the *.mod files for input.
