Displaying 1 result from an estimated 1 matches for "slv2".
Did you mean:
sslv2
2012 Mar 24
1
Solving the equation using uniroot
Hello all,
I was going to solve (n-m)! * (n-k)! = 0.5 *n! * (n-m-k)!
for m when values of n and k are provided
n1<-c(10,13,18,30,60,100,500) # values of n
kx<-seq(1,7,1); # values of k
slv2<-lapply(n1,function(n){
slv1<-lapply(kx,function(k){
lhs<-function(m)
{
fnk<-factorial(n-m)*factorial(n-k)-
0.5*factorial(n-m-k)*factorial(n);
return(fnk);
}
un2<-uniroot(lhs...