В коде эксперимента с квантовой информацией и безопасностью
e_l = (E_nu * Q_nu * exp(nu) - Y_0 / 2) / nu * Y_l;
есть ошибка, должно быть
e_l = (E_nu * Q_nu * exp(nu) - Y_0 / 2) / (nu * Y_l);