Compute divided differencing matrix of order O
Smooths the Y vs.
Smooths Y vs.
Wrapper function for rgdtsmcore in order to minimize over LAMBDA w.r.t.