kayj <kjaja27 <at> yahoo.com> writes: > > > Hi All, > > I was wodering if it is possible to increase the precision using R. I ran > the script below in R and MAPLE and I got different results when k is large. > Any idea how to fix this problem? thanks for your help > > for (k in 0:2000){ > s=0 > for(i in 0:k){ > s=s+((-1)^i)*3456*(1+i*1/2000)^3000 > } > }
(1) see http://wiki.r-project.org/rwiki/doku.php?id=misc:r_accuracy:high_precision_arithmetic (2) consider whether there is more accurate algorithm you could use. I don't recognize the series, but perhaps it has a closed form solution, maybe as a special function? How much accuracy do you really need in the solution? Ben Bolker ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.