Hello,(adsbygoogle = window.adsbygoogle || []).push({});

I am attempting to create a program that will evaluate the Bateman equation for radiactive decay series, given several decay constants stored in a text file. I am not sure where this is going wrong, but no matter what chain member I enter in as the desired nuclide to find the mass remaining after 1 million years, it outputs NAN - even when I choose the nuclide that I am certain to have a perceivable value. The equation I am emulating later in the program is shown below, please help to explain why this is not working! All help is appreciated!

http://img515.imageshack.us/img515/5097/screenshotml6.jpg [Broken]

Thanks again!Code (Text):

#include <iostream>

#include <math.h>

#include <fstream>

#include <string>

using namespace std;

int main(int argc, char** argv)

{

string line; //for file input

int iter; //stores which nuclide will be calculated for

double dcs[23]; //array for decay constants

int anum = 1; //counting var

ifstream deconst;

deconst.open ("decayconstants.txt", ios::in); //file with decay constants

if (deconst.is_open())

{

while (! deconst.eof())

{

getline(deconst, line);

dcs[anum]=strtod(line.c_str(), NULL); //store as double in decay constant array

anum++;

}

}

std::cout << "What nuclide in the chain should be evaluated?";

std::cin >> iter;

//I separated the equation into multiple parts

double consts=1; //part one: the multiplied decay constants at the beginning of the equation

double init = 100; //the intial mass, N(O)

double time = 1000000; //the variable t, one million years

double sum=0; // the second half of the equation, the summation

double totalmass; //the product of consts,init, and sum

for (int i = 1; i < iter; i++)

{

consts=consts*dcs[i]; //multiply all the necessary decay constants

}

for (int q = 1; q < iter; q++) //the summation

{

double expo = exp(-dcs[q]*time); //determine the numerator in the summation

double denom = 1;

for (int w = 1; w < iter; w++) //the product

{

if (w==iter) //k cant be j

{

break;

}

else

{

denom = denom * (dcs[w]-dcs[q]); //multiply together to get the product

}

}

sum = sum + (expo/denom); //add up all of the summations

}

totalmass = sum*consts*init;

std::cout << totalmass;

return 0;

}

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Bateman Equation in C++

**Physics Forums | Science Articles, Homework Help, Discussion**