Any suggestion plz. What could be missing that I donot get function root to be zero at some points?
On Tue, 24 Dec 2019, 1:45 PM Shaami, <nzsh...@gmail.com> wrote: > Dear Friends > > Could you please help with the following example? I do not get function > value at zero using uniroot() in R. Could you suggest something? I want > function value zero at any point like 2.4, 2.5, 2.6, 4.5, 5.9 etc. > > Thank You > > x = c(-100, -20, 100, 230, 500, 920, 1000) > a = c(1000, 12, 300, 0, 100, 300, 750) > Fn = function(t, x, a){ > sum = 0 > for(i in 1:length(x)){ > sum = sum + punif(t, min=x[i]-a[i], max=x[i]+a[i]) > } > return(sum) > }uniroot(function(t){Fn(t, x, a) - 2.4}, interval=c(min(x-a), max(x+a))) > uniroot(function(t){Fn(t, x, a) - 2.5}, interval=c(min(x-a), max(x+a))) > uniroot(function(t){Fn(t, x, a) - 2.8}, interval=c(min(x-a), max(x+a))) > >
_______________________________________________ Rcpp-devel mailing list Rcpp-devel@lists.r-forge.r-project.org https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel