plot Symmetric Information Rate for a channel in terms of p

Hi, I want to plot the symmetric information rate(SIR) for the channel based on the below expression:

SIR =(1/n)H(y1,y2,...yn)-h(p)/2

where H determines the entropy function and h(.) is the binary entropy function. note that input is an i.u.d process and z is the Bernoulli process with parameter 'p' and the relation between input,output and state z is as following:

I would be so pleased if some one helps me to do that. Thanks