import tensorflow as tf

    sigma -- Fermi-Dirac distribution

    """

    if activation == "esp" or activation == "softmax":
         A = tf.multiply(h, tf.nn.sigmoid(tf.multiply(beta,h)) )

    elif activation == "sigmoid":
         A = tf.nn.sigmoid(tf.multiply(beta,h))

