[С++ Builder 6] pow: domain error Опции

[18192123]

Здравствуйте!

Столкнулась вот с такой проблемой:
В программе генерируются случайные числа, распределённые по нормальному з-ну. Но суть не совсем в этом.
Когда требуется сгенерировать 500-1000 значений - всё нормально, когда же я захотела получить 10 000 значений - ошибка pow: DOMAIN error..

double xi = 0, xiStar = 0, yi = 0;
float sigma1 = 1, sigma3 = 1;
float m1 = 5, m3 = 3;
..............
double GaussLow(int var)
{
  randomize;
  xi = static_cast<double>(rand()) / static_cast<double>(RAND_MAX);
  xiStar = static_cast<double>(rand()) / static_cast<double>(RAND_MAX);
  switch (var)
  {
    case 1:
        yi = sigma1 * cos(2*M_PI*xi) * pow(-2*log10(xiStar),1/2) + m1;
        break;
    case 2:
        yi = sigma3 * cos(2*M_PI*xi) * pow(-2*log10(xiStar),1/2) + m3;
  }
  return yi;
}
............
void __fastcall TForm2MS::ButGenClick(TObject *Sender)
{
  double tempG = 0;
  int test = 10000;
  StrGrGen->RowCount = test+1;
  StrGrGen->Visible = true;
  StrGrGen->Cells[1][0] = "Íîðìàëüíûé ç-í";
  StrGrGen->Cells[2][0] = "Íîðìàëüíûé ç-í";
  for (int i=0;i<test;i++)
  {
    StrGrGen->Cells[0][i+1] = i+1;
    tempG = GaussLow(1);
    StrGrGen->Cells[1][i+1] = FloatToStrF(tempG,ffFixed,8,5);
    tempG = 0;
    tempG = GaussLow(2);
    StrGrGen->Cells[2][i+1] = FloatToStrF(tempG,ffFixed,8,5);
  }
}

Объясните пожалуйста, в чём тут дело? Как избежать ошибки..?
Подозреваю, что проблема можеть быть связана с областями определения (значения) pow или log10..

[volvo]

Может рандом выдаёт какое-то недопустимое значение аргумента...?

Вполне возможно... Я тебе даже скажу какое значение должно выдаваться, чтобы ты получала то, что получаешь: НОЛЬ... Если rand() вернул 0, то xiStar (естественно) тоже будет = 0, а это некорректные входные данные, это оговорено в описании функции. Так что придется тебе гарантировать, что аргумент log10 не будет нулевым. Ну, например, так:

int sure_rand()
{
	int rnd = rand();

	if(!rnd) rnd += 1;
	return rnd;
}
...
xiStar = static_cast<double>(sure_rand()) / static_cast<double>(RAND_MAX);

Если сделать просто:

xiStar = static_cast<double>(rand() + 1) / static_cast<double>(RAND_MAX);

, то получишь сбой, когда rand() вернет RAND_MAX... Тогда xiStar примет значение > 1, следовательно log10(xiStar) станет больше 0, и ты попытаешься вычислить квадратный корень из отрицательного числа. Поверь, ни к чему хорошему это не приведет...