・3次元データ
＞library(scatterplot3d)
＞x <- (1:9999)
＞y <- (1:9999)
＞xy <- x*y
＞scatterplot3d(x,y,xy)
a <- (1:99999)
b <- (1:99999)
c <- a^2/(a^2+b^2)
scatterplot3d(a,b,c)
和達三樹[1988],p115の検証
>#例題5.1（１）
> a2 <- (1:99999)
> b2 <- (1:99999)
> c2 <- (a^3+b^3)/(a^2+b^2)
> scatterplot3d(a2,b2,c2) ＃（0,0,0）で連続
>#例題5.1（2）
> a3 <- (1:9999)
> b3 <- (1:9999)
> c3 <- a3*b3/(a3^2+b3^2)
> scatterplot3d(a3,b3,c3)
>#（0,0）で連続とならない。
>#a3から原点に接近するときとb3からでは異なる値に近づく
> #問題5-1 2-3
> x1 <- (0:999)
> y1 <- (0:999)
> z1 <- sin(x1+y1)/(x1+y1)
> scatterplot3d(x1,y1,z1)
> #計算上：原点に近づけると1に近づく
> #問題5-1 2-3 99までの数字で見る
> x2 <- (0:99)
> y2 <- (0:99)
> z2 <- sin(x2+y2)/(x2+y2)
> scatterplot3d(x2,y2,z2)
> #問題5-1 2-2
> x3 <- (1:99999)
> y3 <- (1:99999)
> z3 <- (x3^2-y3^2)/(x3^2+y3^2)
> scatterplot3d(x3,y3,z3)
> #原点への接近。
> #y3軸だと-1、x3軸だと１にに近づく。
> #99までで確認
> x4 <- (1:99)
> y4 <- (1:99)
> z4 <- (x4^2-y4^2)/(x4^2+y4^2)
> scatterplot3d(x4,y4,z4)
