/* -*- mode: C -*- */ /* $Id$ */ /* Horn's G1 */ load("yang_D.rr")$ V=[x,y]$ yang.define_ring(V)$ P=dx/x$ Q=dy/y$ R=(dx^2-dx)/x^2$ S=dx*dy/(x*y)$ T=(dy^2-dy)/y^2$ U1 = x*(1+x)*R - y*S -y^2*T +(1-b1+(a+b2+1)*x)*P +(b2-a-1)*y*Q + a*b2$ U2 = y*(1+y)*T - x*S -x^2*R +(1-b2+(a+b1+1)*y)*Q +(b1-a-1)*x*P + a*b1$ /* L1 = base_cancel(x*U1)$ */ /* L2 = base_cancel(y*U2)$ */ L1 = yang.mul(dx,dx-dy-b1) + x*yang.mul(dx+dy+a,dx-dy+b2)$ L2 = yang.mul(dy,dy-dx-b2) + y*yang.mul(dx+dy+a,dy-dx+b1)$ L = [L1,L2]$ Gr=yang.gr(L); Base=[1,dx,dy]$ Pf = yang.pf(Base,Gr)$ vtol(map(matrix_matrix_to_list,Pf)); /* print(Pf[0])$ print("")$ print(Pf[1])$ */ end$