Message Boards Message Boards

0
|
2426 Views
|
3 Replies
|
3 Total Likes
View groups...
Share
Share this post:

Solve equation with Poisson Sum?

Posted 6 years ago

Consider the following code:

In[263]:= Clear["Global`*"]
            m = 2;
            a = 1/20;
            k = 0;

        Solve[Sum[PDF[PoissonDistribution[p], i], {i, k, m}] == a, p]

        During evaluation of In[263]:= Solve::nsmet: This system cannot be solved with the methods available to Solve.

    Out[267]= Solve[E^-p + E^-p p + 1/2 E^-p p^2 == 1/20, p]

I've tried many ways but haven't succeed. So how to solve this equation?

POSTED BY: Alex Graham
3 Replies

Sorry, I had misread the equation. Simply enlarge the interval:

Solve[Sum[PDF[PoissonDistribution[p], i], {i, k, m}] == a && 
  0 < p < 7, p]
POSTED BY: Gianluca Gorni
Posted 6 years ago

https://www.wolframalpha.com/input/?i=Solve%5BExp+(-p)+%2B+Exp+(-p)+p+%2B+Exp+(-p)+p%5E2%2F2+%3D%3D+0.05,+p%5D Internet version of mathematica finds the solution p=6.29 as it should be. But it writes that this solution is over reals, but 6.29 is real number, not comlex. So my question is how to find this 6.29 solution?

POSTED BY: Alex Graham

Your function has negative derivative and it vanishes at +infinity. There is no real solution. You can find complex solutions if you give a bound for p:

Solve[Sum[PDF[PoissonDistribution[p], i], {i, k, m}] == a && 
  Abs[p] < 2, p]
POSTED BY: Gianluca Gorni
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract