前に作ったはずなのに見つからなかったドラゴンカーブのプログラムとほんのちょっとの違いで，結果はかなり違うものになるdrawC <- function(a.x, a.y, b.x, b.y, n){ x <- b.x-a.x y <- a.y-b.y c.x <- a.x+(x-y)/2 c.y <- b.y-(x-y)/2 if (n == 0) { lines…

Taglibro de Hモンテカルロ法http://ito-hi.blog.so-net.ne.jp/2006-09-05短くすればよいと言うものではないが，二次元配列を利用すると良い配列の二乗，列和に colSums，列和が 1 以下になる平均値（！）を求め（sum(p)/n って，mean(p) のこと）4倍するset…

これも以前書いたものdrawTree <- function(a.x, a.y, b.x, b.y, n, col=1){ STEM.RATIO <- 0.25 BRANCH.RATIO <- 0.6 xx <- b.x-a.x yy <- a.y-b.y angle1 <- atan(yy/xx)+pi/6 angle2 <- atan(yy/xx)-pi/6*1.5 center.length <- sqrt(xx^2+yy^2)*(1-STEM.R…

これも，以前書いたものdrawKoch <- function(a.x, a.y, b.x, b.y, n){ c.x <- (2*a.x+b.x)/3 c.y <- (2*a.y+b.y)/3 d.x <- (a.x+2*b.x)/3 d.y <- (a.y+2*b.y)/3 x <- b.x-a.x y <- a.y-b.y d <- sqrt(x^2+y^2)/sqrt(3) if (x >= 0) { a1 <- atan(y/x)+pi/6 …

Taglibro de Hドラゴン曲線http://ito-hi.blog.so-net.ne.jp/2008-09-06前に作っていたものだけどdrawDragon <- function(a.x, a.y, b.x, b.y, n){ x <- b.x-a.x y <- a.y-b.y c.x <- a.x+(x+y)/2 c.y <- b.y+(x+y)/2 if (n <= 0) { lines(c(a.x, c.x, b.x),…

Taglibro de Hhttp://ito-hi.blog.so-net.ne.jp/2009-11-05ついでにスピログラフのプログラムを while などを使わずに書いてみる。ずいぶん短くなるし，描画もあっという間であるけれど。while で描く方が，とろとろ描いていておもしろいのだけど（笑）。spi…

良くあるパターンであるが，引数を入れ替える必要はないTaglibro de Hhttp://ito-hi.blog.so-net.ne.jp/2009-11-05 ## 最大公約数 gcd <- function(a, b) { if (a < b) { t <- a a <- b b <- t } r <- a %% b return(ifelse(r == 0, b, gcd(b, r))) }a が b …

Taglibro de Hhttp://ito-hi.blog.so-net.ne.jp/2008-12-02?comment_success=2010-09-01T15:26:42&time=1283322402R: order()で並べ替え [統計] [編集]こういうデータがあったとする。m <- matrix(c( 261, "コナラ","Quercus serrata", 65, "ウリハダカエデ"…