The R Project SVN R

Rev

Rev 35253 | Rev 42034 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 35253 Rev 42031
Line 29... Line 29...
29 
abs(Im(asin(sin(1i))) -	 1) < 2*Meps
 29 
abs(Im(asin(sin(1i))) -	 1) < 2*Meps
30 
cos(1i) == cos(-1i)# i.e. Im(acos(*)) gives + or - 1i:
 30 
cos(1i) == cos(-1i)# i.e. Im(acos(*)) gives + or - 1i:
31 
abs(abs(Im(acos(cos(1i)))) - 1) < 4*Meps
 31 
abs(abs(Im(acos(cos(1i)))) - 1) < 4*Meps
32 
 
 32 
 
33 
 
 33 
 
34 
.Random.seed <- c(0, 629, 6137, 22167) # want reproducible output
 34 
.Random.seed <- c(0L, 629L, 6137L, 22167L) # want reproducible output
35 
Isi <- Im(sin(asin(1i + rnorm(100))))
 35 
Isi <- Im(sin(asin(1i + rnorm(100))))
36 
all(abs(Isi-1) < 100* Meps)
 36 
all(abs(Isi-1) < 100* Meps)
37 
##P table(2*abs(Isi-1)	/ Meps)
 37 
##P table(2*abs(Isi-1)	/ Meps)
38 
Isi <- Im(cos(acos(1i + rnorm(100))))
 38 
Isi <- Im(cos(acos(1i + rnorm(100))))
39 
all(abs(Isi-1) < 100* Meps)
 39 
all(abs(Isi-1) < 100* Meps)