gamma=[0:0.01:1];delta=sqrt(1-gamma.^2);p=(2*(abs(gamma.*delta)).^2)/(1+2*gamma.*delta*exp(-2)).^(-2);plot(gamma,p)