Discussion :

[thasseda]

Bonjour tout le monde,

j'ai ce programme pour le calcule de la somme "s" et la somme carrée "v" d'une variable aléatoire x, lors de l'exécution pour N grand il me donne #IND ou #INF. qlq peut m'aider SVP

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
#include <cmath>
#include <cstdlib>
#include <iostream>
#include <fstream>
#include <math.h>
#include <string>
#include <ctime>
 
// monte carlo pour une variable aléatoire exponenetielle de paramètre lambda
 
using namespace std;
 
int main(int argc, char *argv[])
{
float x,U,s,v;
int N=100000;  
float lambda=2;
 
srand(time(NULL));                         
s=0; v=0;    // initialisation de la somme de la variable aléatoire (somme de x et x^2)
            for(int i=0;i<N;i++)
 
                {
                    U=(float)rand()/((float)RAND_MAX+1.0);
 
                     x= (-1/lambda)*log(U);
 
                     s=s+x; 
                     v=v+x*x;
                     cout<<"s="<<s<<endl;
                     cout<<"v="<<v<<endl;
 
 
                 } 
 
 
           system("pause");
return EXIT_SUCCESS;
}

[Altess]

Bonjour

Un variable int va de -32 768 à 32 767, or toi tu l'initialise à 100000. Ce qui est trop grand. Pour palier, tu peux utiliser un d'autres type, ou bien limiter ta valeurs N à quelque chose de raisonnable.

Je pense que ca vient de là, après, il se peut que je soit à coté de la plaque.

Cordialement,
PKO

[thasseda]

en fait N est le nombre pour lequel la boucle for tourne, donc pour une bonne précision de mes calculs, il faut que je mets un nombre très grand.

[screetch]

N va trés bien, les entiers vont jusque 2 milliards.
mais log(0) ca n'est pas défini

si U se rapproche de 0, log(U) tend vers -infini, et la variable flottante peut alors atteindre de très (trop) grandes valeurs. cette probabilité augmente lorsque N augmente.

[Altess]

N va trés bien, les entiers vont jusque 2 milliards.
mais log(0) ca n'est pas défini

si U se rapproche de 0, log(U) tend vers -infini, et la variable flottante peut alors atteindre de très (trop) grandes valeurs. cette probabilité augmente lorsque N augmente.

Mea culpa

[thasseda]

comment je fais alors SVP