gauss_der = (
        tf.exp(-((t - t_final / 2) ** 2) / (2 * sigma**2))
        * (t - t_final / 2)
        / sigma**2