####################################################################
#  Compute curve neighborhoods of a point in B2/B.                 #
#  This reproduces Example 5.3 in the paper                        #
#  "Curve neighborhoods of Schubert varieties"                     #
#  by A. Buch and L. Mihalcea.                                     #
####################################################################


read equivcalc:
with(equivcalc):


printf("B2/B:\n");
OF(5):
for al in pos_roots() do
  zalv := curve_nbhd(weyl_id(), root_degree(al)):
  printf("al = %a ; z_{al^v} = %a\n", al, weyl_redexp(zalv)):
od:


###
