tadpoles

   tp1 =
       + [d-4]^-1 * ( 16/9*mp^2*sh*[sh-mk^2]^-1 - 16/9*mp^2 + 52/9*mp^4*
         [sh-mk^2]^-1 + 68/9*mk^2 + 68/9*sh - 68/9*sh^2*[sh-mk^2]^-1 )

       + ln(mp)*[4Pi^2]^-1 * ( 3*mp^4*[sh-mk^2]^-1 )

       + ln(mk)*[4Pi^2]^-1 * (  - 2*mk^4*[sh-mk^2]^-1 )

       + ln(me)*[4Pi^2]^-1 * ( 8/9*mp^2*mk^2*[sh-mk^2]^-1 - 1/9*mp^4*
         [sh-mk^2]^-1 - 16/9*mk^4*[sh-mk^2]^-1 );

   tp2 = 0;

   tp3 = 0;

unitarity contr.

   uc1 =
       + [d-4]^-1 * ( 224/9*mp^2 - 119/9*mk^2 - 5*ss + 5/3*sh )

       + ln(mp)*[4Pi^2]^-1 * ( 12*mp^2 - 3/2*mk^2 - 9/2*ss + 3/2*sh )

       + ln(mk)*[4Pi^2]^-1 * (  - 7/2*mk^2 + 3/2*ss - 1/2*sh )

       + ln(me)*[4Pi^2]^-1 * ( 4/9*mp^2 - 29/18*mk^2 + 1/2*ss - 1/6*sh )

       + BB(ss,mp) * (  - 6*i_*mp^2 + 3/2*i_*mk^2 + 9/2*i_*ss - 3/2*i_*sh )

       + BB(ss,mk) * ( 3/2*i_*mk^2 - 3/2*i_*ss + 1/2*i_*sh )

       + BB(ss,me) * (  - 2/9*i_*mp^2 + 13/18*i_*mk^2 - 1/2*i_*ss + 1/6*i_*sh
          );

   uc2 = 0;

counterterms

   ct1 =
       + [d-4]^-1 * (  - 164/9*mp^2 + 59/9*mk^2 + 5*ss - 5/3*sh )

       + Nr(10) * ( 16*mp^2 )

       + Nr(11) * ( 40*mp^2 - 16*mk^2 )

       + Nr(21) * (  - 4*mk^2 + 4*ss - 4*sh )

       + Nr(22) * (  - 2*mk^2 + 2*ss + 2*sh );

   ct2 =
       + [d-4]^-1 * (  - 16/9*mp^2*sh*[sh-mk^2]^-1 - 44/9*mp^2 - 52/9*mp^4*
         [sh-mk^2]^-1 - 8/9*mk^2 - 68/9*sh + 68/9*sh^2*[sh-mk^2]^-1 )

       + Nr(10) * ( 16*mp^2*mk^2*[sh-mk^2]^-1 - 16*mk^4*[sh-mk^2]^-1 )

       + Nr(11) * ( 8*mp^2*mk^2*[sh-mk^2]^-1 + 8*mp^4*[sh-mk^2]^-1 - 16*mk^4*
         [sh-mk^2]^-1 )

       + Nr(21) * (  - 8*mp^2*sh*[sh-mk^2]^-1 + 8*mk^2*sh*[sh-mk^2]^-1 );

