/* -*- mode: C -*- */ /* $Id$ */ /* Horn's confluent Eta2 */ load("yang_D.rr")$ V=[x,y]$ yang.define_ring(V)$ L = [ yang.mul(dx,dx+d-1)-x*yang.mul(dx-dy+a,dx+b), yang.mul(dy,dy-dx-a)+y*(dy+c) ]; Gr=yang.gr(L); /* yang.stdmon(Gr);end; */ sprintf("Rank=%a",length(yang.stdmon(Gr))); Base=[1,dx,dy,dx*dy]$ Pf = yang.pf(Base,Gr)$ vtol(map(matrix_matrix_to_list,Pf)); /* print(Pf[0])$ print("")$ print(Pf[1])$ */ end$