diff --git a/data/akima760.rda b/data/akima760.rda
new file mode 100644
index 0000000000000000000000000000000000000000..4ea4ea7d29ff90c5ac585180f6927227c019efeb
Binary files /dev/null and b/data/akima760.rda differ
diff --git a/man/akima760.Rd b/man/akima760.Rd
new file mode 100644
index 0000000000000000000000000000000000000000..e8e14d2620045954ed007094779969273e06090b
--- /dev/null
+++ b/man/akima760.Rd
@@ -0,0 +1,63 @@
+\name{akima760}
+\alias{akima760}
+\title{
+Sample data from Akima's Bicubic Spline Interpolation code (TOMS 760)
+}
+\description{
+\code{akima760} is a list with vector components \code{x}, \code{y} and a matrix \code{z} which
+represents a smooth surface of \code{z} values at the points
+of a regular grid spanned by the vectors \code{x} and \code{y}.
+
+}
+\references{
+Hiroshi Akima, "
+
+", 
+ACM Transactions on Mathematical Software,
+Vol. 22, No. 3, September 1996, pp. 357-361.
+}
+\examples{
+\dontrun{
+library(rgl)
+data(akima)
+# data
+rgl.spheres(akima760$x,akima760$z , akima760$y,0.5,color="red")
+rgl.bbox()
+# bivariate linear interpolation
+# interp:
+akima.li <- interp(akima$x, akima$y, akima$z, 
+                   xo=seq(min(akima$x), max(akima$x), length = 100),
+                   yo=seq(min(akima$y), max(akima$y), length = 100))
+# interp surface:
+rgl.surface(akima.li$x,akima.li$y,akima.li$z,color="green",alpha=c(0.5))
+# interpp:
+akima.p <- interpp(akima$x, akima$y, akima$z,
+                    runif(200,min(akima$x),max(akima$x)),
+                    runif(200,min(akima$y),max(akima$y)))
+# interpp points:
+rgl.points(akima.p$x,akima.p$z , akima.p$y,size=4,color="yellow")
+
+# bivariate spline interpolation
+# data
+rgl.spheres(akima$x,akima$z , akima$y,0.5,color="red")
+rgl.bbox()
+# bivariate cubic spline interpolation
+# interp:
+akima.si <- interp(akima$x, akima$y, akima$z, 
+                   xo=seq(min(akima$x), max(akima$x), length = 100),
+                   yo=seq(min(akima$y), max(akima$y), length = 100),
+                   linear = FALSE, extrap = TRUE)
+# interp surface:
+rgl.surface(akima.si$x,akima.si$y,akima.si$z,color="green",alpha=c(0.5))
+# interpp:
+akima.sp <- interpp(akima$x, akima$y, akima$z,
+                    runif(200,min(akima$x),max(akima$x)),
+                    runif(200,min(akima$y),max(akima$y)),
+                   linear = FALSE, extrap = TRUE)
+# interpp points:
+rgl.points(akima.sp$x,akima.sp$z , akima.sp$y,size=4,color="yellow")
+
+}
+}
+\keyword{datasets}
+% Converted by Sd2Rd version 0.2-a3.