%% function [y,yp]=fr(snr)%% Defines rate function f and if necessary f' trough y and yp% as a function of a SNR vector.%function [y,yp]=fr(x)y = log(1+x);if nargout > 1    yp = 1./(1+x);end