Find the capacity of a channel between X and Y by maximizing the mutual information

The gradient based approache uses Newton's method to find a critical point of the Lagrangian

where probabilities are represented on log-scale.