%infa.m

for i=1:6
    a(i)=i*0.001;
end la=length(a); m=200000;s=0.3*m;c=250;

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

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