-- -*- M2-comint -*- {* hash: 856279973 *}

i1 : R = QQ[x_1..x_4,D_1..D_4, WeylAlgebra=>(apply(4,i->x_(i+1)=>D_(i+1)))] 

o1 = R

o1 : PolynomialRing

i2 : stafford ideal (D_1,D_2,D_3,D_4)

                 4        2        3
o2 = ideal (D , x x D  + x x D  + x D  + x D  + D )
             1   1 4 4    1 3 3    1 4    1 3    2

o2 : Ideal of QQ[x , x , x , x , D , D , D , D ]
                  1   2   3   4   1   2   3   4

i3 :