# NOT RUN {
u <- seq(0, 4, by=0.01)
z <- u^2/2 ### this is dfun
res <- compute_wt(z, cfun=1, s=1)
plot(z, res, ylab="Weight", type="l", lwd=2,
main=expression(paste("hcave", " (", sigma, "=1)", )))
# }
Run the code above in your browser using DataLab