sigma function in solvers « Next Oldest | Next Newest »

 ▼ Don Shepherd Unregistered Posts: 1,392 Threads: 142 Joined: Jun 2007 03-03-2009, 12:36 PM Yesterday a college student asked me if it was possible to use the solver in the 17bii+ to implement the following equation: I created the following solver formula: ```ERLANG=(R^M/FACT(M)) / (1+SIGMA(I:1:M:1:R^I/FACT(I))) ``` Since many formulae require the summation, or sigma, function, I was wondering why it was not included in the 33s and 35s. Did the designers assume that, since those models allow keystroke programming (so you could "roll your own" sigma function), it just was not necessary? I would think it would have been nice to include the sigma function in the solver without requiring the user to learn keystroke programming. ▼ Peter A. Gebhardt Unregistered Posts: 174 Threads: 20 Joined: Sep 2006 03-03-2009, 04:18 PM Don, how do you solve for M (the # of servers in the ERLANG B formula, if given R "Erlangs")??? Just tried it on my 200LX - and get "Bad guesses". Best regards, Peter A. Gebhardt ▼ Don Shepherd Unregistered Posts: 1,392 Threads: 142 Joined: Jun 2007 03-03-2009, 07:15 PM Peter, the student that asked about this, he knows R and M, and he solves for ERLANG. I don't know enough about it to say more than that, but it seemed to be an easy formula for the 17bii+, given the sigma function. ▼ Peter A. Gebhardt Unregistered Posts: 174 Threads: 20 Joined: Sep 2006 03-04-2009, 06:47 AM Thx. Don. So I repost my above question to the other readers: What causes the HP200LX solver (and I guess the 17b/19b-family solvers too) to refuse to calculate M in Don's above implementation? Is this a problem of the solver 'per se' or a mathematical 'impossibility'? Best regards, Peter A. Gebhardt ▼ Marcus von Cube, Germany Unregistered Posts: 3,283 Threads: 104 Joined: Jul 2005 03-04-2009, 08:27 AM My 19Bii fails early: It does not seem to allow the FACT function in a solver equation. :( EDIT: That must have been a typing mistake, now the formula works. It looks like the problem does not have a solution for certain combinations of R end ERLANG. What does the formula compute and what would be reasonable input? Edited: 4 Mar 2009, 8:46 a.m. ▼ Peter A. Gebhardt Unregistered Posts: 174 Threads: 20 Joined: Sep 2006 03-04-2009, 09:44 AM Marcus, pls. take a look here: Best regards, Peter A. Gebhardt ▼ Marcus von Cube, Germany Unregistered Posts: 3,283 Threads: 104 Joined: Jul 2005 03-04-2009, 03:08 PM Thanks Peter! Now I know that the function is Erlang B, just an intermediate step in calculating Erlang C which is the probability that a customer has to wait for a server. But the article doesn't help with the proper values for Erlang B so that I can test the formula for solvability? :( Can somebody shed some more light on the subject? ▼ Don Shepherd Unregistered Posts: 1,392 Threads: 142 Joined: Jun 2007 03-04-2009, 09:05 PM Marcus, the student that asked me about this provided a copy of a page of his textbook which showed that, for r=0.67 and m=3, the probability is 0.0255 (from a lookup table in the textbook). My formula returned probability 0.0258. I can send that page to you if you want, because is does describe this function with a bit more detail. ▼ Marcus von Cube, Germany Unregistered Posts: 3,283 Threads: 104 Joined: Jul 2005 03-05-2009, 02:58 AM Just offer a scan for download here. I'm hopefully not the only one who is entertaining his brain with some "new" math. :) With the given values, the 19BII solves happily for M. If I change the value of ERLANG to the result given in the original paper, the solver approaches 3 but cannot reach the exact value for ERLANG, due to the fact that FACT returns an error for non integer arguments which makes the solver retry. I've modified the equation slightly: ```ERLANG=(R^M/FACT(IP(M))) / (1+SIGMA(I:1:M:1:R^I/FACT(I))) ``` Now I get 3,02707 for ERLANG. It's the continuous R^M part that makes this possible. Edited: 5 Mar 2009, 3:15 a.m. ▼ Don Shepherd Unregistered Posts: 1,392 Threads: 142 Joined: Jun 2007 03-05-2009, 04:12 AM Here is that little textbook writeup on this function, along with the lookup table. ▼ Marcus von Cube, Germany Unregistered Posts: 3,283 Threads: 104 Joined: Jul 2005 03-05-2009, 09:51 AM OK, I got it. Erlang B is the probability of loss (all servers busy) where R is the ratio of service time : average arrival rate and M is the number of servers. Thanks for the link. Peter A. Gebhardt Unregistered Posts: 174 Threads: 20 Joined: Sep 2006 03-05-2009, 05:11 AM Marcus, Thx. a lot! From your findings an advice given (independently of SIGMA) could be: "One should check the usage of variables that are supposed to stay integers during the whole calculation." Best regards, Peter A. Gebhardt Katie Wasserman Unregistered Posts: 1,477 Threads: 71 Joined: Jan 2005 03-03-2009, 06:11 PM Don, We explored the possibilities of the 35S solver here and concluded the it needs a G() function to solve equations like this but has a functional L() and can iterate. I don't know why they allowed for STO in equations but not RCL, it make no sense at all. -Katie Edited: 3 Mar 2009, 6:16 p.m. Allen Unregistered Posts: 562 Threads: 27 Joined: Oct 2006 03-03-2009, 08:29 PM Quote: Since many formulae require the summation, or sigma, function, I was wondering why it was not included in the 33s and 35s. Don, Greetings! Excellent post! Two quick observations: 1. The sigma function can be implemented in loops if not explicitly accessible. This is often faster since many of the Calcs that have a SIGMA function also incorporate single/double/x^2/y^2 etc into the Sigma+ key (e.g. 42s). I think most of the time these are not needed, so it's faster to run a summation loop to sufficient accuracy. 2. For many (I think not all..) summations, there is a general closed-form solution. (e.g the Change for a Dollar solution that Egan Provided) David Hayden Unregistered Posts: 528 Threads: 40 Joined: Dec 2008 03-03-2009, 08:47 PM The hard part about calculating the Erlang formula is that it's easy to exceed the precision of the calculator if you do the calculations in the wrong order. Dave Karl Schneider Unregistered Posts: 1,792 Threads: 62 Joined: Jan 2005 03-04-2009, 03:05 AM Hi, Don -- Quote: Since many formulae require the summation, or sigma, function, I was wondering why it was not included in the 33s and 35s. The crux of the matter is that the HP-33s/35s equation editor is fundamentally the one ported from the algebraic HP-22S into the RPN HP-32SII in 1991 and subsequently carried forward; it is not the more-advanced implementation of the HP-17B/27S. The Sigma operator for the HP-17B/27S equation editor will, by looping, compute finite series using an expression that defines every term. (Sigma+, by contrast, only adds a single datum to a summation.) The HP-17B/27S does not use summation for statistics, so there is no Sigma+ function. The HP-22S equation editor allowed only programmable calculating operations -- or most of them -- defined on its keyboard. Sigma+ is one that was not supported, because it is not very practical for an equation. Sigma as a construct used in equations, however, might have caused confusion with Sigma+. -- KS Edited: 4 Mar 2009, 3:19 a.m. ▼ Don Shepherd Unregistered Posts: 1,392 Threads: 142 Joined: Jun 2007 03-04-2009, 06:35 AM Thanks for the background, Karl. It makes me apprerciate the 17bii solver all the more. What a gem!

 Possibly Related Threads… Thread Author Replies Views Last Post HP50g: Writing a function that returns a function Chris de Castro 2 2,072 12-10-2013, 06:49 PM Last Post: Han [Sigma]DAT on the 50G Matt Agajanian 2 1,142 09-06-2013, 09:46 PM Last Post: Matt Agajanian [WP 34s] Summation Function (Sigma) Paul C 11 2,823 01-29-2013, 07:42 AM Last Post: C.Ret The Sigma Function on the WP34 Eddie W. Shore 7 1,985 04-26-2011, 07:47 AM Last Post: Namir Stupid question re HP-41 Sigma key... Jim L (OHIO) 5 1,681 07-19-2010, 06:04 PM Last Post: BruceH A Sigma Function in the 35s Solver! Katie Wasserman 16 3,872 08-21-2007, 10:10 AM Last Post: Gene Wright Solvers in 27S and 42S Anders 7 2,024 10-26-2006, 03:52 PM Last Post: Anders How to compare HP solvers ? jose goncalves (Brasil) 9 2,239 07-26-2005, 01:07 AM Last Post: Vieira, Luiz C. (Brazil) SIGMA(N:0:100:1:FACT(N))-X speed test on the 17BII+ Bob Wang 0 666 06-04-2005, 04:43 PM Last Post: bobwang EQ for algebraic solvers Rick Mave 8 1,973 12-30-2003, 09:09 PM Last Post: Katie

Forum Jump: