R:= RootSystem( "B",3);
 U:=QuantizedUEA(R); 
 g:= GeneratorsOfAlgebra( U ); 
 PositiveRoots(R); 
 PositiveRootsInConvexOrder(R); 
########## generators of U ###############
 f1:=g[1]; f2:=g[3]; f3:=g[9];
 e1:=g[16]; e2:=g[18]; e3:=g[24];
 k1:=g[10]; k1m:=g[11]; k2:=g[12]; k2m:=g[13]; k3:=g[14]; k3m:=g[15];
########## generators of Uq'(k) ###############
 b1:=f1-k1m*e1;
 b2:=f2-k2m*e2;
 b3:=f3-k3m*e3;
########## q-numbers #################### 
 2q:=_q+_q^-1;
 3q:=_q^2+1+_q^-2;
######### relations in Uq'(k) as functions #########
 rel0:=function(a,b)
   return a*b - b*a;
 end;
 rel1:=function(a,b)
   return a^2*b - (_q^2+_q^-2)*a*b*a + b*a^2 +_q^-2 *b;
 end;
 rel2:=function(a,b)
   return b^3*a - 3q *b^2*a*b + 3q*  b*a*b^2 - a*b^3 +_q^-1*2q^2 *(b*a-a*b);
 end;
###############################################################################
########### We check that the generators bi really satisfy the relations ######
########### given in Proposition 3.1.                                    ######
###############################################################################
 rel0(b1,b3);
 rel1(b1,b2);
 rel1(b2,b1);
 rel1(b2,b3);
 rel2(b2,b3);
######## Definition of \tau_1^- ##########
 tau1mb1 := b1;;
 tau1mb2 := b1*b2 - _q^2 *b2*b1;;
 tau1mb3 := b3;;
####### \tau_1^- is an algebra homomorphism #####
 rel0(tau1mb1,tau1mb3);
 rel1(tau1mb1,tau1mb2);
 rel1(tau1mb2,tau1mb1);
 rel1(tau1mb2,tau1mb3);
 rel2(tau1mb2,tau1mb3);
######## Definition of \tau_1 ##########
 tau1b1 := b1;;
 tau1b2 := b2*b1 - _q^2 *b1*b2;;
 tau1b3 := b3;;
####### \tau_1 is an algebra homomorphism #####
 rel0(tau1b1,tau1b3);
 rel1(tau1b1,tau1b2);
 rel1(tau1b2,tau1b1);
 rel1(tau1b2,tau1b3);
 rel2(tau1b2,tau1b3);
###### \tau_1 is the inverse of \tau_1^- #####
  tau1b1*tau1b2 - _q^2 *tau1b2*tau1b1 - b2;
  tau1mb2*tau1mb1 - _q^2 *tau1mb1*tau1mb2 -b2;
######## Definition of \tau_2^- ##########
 tau2mb1 := b2*b1-_q^2 * b1*b2;;
 tau2mb2 := b2;;
 tau2mb3 := b2*b3-_q^2 * b3*b2;;
####### \tau_2^- is an algebra homomorphism #####
 rel0(tau2mb1,tau2mb3);
 rel1(tau2mb1,tau2mb2);
 rel1(tau2mb2,tau2mb1);
 rel1(tau2mb2,tau2mb3);
 rel2(tau2mb2,tau2mb3);
######## Definition of \tau_2 ##########
 tau2b1 := b1*b2-_q^2 * b2*b1;;
 tau2b2 := b2;;
 tau2b3 := b3*b2-_q^2 * b2*b3;;
####### \tau_2 is an algebra homomorphism #####
 rel0(tau2b1,tau2b3);
 rel1(tau2b1,tau2b2);
 rel1(tau2b2,tau2b1);
 rel1(tau2b2,tau2b3);
 rel2(tau2b2,tau2b3);
####### \tau_2 is the inverse of \tau_2^- #########
 tau2b2*tau2b1-_q^2 * tau2b1*tau2b2 -b1;
 tau2b2*tau2b3-_q^2 * tau2b3*tau2b2-b3;
 tau2mb1*tau2mb2-_q^2 * tau2mb2*tau2mb1 -b1;
 tau2mb3*tau2mb2-_q^2 * tau2mb2*tau2mb3-b3;
######## Definition of \tau_3^- ########################## 
 tau3mb1 := b1;;
 tau3mb2 := 2q^-1*b3^2*b2-_q*b3*b2*b3+_q^2*2q^-1*b2*b3^2+  b2;;
 tau3mb3 := b3;;
####### \tau_3^- is an algebra homomorphism #####
 rel0(tau3mb1,tau3mb3);
 rel1(tau3mb1,tau3mb2);
 rel1(tau3mb2,tau3mb1);
 rel1(tau3mb2,tau3mb3);
 rel2(tau3mb2,tau3mb3);
