#include "kin_funcs.h"
#include <math.h>


// This is from Byukling Kayanti Formula (6.3)
double Lambda( double x, double y, double z )
{
  return (x - y - z)*(x - y - z) - 4*y*z;
}

//From Byukling Kayanti Formula (5.14)
double T_min( double ma, double mb, double m1, double m2, double s)
{
  return ma*ma + m1*m1 - (1/(2*s))*( (s + ma*ma - mb*mb)*(s + m1*m1 - m2*m2) - sqrt( Lambda(s, ma*ma, mb*mb)*Lambda(s, m1*m1, m2*m2) ) );
}

//From Byukling Kayanti Formula (5.14)
double T_max( double ma, double mb, double m1, double m2, double s)
{
  return ma*ma + m1*m1 - (1/(2*s))*( (s + ma*ma - mb*mb)*(s + m1*m1 - m2*m2) + sqrt( Lambda(s, ma*ma, mb*mb)*Lambda(s, m1*m1, m2*m2) ) );
}