2013-08-12から1日間の記事一覧

サインカーブの長さ

参考記事 L.sin <- function(a,b,x0=0,x1=2*pi,N=100000){ x <- seq(from=x0,to=x1,length=N) x <- x[-N] sum(sqrt(1+a^2*b^2*cos(b*x)^2)) * (x1-x0)/(N-1) } L.sin(1,1) as <- seq(from=0,to=20,length=50) Ls <- rep(0,length(as)) for(i in 1:length(as)…