Moi, propriétaire des droits d’auteur sur cette œuvre, la place dans le domaine public. Ceci s'applique dans le monde entier. Dans certains pays, ceci peut ne pas être possible ; dans ce cas : J’accorde à toute personne le droit d’utiliser cette œuvre dans n’importe quel but, sans aucune condition, sauf celles requises par la loi.

z <- (-1000:1000)/100
# Hubert k
k  <- 1.339
wh <- 1*(abs(z)<=k)+
      k*sign(z)/z*(abs(z)>k)

# Hampel
a   <- 1.7
b   <- 3.4
c   <- 8.5
wha <- 1*(abs(z)<=a)+
       a*sign(z)/z*((abs(z)>a) & (abs(z)<=b))+

       a*(c-abs(z))/(c-b)*sign(z)/z*((abs(z)>b) & (abs(z)<=c))
# Andrews wave
a   <- 1.34*pi
wav <- a/(pi*z)*sin(pi*z/a)*(abs(z)<=a)

# Tukey's biweight
a   <- 4.685
wtb <- (1-z^2/a^2)*(1-z^2/a^2)*(abs(z)<=a)

plot(z, wh, type="l", main="Hubert (k=1.339)")
plot(z, wha, type="l", main="Hampel (a=1.7, b=3.4, c=8.5)")
plot(z, wav, type="l", main="Andrews wave (a=1.34*pi)")
plot(z, wtb, type="l", main="Tukey's biweight (a=4.685)")

 Ce graphique a été créé avec R.