search for: quadritur

...gral(jac*finv(s2)*L(X|finv(s2))*p(finv(s2))))/(integral(L(X|finv(s2))*p(finv(s2)))) . the bounds on both integrals is (-1,1) jac = the jacobian and finv is the inverse function of f(s2) once again the integrals are calculated using the intergrate function Method 3 ========= legendre gaussian quadriture is used to calculate each of the integrals using Method 2. Some results and comments are below. method 3 seems to give the correct solution for n<=343 differences between the different methods can be detected from as early as n=3 for method 1 IS THERE A WAY TO MODIFY THE PROBLEM SO THAT...

hello R users strange behavior of the integrate function! i assume this occurs because of the way in which the quadriture is set up! (any comments.) f=function(x){exp(-exp(-x)-5*x)/gamma(5)} xx=seq(from=-20, to=20, length.out=1000) plot(xx,f(xx),type="l") integrate(f, lower=-Inf, upper= 1) integrate(f, lower=-Inf, upper= 10) integrate(f, lower=-10, upper= 11) integrate(f, lower=-Inf, upper= 11) integra...