function diff_operator(nf, order) creates a [nf-order]x[nf] differential operator matrix D so that D*x = diff(x, order)