def integral_s(K, sX):
    r = K-1
    f = symbols('f')
    s = symbols('s')
    const = 1/((2**(r/2))*gamma(r/2))
    f = (s**(r/2-1))*exp(-s/2)
    f = integrate(f,(s,sX,float('+inf')))
    return f*const

# проверка хи2
#array-эмпирический массив значений, K - кол-во интервалов
def hi2(array, K, N, lya_cur, tetha_cur, alpha):
  #  for i in range(len(array)):
  #      array[i] = (array[i]-tetha_cur)/lya_cur

    intervals = SelectIntervals(array, K)

    #эмпирическая частота попадания
    frecuency0 = np.zeros(len(intervals)-1)
    for item in (array):
        for i in range(len(frecuency0)):
            if(intervals[i] <= (item-tetha_cur)/lya_cur and (item-tetha_cur)/lya_cur <intervals[i+1]):
                frecuency0[i]+=1
                break

    # теоретическая частота попадания
    frecuency1 = np.zeros(len(intervals)-1)
    for i in range(len(frecuency1)):
        frecuency0[i] /= N
        if(intervals[i] < -1 and intervals[i+1] >= -1 and intervals[i+1] <= 1):
            frecuency1[i] = 0.1 - CosExpD(v, b1, intervals[i], 1, 0) + CosExpD(v, b1, intervals[i+1], 1, 0) - CosExpD(v, b1, -1, 1, 0)
        else:
            if(intervals[i+1] > 1 and intervals[i] >= -1 and intervals[i] <= 1):
                frecuency1[i] = CosExpD(v, b1, 1, 1, 0) - CosExpD(v, b1, intervals[i], 1, 0) - CosExpD(v, b1, intervals[i+1], 1, 0) + 0.1
            else:
                if(intervals[i+1] > 1 and intervals[i] > 1):
                    frecuency1[i] = - CosExpD(v, b1, intervals[i + 1], 1, 0) + CosExpD(v, b1, intervals[i], 1, 0)
                else:
                    frecuency1[i] = CosExpD(v, b1, intervals[i+1], 1, 0)- CosExpD(v, b1, intervals[i], 1, 0)

    result = 0
    for i in range(K):
        result += ((frecuency0[i] - frecuency1[i])**2)/frecuency1[i]
    result *= N

    print("Критерий Хи-Квадрат Пирсона")
    if (integral_s(K, result) > alpha):
        print(integral_s(K, result), '>', alpha)
        print("Гипотеза не отвергается")
    else:
        print(integral_s(K, result), '<', alpha)
        print("Гипотеза не принимается")

    DrawHistogram(intervals[0:K], frecuency0,intervals[0:K], frecuency1)
    return