/* -*- mode: C -*- */ /* $Id$ */ load("yang_D.rr")$ C=c$ V=[x,y,z]$ yang.define_ring(V)$ S=dx+dy+dz$ L1 = yang.mul(dx,S+C-1) - x*yang.mul(dx+a1,dx+b1)$ L2 = yang.mul(dy,S+C-1) - y*yang.mul(dy+a2,dy+b2)$ L3 = yang.mul(dz,S+C-1) - z*yang.mul(dz+a3,dz+b3)$ L = [L1,L2,L3]$ Gr=yang.gr(L)$ sprintf("Rank=%a",length(yang.stdmon(Gr))); Base = [1,dx,dy,dz,dx*dy,dx*dz,dy*dz,dx*dy*dz]$ Pf = yang.pf(Base,Gr)$ vtol(map(matrix_matrix_to_list,Pf)); /* print(Pf[0])$ print("")$ print(Pf[1])$ print("")$ print(Pf[2])$ */ /* Check = yang.check_pfaffian(Pf,V); */ end$