Piyush3dB

function [C r] = BlahutArimoto(p)

disp('BlahutArimoto')

% Capacity of discrete memoryless channel
% Blahut-Arimoto algorithm

% Input
% p: m x n matrix
% p is the transition matrix for a channel with m inputs and n outputs