grttt nbbfg
2011-Dec-04 03:23 UTC
[R] upper bound in the integrate function depends on a parameter
Sorry for my English, is not my first language.. I have some trouble in terms of using integrate function in R. fx is a function of m and x where m is supposed to be a unknown parameter.>fx=function(m,x){+ x*2/(3*m)*(1-x/(3*m)) + } The problem is in upper bound, it depends on parameter m.>integrate(fx,lower=0,upper=3*m)$valueIs it possible to use the integrate function when bounds depend on a parameter? How to define m to solve this integral? The result should be m. [[alternative HTML version deleted]]
David Winsemius
2011-Dec-04 05:06 UTC
[R] upper bound in the integrate function depends on a parameter
On Dec 3, 2011, at 10:23 PM, grttt nbbfg wrote:> Sorry for my English, is not my first language.. > > I have some trouble in terms of using integrate function in R. > fx is a function of m and x where m is supposed to be a unknown > parameter.R is not an algebraic solver (and R-help is not a homework list.).> >> fx=function(m,x){ > + x*2/(3*m)*(1-x/(3*m)) > + } > > The problem is in upper bound, it depends on parameter m. > >> integrate(fx,lower=0,upper=3*m)$value > > Is it possible to use the integrate function when bounds depend on a > parameter?If it is known, yes. You could also do some investigation: curve(x*2/(3*5)*(1-x/(3*5)), 0, 15)> > How to define m to solve this integral? The result should be m.As I said, no homework answers here. -- David Winsemius, MD West Hartford, CT