 %calcul.m
 
P=zgamma(20,m,s,c,a)-c/a
g=zgamma(20,m,s,c,a);
p=1-normcdf(g-c/a,m,s);
T=(1-p)/p Fa=g*(1+a)-c/a
