def prop_int_ksigma(NS,N,n,k):
    e=0
    S=SerieMoyennes(NS,N,n)
    mu=(1+n)/2
    sigma=EcartType(S)
    for m in S:
        if abs(m-mu)<k*sigma/sqrt(N):
            e=e+1
    return e/N