f=5.11;
 m=200.15;
 a=m*f/(f-1);b=m-a;
 sigma=102.66/1.41;
 lambda=m*f/((f-1)*sigma);
 d=lambda^2;

u=0;y=m;v1=lambda*(y-m)/sqrt(a*(y-b));v2=-lambda*(y+a-b)/sqrt(a*(y-b))


g1=.5*(1+erf(v1/sqrt(2)));
g2=.5*(1+erf(v2/sqrt(2)));
frep=g1+(exp(2*d))*g2