\begin{Schunk}
\begin{Sinput}
 library(scatterplot3d)
 s3d <- scatterplot3d(dat, type="p", color="blue", angle=55, pch=16, main="scatterplot3d")
 s3d$plane3d(m1)
\end{Sinput}
\end{Schunk}
