%infc.m

for i=1:200
    c(i)=10*i;
end lc=length(c); m=200000;s=0.3*m;a=0.001;

Fa=zeros(1,lc); for i=1:lc
    Fa(i)=zgamma(200,m,s,c(i),a)*(1+a)-c(i)/a;
end

T=zeros(1,lc); for i=1:lc
    g(i)=zgamma(20,m,s,c(i),a);
    p(i)=1/2*erfc((g(i)-m-c(i)/a)/sqrt(2*s^2));
    T(i)=(1-p(i))/p(i);
end
