### R code from vignette source 'contfrac.Rnw' ################################################### ### code chunk number 1: contfrac.Rnw:132-134 ################################################### library("contfrac") as_cf(pi, n=7) # calculate the first 7 terms ################################################### ### code chunk number 2: contfrac.Rnw:142-144 ################################################### (jj <- convergents(as_cf(pi,n=7))) jj$A/jj$B - pi ################################################### ### code chunk number 3: contfrac.Rnw:154-157 ################################################### (phi <- (1+sqrt(5))/2) as_cf(phi,n=50) rle(as_cf(phi,n=50)) ################################################### ### code chunk number 4: contfrac.Rnw:179-182 ################################################### tan_cf <- function(z,n=20){GCF(c(z,rep(-z^2,n-1)),seq(from=1,by=2,len=n))} z <- 1+1i tan(z) - tan_cf(z) # should be small