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