Иногда функция возвращает значение NaN

[Lost7070]

Доброго времени суток.
У меня следующая проблема. Не всегда, но в 70% случаев функция возвращает значение NaN. Никак не могу найти причину этого. Очень буду признателен, если подскажете, чем это может быть вызвано и где стоит искать косяк.
Ниже код самой функции:

double distributions::gaverage(double **matr,double number_v,double number)
{
    double *average = new double [number_v];
    double temp=0;
    for(int j=0; j<number_v; j++)
    {
        for(int i=0; i<number; i++)
            temp=temp+matr[j][i];
        average[j] = temp/number;
        temp=0;
    }
 
    double p_s=0;
    temp=0;
    for(int i=0; i<number_v; i++)
        temp=temp+average[i];
    p_s=(1/number)*temp;
    
    temp=0;
    for(int i=0; i<number_v; i++)
        temp=temp+(pow((average[i]-p_s),2));
    distributions::GAverage=(1/number)*temp;
    return distributions::GAverage;
}

Если недостаточно какой-либо информмации, скажите, я предоставлю.
Заранее благодарен всем отозвавшимся.

[s-andriano]

Естественно, недостаточно.
При каких именно входных значениях функция возвращает NaN?

[Lost7070]

Ну, к примеру при значениях: number_v=100, number=600 и двумерный массив случайных чисел, числа генерятся следующим образом:

double distributions::random(double a, double b)
{
return ((b-a)*((double)rand()/RAND_MAX))+a;
}

double distributions::nr(double p1, double p2)
{
return ((sqrt(-2*log(random(0,1))))*(cos(2*3.14*ra ndom(0,1)))*p1)+p2;
}
p1=1, p2=0.

При этих и любых других возвратиться может и не NaN, но в большинстве случаев он.