Problema Simulación

[Linearz]

Bueno, el detale es este. Tengo que usar Montecarlo en una simulacion que genera numeros random distribuidos Uniforme(0,1), esto lo hago dividiendo el numero aleatorios mas grande de c++ (RAND_MAX) entre el rand().

EL problema en si es el output, no me evalúa bien las cosas. Les pongo lo que he hecho:

Código C++:

#include <iostream>
#include <cstdlib>
#include <ctime>
 
using namespace std;
 
double func_esperanza(double y)
{
    double f=0;
    if (y>200) {
        f=y*10-200*5-2*(y-200);                            //si hay mas demanda que 200 
    }
    else {
        f=y*10-200*5+.5*(200-y);                           //hay menos demanda que 200
    }
    return f;
}
 
double define_dato(double aux)
{
    double y;
    
    if ((aux>0) && (aux<=(45/300)))
    {
        y=100;
    }
    else if ((aux>(45/300)) && (aux<=(95/300)))
    {
        y=150;
    }
    else if ((aux>(95/300)) && aux<=((155/300)))
    {
        y=175;
    }
    else if ((aux>(155/300)) && (aux<=(215/300)))
    {
        y=200;
    }
    else if ((aux>(215/300)) && (aux<=(270/300)))
    {
        y=250;
    }
    else if ((aux>(270/300)) && (aux<=(1)))
    {
        y=300;
    }
    else
    {
        y=0;
    }
    
    return y;
    
}
 
 
 
double genera_dato()
{
    double R=RAND_MAX;                          //genero los aleatorios uniforme(0,1)
    double dato=rand()/R;                       //divido el mas grande aleatorio / rand()
    return dato;
}
 
 
 
int main()
{
    srand((unsigned)time(0));
    double aux,aux1,p=0;
    double suma=0;
    for (int i=0; i<10; i++) {
        aux= genera_dato();                             //genero el numero random
        cout<<aux<<endl;
        aux1= define_dato(aux);                         //lo evaluo en que intervalo cae
        cout<<aux1<<endl;
        p=func_esperanza(aux1);                         //calculo su esperanza
        cout<<p<<endl;
        suma+=p;
        aux=0;
        aux1=0;
        p=0;
    }
    
    cout<<suma/10;                                      //aplico montecarlo
    return 0;               
}
 

Este es el output que me ha salido:

0.0765
300
1800
0.736337
300
1800
0.623635
300
1800
0.43767
300
1800
0.914628
300
1800
0.153492
300
1800
0.737618
300
1800
0.149986
300
1800
0.80858
300
1800
0.797272
300
1800
1800

Cabe mencionar que uso Xcode(Mac, Mountain Lion). ¿Alguna sugerencia?

[Fw190]

Una cosa que veo que puede dar problemas son las comparaciones con aux que haces en la función define_dato. En la primera comparación:

if ((aux>0) && (aux<=(45/300)))

estás comparando aux (definida como double) con 95/300. Esta última expresión, al ser una operación entre constantes, se va a evaluar como división entera en tiempo de copmpilación, ya que ninguno de los operandos tiene punto decimal. Por tanto, 95/300 da cociente 0. Lo mismo ocurre en el resto de comparaciones que estás haciendo con aux en esa función. Prueba poniendo un punto decimal en, al menos, uno de los elementos (numerador o denominador) del cociente que comparas con aux, es decir, bastaría con poner, por ejemplo,

if ((aux>0) && (aux<=(45./300)))

con lo que ya obligas al compilador a evaluar ese cociente en coma flotante y no como entero.

Una sugerencia: En las constantes que vayas a utilizar en expresiones en coma flotante, pon siempre el punto decimal al final aunque sean valores enteros. Por ejemplo, en lugar de 5, pon 5. o, mejor aún, por legibilidad, pon 5.0

Prueba cambiando esto que te he comentado, y si sigue habiendo problemas ya buscaremos la causa del fallo.

Saludos,