######## Definition of \tau_3 ########################## 
 tau3b1 := b1;;
 tau3b2 := 2q^-1*b2*b3^2 -_q*b3*b2*b3+_q^2*2q^-1*b3^2*b2 +  b2;;
 tau3b3 := b3;;
####### \tau_3^- is an algebra homomorphism #####
 rel0(tau3mb1,tau3mb3);
 rel1(tau3mb1,tau3mb2);
 rel1(tau3mb2,tau3mb1);
 rel1(tau3mb2,tau3mb3);
 rel2(tau3mb2,tau3mb3);
####### \tau_3^- is the inverse of \tau_3 ##########
  2q^-1*tau3b3^2*tau3b2-_q*tau3b3*tau3b2*tau3b3+_q^2*2q^-1*tau3b2*tau3b3^2+  tau3b2 - b2;
  2q^-1*tau3mb2*tau3mb3^2-_q*tau3mb3*tau3mb2*tau3mb3+_q^2*2q^-1*tau3mb3^2*tau3mb2 +tau3mb2-b2;
##################################################################
######## Verifying the braid relations ###########################
######## We write t2323 for tau_2 tau_3 tau_2(b3) etc. ###########
##################################################################
t11:= b1;;
t12:= b2*b1 - _q^2 *b1*b2;;
t13:= b3;;
t21:= b1*b2-_q^2 * b2*b1;;
t22:= b2;;
t23:= b3*b2-_q^2 * b2*b3;;
t31:= b1;;
t32:= 2q^-1*b2*b3^2 -_q*b3*b2*b3+_q^2*2q^-1*b3^2*b2 +  b2;;
t33:= b3;;
############
t311:=t31;;
t312:=t32*t31 - _q^2 *t31*t32;;
t313:=t33;;
t131:=t11;;
t132:=2q^-1*t12*t13^2 -_q*t13*t12*t13+_q^2*2q^-1*t13^2*t12 +  t12;;
t133:=t13;;
############
t121:= t11*t12-_q^2 * t12*t11;;
t122:= t12;;
t123:= t13*t12-_q^2 * t12*t13;;
##
t211:= t21;;
t212:= t22*t21 - _q^2 *t21*t22;;
t213:= t23;;
##
t231:=t21;;
t232:=2q^-1*t22*t23^2 -_q*t23*t22*t23+_q^2*2q^-1*t23^2*t22 +  t22;;
t233:=t23;;
##
t321:=t31*t32-_q^2 * t32*t31;;
t322:=t32;;
t323:=t33*t32-_q^2 * t32*t33;;
############
t2121:= t211*t212-_q^2 * t212*t211;;
t2122:= t212;;
t2123:= t213*t212-_q^2 * t212*t213;;
##
t1211:= t121;;
t1212:= t122*t121 - _q^2 *t121*t122;;
t1213:= t123;;
##
t3231:=t321;;
t3232:=2q^-1*t322*t323^2 -_q*t323*t322*t323+_q^2*2q^-1*t323^2*t322 +  t322;;
t3233:=t323;;
##
t2321:=t231*t232-_q^2 * t232*t231;;
t2322:=t232;;
t2323:=t233*t232-_q^2 * t232*t233;;
############
t23231:=t2321;;
t23232:=2q^-1*t2322*t2323^2 -_q*t2323*t2322*t2323+_q^2*2q^-1*t2323^2*t2322 +  t2322;;
t23233:=t2323;;
##
t32321:=t3231*t3232-_q^2 * t3232*t3231;;
t32322:=t3232;;
t32323:=t3233*t3232-_q^2 * t3232*t3233;;
######## braid relations #################################
t132-t312;
t2121-t1211;
t2122-t1212;
t2123-t1213;
t32321-t23231;
t32322-t23232;
t32323-t23233;
######## :-) End :-) #####################################

#############################################################
############ Relation to Lusztig automorphisms       ########
############ Here we check Relation (3.2) for aij=-2 ########
#############################################################
tau:=AntiAutomorphismTau( U );;
T2:=AutomorphismTalpha(U,2);;
T3:=AutomorphismTalpha(U,3);;
T2m:= tau*T2*tau;;
T3m:= tau*T3*tau;;
Image(T2m,b2)-tau2mb2;;
Image(T2m,b3)-tau2mb3;;
Image(T3m,b2)-tau3mb2;;
Em2:=-(_q-_q^-1)*(_q^-1*f2*k3m^2 +_q*(_q-_q^-1)*f2*k3m^2*e3^2+(f3*f2-_q^2* f2*f3)*k3m*e3);;
Image(T3m,b2)-tau3mb2-Em2;
#############################################################