kite v.s. sunset

gistfile1.txt

Maxima 5.41.0 http://maxima.sourceforge.net
using Lisp GNU Common Lisp (GCL) GCL 2.6.12
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) batch("kiteVSsunset.mac")

read and interpret file: #p/home/ray/Documents/latex/IBP/Schwinger/kiteVSsunset.mac
(%i2) kill(all)
(%i1) load("Schwinger.mac")
The list of functions:
  schwinger_matrix(list_of_loop_momenta, list_of_denominators)
  
  symanzik_u(ls,ds) 
  symanzik_f(ls,ds) 
  lee_pomeransky_g(ls,ds)
  homogeneous_lee_pomeransky_g(ls,ds)
As helper functions, the following functions are used:
  det_but_outer_egdes(matrix) 
  lee_pomeransky_g_from_matrix(schwinger_matrix)
 
  schwinger_parameters(ds)
  schwinger_exponents(ds) 
 
  at_zero(ls) 
  c_of_schwinger_matrix(ls,ds) 
  b_of_schwinger_matrix(ls,ds) 
  a_of_schwinger_matrix_(ls,ds)
  a_of_schwinger_matrix(ls,ds) 
 
(%i1) linel:72
(%i2) s_kite:schwinger_matrix([l1,l2],
                              [l1^2+m2,l2^2+m2,(l1-l2+p)^2+m2,
                               (l1+p)^2+m2,((-l1)+l2)^2+m2])
ARRSTORE: use_fast_arrays=false; allocate a new property hash table for $ARRAY_A
               [ x5 + x4 + x3 + x1 ]         [ (- x5) - x3  ]
               [                   ]         [              ]
(%o2)  Col 1 = [    (- x5) - x3    ] Col 2 = [ x5 + x3 + x2 ]
               [                   ]         [              ]
               [    p x4 + p x3    ]         [    - p x3    ]
               [                      p x4 + p x3                      ]
               [                                                       ]
       Col 3 = [                        - p x3                         ]
               [                                                       ]
               [          2               2                            ]
               [ m2 x5 + p  x4 + m2 x4 + p  x3 + m2 x3 + m2 x2 + m2 x1 ]
(%i3) s_sunset:schwinger_matrix([l1,l2],
                                [l1^2+m2,l2^2+m2,(l1-l2+p)^2+m2])
          [ x3 + x1   - x3                p x3              ]
          [                                                 ]
(%o3)     [  - x3    x3 + x2             - p x3             ]
          [                                                 ]
          [                    2                            ]
          [  p x3    - p x3   p  x3 + m2 x3 + m2 x2 + m2 x1 ]
(%i4) s:subst([x4 = 0,x5 = 0],s_kite)
          [ x3 + x1   - x3                p x3              ]
          [                                                 ]
(%o4)     [  - x3    x3 + x2             - p x3             ]
          [                                                 ]
          [                    2                            ]
          [  p x3    - p x3   p  x3 + m2 x3 + m2 x2 + m2 x1 ]
(%i5) s-s_sunset
                              [ 0  0  0 ]
                              [         ]
(%o5)                         [ 0  0  0 ]
                              [         ]
                              [ 0  0  0 ]
(%o5)                       kiteVSsunset.mac