amplitude (sum)

   kss =
       + ln(mp)*[4Pi^2]^-1 * ( 12*mp^2 + 3*mp^4*[sh-mk^2]^-1 - 3/2*mk^2 - 9/2*
         ss + 3/2*sh )

       + ln(mk)*[4Pi^2]^-1 * (  - 7/2*mk^2 - 2*mk^4*[sh-mk^2]^-1 + 3/2*ss - 1/
         2*sh )

       + ln(me)*[4Pi^2]^-1 * ( 8/9*mp^2*mk^2*[sh-mk^2]^-1 + 4/9*mp^2 - 1/9*
         mp^4*[sh-mk^2]^-1 - 29/18*mk^2 - 16/9*mk^4*[sh-mk^2]^-1 + 1/2*ss - 1/
         6*sh )

       + Nr(10) * ( 16*mp^2*mk^2*[sh-mk^2]^-1 + 16*mp^2 - 16*mk^4*[sh-mk^2]^-1
          )

       + Nr(11) * ( 8*mp^2*mk^2*[sh-mk^2]^-1 + 40*mp^2 + 8*mp^4*[sh-mk^2]^-1
          - 16*mk^2 - 16*mk^4*[sh-mk^2]^-1 )

       + Nr(21) * (  - 8*mp^2*sh*[sh-mk^2]^-1 + 8*mk^2*sh*[sh-mk^2]^-1 - 4*
         mk^2 + 4*ss - 4*sh )

       + Nr(22) * (  - 2*mk^2 + 2*ss + 2*sh )

       + BB(ss,mp) * (  - 6*i_*mp^2 + 3/2*i_*mk^2 + 9/2*i_*ss - 3/2*i_*sh )

       + BB(ss,mk) * ( 3/2*i_*mk^2 - 3/2*i_*ss + 1/2*i_*sh )

       + BB(ss,me) * (  - 2/9*i_*mp^2 + 13/18*i_*mk^2 - 1/2*i_*ss + 1/6*i_*sh
          );


sum of all loops

   loops =
       + [d-4]^-1 * ( 16/9*mp^2*sh*[sh-mk^2]^-1 + 208/9*mp^2 + 52/9*mp^4*
         [sh-mk^2]^-1 - 17/3*mk^2 - 5*ss + 83/9*sh - 68/9*sh^2*[sh-mk^2]^-1 )

       + ln(mp)*[4Pi^2]^-1 * ( 12*mp^2 + 3*mp^4*[sh-mk^2]^-1 - 3/2*mk^2 - 9/2*
         ss + 3/2*sh )

       + ln(mk)*[4Pi^2]^-1 * (  - 7/2*mk^2 - 2*mk^4*[sh-mk^2]^-1 + 3/2*ss - 1/
         2*sh )

       + ln(me)*[4Pi^2]^-1 * ( 8/9*mp^2*mk^2*[sh-mk^2]^-1 + 4/9*mp^2 - 1/9*
         mp^4*[sh-mk^2]^-1 - 29/18*mk^2 - 16/9*mk^4*[sh-mk^2]^-1 + 1/2*ss - 1/
         6*sh )

       + BB(ss,mp) * (  - 6*i_*mp^2 + 3/2*i_*mk^2 + 9/2*i_*ss - 3/2*i_*sh )

       + BB(ss,mk) * ( 3/2*i_*mk^2 - 3/2*i_*ss + 1/2*i_*sh )

       + BB(ss,me) * (  - 2/9*i_*mp^2 + 13/18*i_*mk^2 - 1/2*i_*ss + 1/6*i_*sh
          );

sum of all counterterms

   wct =
       + [d-4]^-1 * (  - 16/9*mp^2*sh*[sh-mk^2]^-1 - 208/9*mp^2 - 52/9*mp^4*
         [sh-mk^2]^-1 + 17/3*mk^2 + 5*ss - 83/9*sh + 68/9*sh^2*[sh-mk^2]^-1 )

       + Nr(10) * ( 16*mp^2*mk^2*[sh-mk^2]^-1 + 16*mp^2 - 16*mk^4*[sh-mk^2]^-1
          )

       + Nr(11) * ( 8*mp^2*mk^2*[sh-mk^2]^-1 + 40*mp^2 + 8*mp^4*[sh-mk^2]^-1
          - 16*mk^2 - 16*mk^4*[sh-mk^2]^-1 )

       + Nr(21) * (  - 8*mp^2*sh*[sh-mk^2]^-1 + 8*mk^2*sh*[sh-mk^2]^-1 - 4*
         mk^2 + 4*ss - 4*sh )

       + Nr(22) * (  - 2*mk^2 + 2*ss + 2*sh );