/* -*- mode: C -*- */ /* $Id$ */ /* Horn's H4 */ load("yang_D.rr")$ V=[x,y]$ yang.define_ring(V)$ L1 = yang.mul(dx,dx+c-1)-x*yang.mul(2*dx+dy+a,2*dx+dy+a+1)$ L2 = yang.mul(dy,dy+d-1)-y*yang.mul(2*dx+dy+a,dy+b)$ L = [L1,L2]$ Gr=yang.gr(L); 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$