//
function erf(x) {
var fneg = 0;
if (x<=0) { fneg=1; x=Math.abs(x); }
var a1 = 0.0705230784;
var a2 = 0.0422820123;
var a3 = 0.0092705272;
var a4 = 0.0001520143;
var a5 = 0.0002765672;
var a6 = 0.0000430638;
var den = Math.pow(1 + a1*x + a2*Math.pow(x,2) + a3*Math.pow(x,3)
+ a4*Math.pow(x,4) + a5*Math.pow(x,5) + a6*Math.pow(x,6), 16);
var z = (1 - (1/den));
if (fneg==1){z = -z}
return z;
}
var u=100;
var x=70;
var s=0.15*u; // s = ecart type
var QI130Plus = 0.5 - 0.5 * erf((u-x)/(s*Math.sqrt(2)));
print ("P(QI>130) : " + QI130Plus);