library(e1071)

## density-estimation

     # create 2-dim. normal with rho=0:
     X <- data.frame (a=rnorm (1000), b=rnorm (1000))
     attach (X)

     # traditional way:
     m <- svm (X, gamma=0.1)

     # formula interface:
     m <- svm (~., data=X, gamma=0.1)
     # or:
     m <- svm (~a+b, gamma=0.1)

     # test:
     newdata <- data.frame(a=c(0,4), b=c(0,4))
     predict (m, newdata)

	par(mfcol=c(2,1))
     # visualization:
     plot (as.matrix(X))
     plot(as.matrix(m$SV), col=2)