Calculation Radioactive Isotope decay rate



#5

I hope someone can help me with this. I have programmers working on something for work, and I need to check their work. I know how to calculate the radioactive isotope decay rate, but I hoping someone can help me program it into my calculator so I can solve for faster. I have 48gx and 35s, 15c, and 28C that I can use. I prefer not use 28c because I have not fixed battery problem yet and I do not feel like losing program after I power it off.

I define formulas for you, but first let me define variables:

Ai = initial activity/radioactivity of isotope
t = elapsed time (hours)
A = current radioactivity (what we solve for)
H = half-life of radioisotope (hours)

Formulas:

A = Ai(1/2)^t/H

t = (H ln (A/Ai))/(ln (1/2)) <<- ln be natural logarithm

I sorry if this overly complicated, as I can do manually, but I would like to program and I am not sure how.


#6

Hello Vincze,

You can use the 48's solver. Press [Right Shift] [Enter] to activate the Equation writer. Enter the first equation; When you're done, press [Enter] to save it to the stack. Then save it as a function, typing a name, such as Activ:
['] [Alpha] [Alpha] ACTIV [STO]

Repeat the process for the second equation, and save it:
['] [Alpha] [Alpha] TIME [STO]

Now go to to the Solver menu, using [Right Shift] [SOLVE] You will see a menu: use the cursor keys to select SOLVE EQUATION and press [Enter] You will get the HP Solve screen, which will prompt you to select an equation. Look for ACTIV and press [Enter] The variables will be created automatically and you can start entering values. After you complete data entry, use the arrow keys to set the cursor in the variable you want to solve and press SOLVE.

The procedure is the same for the other equation.

Hope this helps.


#7

Juan, thank you so much my friend for the help. I know how to program solver but I wonder how it might be done with program to run. I still try and figure out program of calculator, and sometimes it confuse me. I guess solver is easier and for most case, do just fine, as I can solve for any variable, but it would be nice to see how to program the same problem.

Thank you again though.

#8

Vincze

One of the beauties of RPN and keystroke programming is that you program just as you would work out the problem. For Radioactive Decay, I prefer to use the integral form of first order reactions. All radioactive decay follows first order kinetics, which means that the following equations apply...

rate = k*[A]

upon integration you obtain

[A]t = [A]o Exp(-k*t)

where [A]o is the initial concentration
k is the rate constant which can be found by Ln (2)/half life
t is the time of reaction and
[A]t is the amount left after time t has passed by.

For the HP 28c and the 48gx the programming is the same. Using local variables "a" for [A]o, "b" for half life and "c" for time you can enter the program as follows.

<< -> a b c << 'a * EXP ( - ( LN(2) / b ) * c ) ' >> >>

save it as 'decay' then enter the values into the stack and press 'decay' from the user menu.

for this program the values have to be [A]o , half life and time in that order into the stack.

if you do not want to use local variables a shorter program would be

<< 2 LN SWAP / * NEG EXP * >>

save it as 'dacay' or any name you want to give it then enter [A]o, time and half life into the stack in that order. The program will return the amount left over. This last program should also work for the 35s, even though I do not have one, I believe you just need to change SWAP for x<>y and NEG for +/-.

Good luck

Sorry I had to edit to fix the formatting

Edited: 12 Aug 2007, 9:37 p.m.


Possibly Related Threads...
Thread Author Replies Views Last Post
  HP35s Program Four Slings Lift Calculation Jean-Marc Biram (Australia) 2 282 12-16-2013, 07:21 PM
Last Post: Jean-Marc Biram (Australia)
  [HP-Prime] AMBIGUITY between Numerical Calculation (HOME) and Numerical/Symbolic Calculation (CAS mode) CompSystems 2 187 08-18-2013, 07:06 PM
Last Post: CompSystems
  slow calculation Bill Plaxton 3 223 08-15-2013, 11:57 AM
Last Post: Bill Plaxton
  hp 17bii+ dates calculation... JoePaul 6 348 07-14-2013, 11:32 AM
Last Post: Katie Wasserman
  HP35 Ad. from 1972 with calculation error reused in 2009? Michael Kathke 46 1,708 05-14-2013, 03:52 PM
Last Post: Dieter
  Good (simple?) calculation for benchmarking? Jedidiah Smith 28 1,095 03-01-2013, 05:13 PM
Last Post: Harald
  [WP34s et al.] Solving the TVM equation for the interest rate Dieter 24 907 12-01-2012, 05:53 AM
Last Post: Paul Dale
  Complex number calculation on wp34s wildpig 13 438 08-30-2012, 09:46 AM
Last Post: Ángel Martin
  HP-10BII+ Price calculation bug? Thomas Radtke 0 121 02-15-2012, 06:53 AM
Last Post: Thomas Radtke
  HP 30b programs Interest Rate conversions Negative Periods per Year and Fractional Periods per year Gene Wright 1 116 05-03-2011, 11:25 AM
Last Post: Bart (UK)

Forum Jump: