reflect(R,i,p)
i1 : R = rootSystemB(4) o1 = RootSystem{...8...} o1 : RootSystem
i2 : p = halfSumOfRoots(R) o2 = | 1 | | 1 | | 1 | | 1 | 4 o2 : ZZ
i3 : reflect(R,4,p) o3 = | 1 | | 1 | | 2 | | -1 | 4 o3 : ZZ